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ABSTRACT 


Several approaches for an attitude control system are 
studied and developed for a large space construction 
base tha-t is structurally flexible. Digital simulations 
were .obtained using the following techniques: (a) the 
multivariable Nyquist array methqd combined with closed 
loop pole allocation, (b) the linear quadratic regulator 
method. Equations for the three-axis simulation using 
the multilevel control method were generated and are 
presented. Several alternate control approaches are also 
described. A technique is demonstrated, .for obtaining the 
dynamic structural properties of a vehicle, which is con- 
structed' of two or more submodule's of known dynamic 
characteristics. 
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SECTION 1 


1 . 0 INTRODUCTION 

This report is submitted in compliance with the statement 
of Work under contract NAS8-32660. The period of per- 
formance covered by the report is the calendar year of 
August 26, 1978 to August 26, 1979. The submission and 
approval of this report constitute the successful com- 
pletion of the Exhibit "A" and Exhibit "B" portions of 
the contract. 

This report is the second of a series under the same 
contract number. The previous report (Reference 1-1) 
was submitted in October 1978 and covered the period 
from July 27, 1977 to July 27, 1978 in compliance with 
Exhibit "A" of contract NAS8-32660 . 

1 . 1 OBJECTIVES 

The sections that follow summarize the effort expended 
on the Space Construction Base (SCB) Control System 
Study contract. The primary objective of the study 
was to develop a control system and flexible control tech- 
niques that will stabilize a large and growing space 
station of the future. The study emphasized control 
methods, but additional topics include dynamic modeling 
of the vehicle and a technique for combining the dynamics 
of two or more submodules into one overall model . 

1 . 2 SCOPE 

Study effort was concentrated in two major areas: 

a. Application of different methods for control of 
large flexible spacecraft . 


1-1 



b. Generation of flexible vehicle dynamic models. 

Control equations were developed for three distinct 
approaches : 

a.. Multilevel control 
'b. Linear quadratic regulator (LQR) 

c. Multivariable Nyquist array (MNA) method and Pole 
Placement Techniques. 

Single axis simulations were generated for control of 
the SCB Configuration 1 using the LQR and MNA techniques. 
Control equations were generated for three-axis control 
of the same vehicle using the Multilevel conti*ol method. 

A heavy programming load resulting from one of the 
vehicle modeling tasks pre-empted the completion of 
this three-axis multilevel control simulation. 

Equations were generated for this three-axis vehicle 
model- and are presented .al ong with control equations for 
the multilevel approach. Single axis vehicle models 
also were derived for application with -the other control 
approach simulations. One additional vehicle modeling 
task was a structural analysis study in which the dynamic 
models of two submodules were combined to form a single 
structural model. Vibration modes for .the individual sub- 
modules and those for the composite structure were 
obtained using computerized matrix techniques for finite 
elements of the bodies. Analytical techniques also were 
applied to verify the results. 
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Some assumptions used throughout the report are that 
• sensor and actuator dynamics are negligible in the control 
analysis because of their relatively broad bandwidths, 
computer quantization and other non-linearities were assumed 
to be negligible for the time being. 

A' number of other items related to, SCB control are not 
discussed here since they .were completed in the previous 
study (Reference 1—1): 

a. . Mission requirements 

b. Pointing and orientation requirements 

c. Actuator sizing requirements 

d. Momentum desaturation requirements-. 

1 . 3 GENERAL 

This report is comprised of nine sections, Section 
2 contains an SCB vehicle model (Configuration 1), which 
is used in Section 4 in equation form and in Section 6 as 
numerical data. Sections 3 through 7 treat control 
' techniques for attitude stabilization of flexible vehicles. 
Section 3 contains a preliminary description of -each of 
the control techniques utilized in Sections 4, 5 and 6. 
Section 4 describes equations for three axis simulation 
of the multilevel control method. Section 5 presents a 
single axis simulation of control using the linear 
quadratic regulator (LQR) method. Section 6 describes 
a simulation of one axis of the SCB using the multi- 
variable • Nyquist array (MNA) method of control, and 
also applies a technique for closed loop pole placement. 
Section 7 discusses several alternate control approaches 
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which may be considered in future efforts. Section 8 
treats a potentially important tool for flexible vehicle 
modeling: a structural analysis study for combining the 

dynamics of two or more distributed parameter submodules 
into a single model when the submodules are securely 
fastened together. Section 9 is a discussion of con- 
clusions and recommendations. A number of Appendices 
are placed at the end of the report. References are 
cited liberally and are listed at the end of each 
section . 

The original RFQ requested that the International System 
of units (designated as SI) be used in the program and 
in any reporting. Torques, moments, angular momentum, 
moments of inertia and distances, however, are stated in 
English units since this was the system of units used 
in presenting all of the vehicle inertia data in the RFQ. 

A liberal approach was also followed in the numbering of 
SCB configurations. The original RFQ initially used 
Roman numerals, but the text here used Arabic numerals 
in most cases — except where material is reprinted from 
the RFQ. 

1 . 4 REFERENCES 

1-1 Guidance Systems Division, The Bendix Corporation, 
Space Construction Base Control System, Final Report, 
Contract NAS8-32660 for Geroge C. Marshall Space Flight 
Center, Oct. 27, 1979. 
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SECTION 2 


2.0 DYNAMIC MODELING OE THE VEHICLE 

Configuration 1 of the SCB consists of two solar wings 
connected to four contiguous modules which are considered 
to be rigid relative to the solar wings. Each solar 
wing is converted to two rigid bodies as shown in 
Figure 2-1. The five individual bodies are designated 
as in the diagram of Figure 2-2.- They are connected 
to each other via the three degree of freedom hinge 
suspensions, designated as 2 through 5 as shown. 

2.1 MATHEMATICAL MODEL 

The mathematical model for the SCB Configuration 
1 spacecraft was obtained using the discrete, coordinate 
method described in Reference 2-1. The system model 
consists of five rigid bodies connected by joints or 
hinges and subjected to arbitrary external forces and 
torques. The modeling procedure assumes that the 
bodies are connected in the form of a topological tree 
and hence contains no "closed loops". The relative 
motion about each joint is assumed to be rotational, and a 
three axis model which includes all cross-coupled axis 
interactions was developed. 
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(The individual bodies are connected to each other 
via three degree of freedom spring hinge suspensions.. 
Hinge numbers are circled.) 

TOPOLOGICAL DIAGRAM OF' SCB CONFIGURATION 1 


FIGURE 2-2 





The linearized model, derived based on small relative 
assumptions using the discrete coordinate approach 
(References 2-2 and 2-3), is expressible in the form. 


where : 


defined 


Pw = Q 


( 2 - 1 ) 


P is 15 x 15 nonsingular matrix, function of the 
model parameters (inertias, masses, etc.) 

w is 15 x 1 angular velocity vector, defined as 
w 


= [o)i X , w ov , w~,,, w— • w 


2X 

U 4X’ UJ 4Y’ '"4Z’ ^5X 
Q is 15 x 1 .torque vector. 


(0 


J 1Y ’ "lZ- 
W 4Y’ W 4Z * 


2Y’ 2Z’ 3X ’ 

T 


3Y ; 


3Z 


“ 5 Y> w 5z1 

i 


as : 


Q - [q 2 , Q 2 > %> %, Q 5 ] 


and 


T- 


Q i“ (Q ix> Q iy’ Q iz 5 for 1 
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The vector components of the Matrix Q are given as follows 


Ql = + T o - T s - T 


a. 


S 


12 


13 


q =T + T 4 - T 0 -v T 

2 eg a 2 S 12- 


^3 = T e + T a + T S * T S 
3 e -3 3 13 35 


( 2 - 2 ) 


Q4 = T e + T a + T S 

4 e 4 a 4 ^24 


Q 


T + T + T q 
e 5 a 5 b 35 


where T_ - is an external torque vector on body i 
e i given by T = (T , T ,T ) T 
■ e i e iX e iY e iZ 


T is an actuator, torque vector on body i 
• a i given as T q = (T„ T a .„> T a ,„) T 


a- 


a, 


iX a iY a iZ 


To is a suspension torque vector between 
13 bodies i and j given as 


To . - (To , T , T )' 

S ij S ijX S ijY S ijZ 
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Define thirty state variables as follows: 


X 1 - 

W 1X 



CO 

■H 

X 

= *i 


X 2 ‘ 

“lY 

(2-3) 


X 17 

- e i 

(2-8) 

X 3 “ 

“lZ 



X 18 



X 4 “ 

W 2X " 

to 

IX 

X 19 

" ^2 

- 


X 5 - 

CO — 

2Y 

“lY (2-4) 

X 20 

= 0 2 

- 6 1 

(2-9) 

X 6 “ 

co - 

2Z 

to 

1Z 

X 21 

= h 

" *1 



X 7 = 

W 3X “ 

w ix 



X 22 


^3 

- 4> x 


X 8 = 

W 3Y 

W 1Y 

(2- 

-5) 

X 23 


°3 

“ 9 1 

(2-10) 

X 9 ‘ 

W 3Z “ 

W 1Z 



X 24 

=5 

^3 

- 


X 10 

= “4X 

W 2X 



X 25 


*4 

4> 2 


X 11 

“4Y 

“ W 2Y 

(2- 

-6) 

X 26 

- 

0 4 

“ 0 2 

(2-11) 

X 12 

W 4Z 

" W 2Z 



X 27 


^4 

- V 


X 13 

W 5X 

~ “3X 



X 28 


^5 

<P 3 


X 14 

u 5Y 

W 3Y 

(2- 

-7) 

X 29 

SS 

0 5 

- 0 3 

(2-12) 

X 15 

“5Z 

~ W 3Z 



• X 30 

SS 

^5 

~ ^3 



Where (Jk, 9^ i]k are the Euler angles about X, Y and 
Z axis, respectively, for body i. 


d> . . = <j> . - d> . represents a relative angle of body i 
ij i J 

with respect to body j in the X direction* 


Similarly , 


G ij = 0 i"V “iJ x Wix 


w ijy w fy w iy 


u 

and to. 


ijz “iz “jz * 
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Also define: 


X Hl ^ X l’ X 2 ’ X 3 ’ X 4’ X 5 J X 6’ X 7 ’ X 8’ X 9’ X 10 5 
X ll» X 12’ X 13’ X 14 ’ X 15 ^ 2 “ 13 ? 


H2 


X 24’ X 25’ X 26 > X 27 1 X 28’ X 29’ X 30^ 


'21 ’ 

X 22 ’ X 23 

29 ’ 

y 

dU ( 2 - 

X 6’ 

Xy , Xg , 

‘13 J 

X 14 ’ X 15 


X 16 ’ X 17 ’ X 18 ’ X 19 ’ X 20 ’ X 21 ' X 22 ' 
X 23’ X 24 * X 25 ' X 26 ’ X 27’ X 28’ X 29’ 


X 30>‘ 


(2-15) 


Prom the definition of the state variables: 


X H2 X H1 


(2-16) 


i.e., X 16 = X 15 X 1? = X 2 , etc. 

Equations (2-3) through (2-7) can easily be rearranged 
so that to can be expressed in terms of X H1 by the 
following relationship: 

a) = N X R1 (2-17) 

Where N is 15 x 15 constant, nonsingular matrix. 


Differentiating (2-17) with respect to time yields; 


co = N X T 


(2-18) 
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Multiplying both sides of equation (2-1) by P 1 and 
substituting for angular acceleration vector u from 
equation (2-18) yields; 

N X„. = P -1 Q (2-19) 

HI 

Multiplying both sides of equation (2-19) by N" 1 gives; 

= N" 1 P _1 Q (2-20) 

HI 

Assuming no external torques, equation set (2-2) can be 
written as; 

Q = T a + T g (2-21) 


where T is 15 x 1 actuating torque vector. 

£L 

To is 15 x 1 suspension torque vector, given by; 


T Q = (-To - To , To - To 

S S 12 S 13 S 12 S 24 

To - To , To , T ) T 

S 13 S 35 S 24 S 35 


( 2 - 22 ) 


(-To - T„ ), etc., are 3x1 vectors. 
b 12 b 13 

The scalar suspension torques can be expressed in terms 
of rate and displacement constants and inertial rate 
difference and relative displacements, for example; 


To 

12X 

= - Ko 

b 12X 

^21 

- Co 

°12X 

(w 

2X 

T 

S 12Y 

# 

= - Ko 

* 12Y 

0 21 

• 

C S - 
b 12Y 

( W 2Y 

• 

• 

• 

• 

• 

• 

T g - 

°35Z 

“ - Kg 

°35Z 

^53 

C o 

S 35Z 

(“5Z 


0) ) 
IX 


w-j^y ) (2—23 ) 

) 
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Substituting equation set (2-23)for all scalar suspension 
torques in (2-22), using definitions of the state vari- 
ables (2-3 through 2-12) and re-arranging the terms, 
equation (2-22) can be expressed as 

T s “ KD X HL + KS *H2 (2 ' 24) 

where KD is 15 x 15 matrix of damping constants and KS 
is 15 x 15 stiffness matrix. 


Substituting (2-24), into (2.-21) gives 

Q = T a + KD X^ + KS X.jj2 (2-25) 

substituting Q given by (2-25) into (2-20) and 
multiplying by N 1 yields; 

X hl = N “ 1 P ' 1 K + “ X H1 + KS X H2 1 (2 - 28) 

Rewriting (2-16) in the following form: 


X 


H2 


= I 


X H1 + 0 X H2 


(2-27) 


where I is 15x15 identity matrix 
and 0 is 15x15 null matrix. 
Combining (2-26) and (2-27) Yields 


r * 

^1 

‘ X H2 " 


N 1 P 1 KD 
I 


N 1 P 1 KS 
0 



1 

1 — 1 

1 

J 

*P 


+ 


N 1 P 11 


T 


a 


0 


(2-28) 
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This can be further expressed as: 


X(t) = AX(t) + Bt ( t ) 


where X = £x 1} X 3( £j 

X ( t ) = T a = [ T alx> T aly) T alZj 


(2-29) 


T 

’ a5z 


f 

J 15x1 


r -i -l 
1 N P KD 


A = 


N 1 P 1 KS 


0 


30x30 


B = 


r N 1 P“ 1 


0 


30x15 


Assuming no torque applied to body 4 and 5 , t can be writter 
as 


-[ 


lx, 


,T 


3Z. 


;] 


9x1 

and B can be reduced to 30x9 matrix. 


(2-30) 


Numerical values of matrices P, N, KD, KS, A and B 
used for analysis and design purposes are shown in 
Appendix A. 
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2.2 Model Reduction Eigenvalue Analysis 

The eigenvalues for the five-body three-axis models were 
obtained as follows: (Set A) 


(1) 

-.0149042 

+ 

3 

1.61213 

(2) 

-.0268672 

+ 

3 

1.20486 

(3) 

-.0153772 

+ 

3 

0.642112 

(4) 

-.00228948 

+ 

3 

0.312132 

(5) 

-.43191 

+ 

J20.9923 . 

(6)- 

-.462163 

+ 

j20 . 4238 

(7) 

-.0653141 

+ 

3 

7.13316 

(8) 

-.0356952 

“I- 

3 

5.55713 

(9) 

-.00502089 

+ 

3 

1.00192 

(10) 

-.00500057 

4- 

3 

1.00005 

(ID 

-.00172613 

+ 

3 

0.322123 

(12) 

-.00160014 

+ 

3 

0.32009 


There are also six eigenvalues at the origin. 

Setting all cro'ss-axis coupling terms to zero, three 
five-body single-axis models can be developed. Under 
this condition the eigenvalues were computed as follows 
(Set B): 
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X-axis 


• Y-axis 


0.0 0.0 

0.0 0.0 


-.014904 

+ 

3 

1.61213 

-.431645 

+ 

j 20 . 989 

-.0266921 

+ 

3 

1.18467 

-.462423 

+ 

J20.427 

-.0153771 

+ 

3 

0.642111 

-.0654125 

+ 

3 7.13772 

-.00222516 

+ 

3 

0.30449 

-.0356052 

+ 

3 5.55064 


Z-axis 

0.0 

0.0 

-.00502089 + j 1.00192 
-.00500057 + j 1.00005 
-.00172613 + j 0.322132 
-.00160014 + j 0.32009 

A comparison of the single axis eigenvalues with 

the three axis eigenvalues leads to the conclusion 
that, from an eigenvalue assignment viewpoint, the three 
axis vehicle control system can be subdivided into three 
single axis control systems. Further comparison 
suggests that cross-axis coupling from an eigenvalue 
perspective is minimal, and thus an increase in the 
damping and/or stiffness factors of the system by 
active control can be accomplished by single axis 
measures . 

To examine the eigenstructure further, consider the 
single axis representations of Figure 2-2. The 
eigenvalues at the origin represent the rigid body 
dynamics of the main unit while the flexible modes are 
attributed to the solar wings. If the main unit were 
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sufficiently larger than the solar wing, the dynamics 
would not interact. The eigenvalues would then appear 
in double complex pairs. As the mass of the main unit 
decreases the eigenvalues of the open loop single axis 
system will shift in a predetermined way. These 
eigenvalues’ will then appear -as single complex pairs. 

Assume the state model for a single axis to be represented 
-by the form 

x = Ax + Bx (2-31) 


with the state vector as 

x(t ) = [ w 2 1^4 2^21^42 1 U 31 W 53^ 31^53 : W l^l] 


(2-32 ) 


.T 


or in partitioned form 
x ( t ) = [x 1 (t) :x 2 (t) :x 3 (t)] 

Substituting (2-33) into (2-31) yields 


! H* 


A 11 

A 12 

A 13 

X 2 

= 

A 21 

A 22 

A 23 

_*3_ 


_ A 31 

A 32 

A 33. 


x-, 


x. 


X, 


+ 


(2-33) 


B-, 


B, 


L B 3 


(2-34) 


with A 11 = A 22 a 4x4 matrix representing individual 
solar wing dynamics, A^ 2 = A 2 ^ a 4x4 matrix representing 
the dynamic interaction of solar wing A on solar wing B, 
A 13 = 0, A 23 = 0, A 31 and A on are 2x4 matrices reflecting 


32 
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the interaction of the solar wings on the rigid body and 



contains the rigid body dynamics. Thus (2-34) becomes 


A 11 

A 12 

0 



A 12 

A 11 

0 

x + Bt 

(2-35) 

A 31 

A 32 

A 33- 




The open loop eigenvalues of (2-35) consist of the 
eigenvalues of A g3 (designated as A(A 33 )) plus the 
eigenvalues of the partitioned block 


"A 

A 

_ 


11 

A 12 

12 

A ll. 


(2-36) 


With A^ and A^ square and of the same dimension 
(i.e., 4x4), row and column manipulations on (2-36) yield 



0 




L A 12 A 11 + A 12 J 


(2-37) 


Since eigenvalues of (2-36) are unaffected by row and 
column manipulations, the eigenvalues in (2-36) equal 
the eigenvalues of (2-37). But (2-37) is in block 
diagonal form, thus the eigenvalues of A in (2-34) equal 
the eigenvalues of A 33 , plus the eigenvalues of A^ - A^g. 
plus the eigenvalues of A^ + A^J i.e., 
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(2-38) 


A(A) = MA 33 ), A ( A h +A 12^ X(A 11~ A 12 ) 


If A^^ = 0, then the eigenvalues will appear in double 
complex pairs. Hence any departure of A^g from the zero 
matrix will cause an eigenvalue shift in the manner 
described by (2-38). 


Re-ordering the state vector for each of the single axis 
models as in (2-32), the single axis eigenvalues pre- 
viously computed were similarly obtained using (2-38). 

Thus the significant eigenvalue shifts in both the X and Y 
axes result from the cross-coupled interactions of the 
solar wings in the respective axis and not by a-ny cross- 
axis coupling. 


To obtain a qualitative measure of the interactive 

eigensh'ifts, assume ] l> 2 ± ~ ^31 anc * ^42 = ^53 an koth‘ 
magnitude and phase for all time. This is equivalent 
to a symmetry assumption about the main unit. The 
eigenvalues for each axis are subsequently computed as 
(Set C) : 


X-axis Y-axis 



0.0 



0.0 



o 

o 



o 

o 


-.0153606 

+ 

3 

0.641803 

-.0697786 + 

j 7.13515 

-.0149224 

+ 

3 

1.61229 

-.482064 + 

j 20. 9874 

-.0153606 

+ 

3 

0.641803 

-.0697786 + 

j 7.13515 

-.0149224 

+ 

3 

1.61229 

-.482064 + 

j20 . 9874 
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Z-axis 



0.0 


| 


0.0 



-.00500564 

+ 

3 

1.00112 

-.00160836 

+ 

3 

0.320495 

-.00500564 

+ 

3 

1.00112 

- .00160836 

+ 

3 

0.320495 


Comparing eigenvalue Set C with eigenvalue Set B we 
conclude that the symmetry assumption for the z axis 
(and possible the y axis) is worthy of further investiga- 
tion. The eigenvalue shift in the x-axis as caused by 
the solar panel interations suggests that this axis 
model should be treated in its entirety. 
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3 . 0 CONTROL APPROACHES 

Three diverse techniques are, in the succeeding sub- 
sections, described which are applicable to the control 
of large, flexible vehicles. The first technique 
(multilevel control) solves a complex control problem 
by reducing the solution to a large number of simpler 
problems. The second approach (LQR) utilizes some of 
the basic concepts of modern control theory for determin- 
ing the optimal control loop feedback. The third scheme 
(MNA) reduces the effective cross-coupling of a complex 
"plant" to an extent that the closed loop design reduces 
to a set of independent single loop design problems; 
techniques for obtaining closed loop pole placement are 
also utilized in this approach. 

The theoretical background and references for each control 
approach are presented in this section. The actual 
equations and some simulation results as applied to 
Configuration 1 of the Space Construction Base are 
offered in Sections 4, 5 and 6 for the multilevel, 

LQR and MNA control approaches respectively. 

3 . 1 MULTILEVEL CONTROL 

The control of large-scale multivariable systems is a 
class of control problems that arises in many areas of 
practical application. Direct application of control to 
a large-scale system often severely taxes or exceeds 
available computer capacity due to the high dimensions 
involved in the overall control problem to be solved. 
Application of multilevel systems analysis techniques 
to the large overall control problem decomposes it into 
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a multilevel hierarchy of subproblems of smaller dimension. 
An example of a two level subproblem hierarchy is 
depicted in Figure 3-1. 


In general a subproblem hierarchy resulting from the 
application of multilevel control techniques may contain 
several levels of subproblems. The roles of the sub- 
problems are correlated with the levels that they occupy 
in the hierarchy. Each of the subproblems on the lowest 
(first) level of the hierarchy' pertains to a portion of 
the original system to be controlled. Each of the sub- 
problems on a level above the lowest one pertain to 
coordination of the solutions of the subproblems on the 
next lower level . The number of subproblems per level 
decreases as the level becomes higher in the hierarchy 
until the top or supremal level is occupied by a single 
overall coordination subproblem. 

Each subproblem comprising the hierarchy transmits its 
solution (s) to other subproblems in the hierarchy in the 
course of the solution of the overall control problem. 
Such a solution proceeds iteratively. For example, 
the subproblems of the two level hierarchy may be 
solved in the following sequence. The second level 
coordination subproblem provides coordination variables, 
C, £, which are held fixed while the optimal solutions 
for x> A. and u are determined for all of the first 
level subproblems independently. After the optimal 
solution is obtained for each first level subproblem, 
numerical responses are sent to the coordination sub- 
problem on the second level. The coordination subproblem 
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COORDINATION EQUATIONS 


fM | 

r~ 

c 


X 


STATE 

u 

CONTROL 

X 

COSTATE 

EQUATIONS 

-<-= 

EQUATIONS 

"< ~ ■ 

EQUATIONS 


J 


FIGURE 3-1 Two Level Subproblem Hierarchy 
For Multilevel Control 
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now adjusts the coordination variables and sends their 
adjusted values to the subproblems on the first level. 

The first level subproblems are again solved independently 
with the adjusted values of the coordination variables. 
This procedure continues iteratively until the overall 
coordination problem at the top of the hierarchy is 
satisfied. The result sought is- an optimized collection 
of subsystems, with interconnections- restored , which is 
equivalent to the original controlled system optimized. 
This procedure may be readily extended to subproblem 
hierarchies with more than two levels. 

3.1.1 General Approach 

The overall procedure to be utilized in applying multi- 
level control techniques to a mathematical model of a 
system to be controlled is outlined as follows: 

1. Express mathematical model in state variable form. 

2. Decompose mathematical model into set of decoupled 
equations . 

3. Construct decomposed performance index. 

4. Form decomposed Hamiltonian. 

5. Develop costate equations with associated coordination 
equations . 

6. Develop control algorithm. 
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7. Construct subproblem hierarchy. 


8. Simulate subproblem hierarchy on digital computer. 
3.1.2 State Variable Model 


x = Ax+Bu 


(3-1) 

where:' 



x = (x 1jX2j x n ) T 

(state vector) 

(3-2) 

u = C U 1, X 2, V T 

(control vector) 

(3-3) 

A = nxn matrix 


(3-4) 

B = nxm matrix; m<n 


(3-5) 


Equation (3-1) is the canonical state variable form of 
the mathematical model of the system to be controlled. 

It represents a set of n simultaneous scalar ordinary 
linear differential equations which; in general, are 
coupled . 

Typical scalar state equation: 

n m 

x. = f.(x,u) = E a . .x. +• E b. .u. i = 1, 2, ,n (3-6) 

1 1 j=l U 3 J=1 3 

where f - ,l n ) T (3-7) 

3.1.3 Decomposed Model 

Define state coordination variables with following 
coordination equations. 
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(3-8) 


C. 

1 


Sa - 


n 

E 


m 


a . . x . + E b . . u . 
j=l 1J 3 j=l 3 3 


m 




Then the scalar state equations may be written in 
decomposed (decoupled) form as follows: 


x. - f. = a. .x.+b. .u.+C, 

1 x XI 1 11 1 1 


with the initial condition, 


(3-9) 


V x ± (o) = X lo 


(3-10) 


3.1.4 Performance Index 

Original performance index to be minimized: 


■ / 


( x^Qx+u^Ru ) dt 


(3-11) 


o 


where 


Q = nxn weighting matrix 

R = mxm weighting matrix 

t = initial time in interval considered 
o 

t f = final time in interval considered 

If the matrices Q and R are- diagonal, the performance 
index is said to be saparable and it may be written in 
decomposed (decoupled) form as follows: 
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P = EP. 
1 


i = 1,2; , n ■■ 


' (3-12) 



t 

If 



where: P. = 

l 

/ p * dt 
t 

o 


(3-13) 

2 2' 

Pi = + r .u. 

for x = 

1,2, m 

(3-14) 

2 

Pi = q^ 

for i = 

m+1 , , n 

(3-15) 

Hamiltonian 

The Hamiltonian is proportional to the energy 

of the 

system and may be 

written 

in decomposed form 

as follows 

H = EH. 

l 

i = 

1,2, , n 

(3-16) 


where: 


H i = Pi +X i f i + P 1 (g i -C i ) (3-17) 

= ith scalar costate variable 

p i = ith Lagrange multiplier associated with the ith 
state coordination equation. 

According to Pontryagin ' s minimum principle minimization 
of the Hamiltonian corresponds to optimization of the 
controlled system with respect to the performance index 
given in (3-11). Application of the necessary conditions 
for minimization of the Hamiltonian of the controlled 
system generates the costate, control and costate 
coordination equations for the system to be optimized. 
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Costate Equations 


X . 
1 



n 


-a . 


11 


A.-2q.x. - E 
1 1 1 3=1 

m 


a ji P 3 


1,2, ,n 


(3-18) 


with the final condition, 


A i (t f ) = 0 (3-19) 

Costate - Coordination Equations 

|§ = 0 -v p ± - A.; i = 1,2, , n (3-20) 

i 

Control Equations 

For a system that is operating suboptimally with respect 
to its performance index, some of the necessary optimality 
conditions applied to its Hamiltonian will not be 
satisfied at the outset. The selection of the specific 
set of necessary conditions to be temporarily violated 
at the beginning of operation of the system determines 
the nature of the decomposition of the system. In' 
this case it is assumed that the system decomposition 
is such that all of the necessary conditions listed so 
far, (3-18), (3-20), and the necessary conditions leading 
to the state and state coordination equations are 
satisfied throughout the time interval of interest, 
(t 0 ,t f ). Given these conditions, for suboptimal opera- 
tion at the outset, the necessary condition associated 
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with the control variables cannot be satisfied at the 
outset. It is therefore approached during the operation 
of the system by means of the following gradient 
formulation . 


(U i- l )l+l 


- <-Va - s i 


(M ) 

<8u iV 


i = 1 , 2 ,- — , in 


( 3 - 21 ) 


where: 

% = the iteration subscript 

s^ = constants to be chosen on the basis of the rate 
of approach of the controlled system to optimal 
operation. 


From (3-8), (3-9), (3-14) and (3-17) an expression 
may be written for the gradient of the Hamiltonian with 
respect to each scalar control variable as follows. 



2r . u . + 


x 1 


n 

l 

3=1 


b . .1 . 
J 


(3-22) 


3.1.9 Subproblem Hierarchy 

Equations (3-8) and (3-20) may be assembled into an 
overall coordination subproblem for the generation . of 
the optimal contrcl contours. Equations (3-9)' with their 
associated initial conditions, equations (3-10), constitute 
the first level state subproblems to be solved. Since 
these equations are temporarily decoupled from each other, 
they may be regarded as n subproblems each of which 
can be solved independently. Similarly, (3-18) with the 
final conditions of (3-19) may be associated into a first 
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level costate subproblem and (3-21) and (3-22) may 
be assembled into an overall control subproblem. These 
four subproblems may then be assembled into a hierarchy 
with the coordination subproblem at its apex as shown 
in Figure 3-1. 

3.2 LINEAR QUADRATIC REGULATOR (LQR) 

With a linear system described in the state variable 

form 

x = Ax + Bu (3-23) 

where x = the system state variable vector 
u = the control vector 
A,B = the system matrices, 

an optimum, feedback control u = Fx can be determined by 
utilizing Pontryagin’s minimum principle. 

An expression for the optimum feedback F can be derived 
(References 3-4 and 3-5) through application of either 
Hamilton-Jacobi equations or utilizing the second method 
of Liapunov with a "quadratic cost" system performance 
index . 

T 

PI - J (x T Qx + u T Ru) dt (3-24) 

o 

where T = upper limit in time (could also be infinite in. 
some derivations) 
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Q,R = weighting factors, where Q must be at least a positive 
semidefinite matrix and R must b.e positive definite 

For an initial state x(0), a matrix Riccati equation is 
generated by either approach from (3-23) and (3-24), 

-P..- PA + A T P - PBR -1 B T P + Q; (3-25) 

where P = P(t) and must be positive definite 
P = P(t)* 

i 

and the equation is subject to the boundary condition 
P(T) = 0. Ordinarily, this equation can be integrated 

backward in time from the known terminal condition over 

# 

the time period of interest, P approaches- zero as T is 
increased towards infinity, and the P matrix becomes a 
constant designated as P. The optimum feedback matrix is then 

F = -R -1 B T P (.3-26) 

as derived from the approaches -cited. Equation (3-23) 
can now be rewritten as a function of the state variables 
only 


x = (A-BF)x = (A+BR 1 B T P)x (3-27) 

Note that this equation implies that all the state 
variables must be fed back to obtain an optimal control 
system. Careful selection of state variables, or 
filtering of a set of measured variables, will be required 
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for feedback control of large flexible vehicles. These 
structures are actually a network of distributed spring 
and mass parameters, rather than a finite number of 
discrete or lumped elements. 

An application of the LQR technique to the control of a 
flexible vehicle, represented by discrete parameters, 
is demonstrated in Section 5 of this report. A solution 
of the system's matrix Riccati equation is given, along 
with an optimal closed loop simulation. 

3.3 MULTIVARIABLE NYQUIST ARRAY (MNA) 

The MNA control method used here is combined with a 
pole placement technique. The basic MNA approach will be 
described first and will then be followed by a brief 
description of the pole placement procedure. 

3.3.1 MNA Design Procedure 

The fundamental objective of the MNA design method is to 
decrease crosscoupled system interaction to such an 
extent that the closed loop system design reduces to a 
set of independent single loop design problems. Although 
simply stated, the actual reduction procedure proposed by 
Rosenbrock (Reference 3-6) and implemented by Munro 
(Reference 3-7) requires a high degree of designer 
intervention and is fundamentally a trial and error 
process . 

In Figure 3-2, G(S) is an mxm transfer matrix represent- * 
ing the coupling of m inputs and m outputs. Restriction of 
G(s) to a square matrix simplifies the analysis, programming 
and computational effort considerably. 
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FIGURE 3-2 MULTIVARIABLE SYSTEM 
CONFIGURATION . 


'The pre- and post— compensator matrices K- and- L, respectively 
are each of dimension mxm. The feedback gain matrix, F, 
is assumed -to be .diagonal and of similar dimensions. 

Clear Ty, if 

Q(S) = LG( S)K ' (3-28) 

is diagonal, loop closure may proceed on an individual 
loop basis with a guarantee of zero loop interaction. 

It is this premise upon which the MNA design philosophy 
is based/ The adherence to strict diagonalization is 
relaxed, however, with the substitution and exploitation 
of the concept of diagonal dominant matrices . 

A matrix Z(S) is diagonal dominant if either or both 
of the following conditions are present for all S: 

m 

a. Z | z. .(S) | / |z. . (S) |<0 -<1 for all i=l,2, m (3-29) 

• -* JL . I 4. -L JL 

3=1 
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m 

b. E |z.,(S)|/Iz. i (S)|<e.<l for all i = 1,2, m (3-30) 

j=l 3 

3fi 

Equation (3-29) defines row dominance while (3-30) defines 
the column dominance condition where 9^ is the level of 
dominance obtained for the ith diagonal element. 

Before the design process can proceed further, Q(S) for 
the direct Nyquist array (DNA) method or Q ^S) in the 
inverse Nyquist array (INA) method must be made dominant 
by manipulation of the elements of the compensator matrices. 
Once dominance is achieved the design process is completed 
using single loop theory to select the diagonal elements 
of E. This selection process is enhanced through applica- 
tion of the Gershgorin and the Ostrowski theorems for 
dominant matrices. 

The Gershgorin theorem (Reference 3-6) states that the 
eigenvalues of a matrix (either Q(S) or Q 3 (S))are 
located in the union of the bands centered about the 
diagonal elements with widths determined by the sum of 
the moduli of the off diagonal elements by row or by 
column. Using the envelope procedure developed by 
Crossley (Reference 3-8) and considering each control 
loop separately, a graphical display of open loop system 
interaction results. The feedback gain selection for 
control loop i is then made in correspondence with the 
generalized Nyquist criterion and the stability theorems 
of Rosenbrock (Reference 3-6). 
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The Ostrowski theorem may be used to further shrink 
the Gershgorin bands, thereby reducing the area of un- 
certainty in each loop. This set of bands is frequently 
referred to as a set of "fuzzy" Nyquist plots (or inverse 
Nyquist plots for the INA). Using the innermost band 
as a conservative estimate of the Nyquist contour in 
each loop, the. design proceeds on a single loop basis. 
Feedback gain selection must be made exterior to. the 
Gershgorin (Ostrowski) band. Thus, phase margin, gain margin an 
dynamic compensation may be used to evaluate and/or 
improve the loop design with a guarantee of low interaction 
from the closure of the remaining loops . 

Diagonal dominance for the DNA method requires the 
selection of pre- and post-compensator matrix para- 
meters so that (3-29) is satisfied when (3-28) represents 
the open loop transfer matrix. For the FNA method Q ^(S) 
is used and the parameters of K and L ^ must be selected. 

An efficient and reliable method for the evaluation of the 
matrix coefficients is described in Reference 3-9. .The 
dominance algorithm uses a conjugate direction function 
minimization algorithm to adjust the parameter set until 
a performance index composed of the dominance definitions 
in (3-29) and (3-30) is minimized. For the INA method 
in a row dominance mode, the optimization problem can be 
separated into three independent optimization efforts; 
one for each row. Here the' performance index by row is 

J, (K. .) = Max E iq.. i (S)l/U..(S)! i=l,2, ,m (3-31) 

J w j=l 

i*i 


3-15 



where §. .(S) is an element of Q 1 (S). For each i, the ■ 
ith row of K is adjusted until J^K^) is minimized. 

In practice, the ratio in (3-31) is computed for each 
discrete frequency point in the range of interest . This 
array is then scanned to identify the maximum ratio. ■ 

Adjusting the elements of row i in K ^ yields a set of 
final dominance levels 

0. = Min J .. (K . . ) (3-32) 

1 K. . 1J 

If the dominance levels in (3-32) are less than unity, 
diagonal dominance has been achieved. In the event that 
some of the dominance levels are greater than one, the 
designer may initiate a dominance sharing search or 
restart the program using new starting values for the 
unspecified compensator parameters. 

The concept of dominance sharing is detailed in Reference 
3-10. It is fundamentally a rescaling of the compensator 
matrices to the extent that low dominance levels may be 
intentionally increased to a point where the previous 
non-dominant levels may be shifted to a range of acceptability. 
This procedure has been automated in the latest version 
of the dominance algorithm and is initiated .by the 
designer after the final set of dominance levels have 
been evaluated. 

Once a set of coefficients has been determined for 
dominance of the open loop transfer matrix, each control 
loop may be treated independently using single-input/ 
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•single-output control theory.. The dominance algorithm 
briefly outlined above shifts the- burden of establishing 
the dominance condition from the designer to an automated 
procedure. Thus designer, intervention, is only required 
during the actual design process and is no longer 
needed to' establish the requisite dominance condition. 
Experience with the algorithm suggests that dominance 
may often be obtained within- several CPU minutes or less 
depending upon the characteristics of G(S), 

Using the appropriate transfer function's for the vehicle 
model, a DNA- design was initiated over a frequency 
range which includes the critical modes.. Recent results 
described in reference 3-11. for finite frequency Nyquist 
array design established this range .as acceptable for 
design and stability purposes-. For this DNA design the 
post compensator matrix was prespecified as the identity 
matrix with the precompensator matrix selected in' . 
accordance with the dominance algorithm in reference 3-9. 
The closed loop diagram for vehicle control., with 
Q(S) = L G(S)K 

representing the open loop transfer matrix, is shown 



Figure 3-3 Closed Loop Spacecraft Control Design 
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3.3.2 


Retallack Pole Placement Method 

The dyadic pole placement technique developed by Simon 
and Mitter, reference 3-12, and further refined by 
Retallack and MacFarlane, reference 3-13, utilizes the 
model structure of a linear system to give the designer 
direct control over the movement of given system modes 
to desired closed loop locations. This method is par- 
ticularly attractive from a ’.design viewpoint whenever 
the open loop poles can be directly associated with 
specific parts of the physical system. In this way 
only those poles directly associated with the subsystem 
will be shifted. Let Aj (j=l,....B) represent the open 
loop poles of the system and let p^(i=l, . . . . ,n) be the 
desired closed loop set. Several algorithms are pre- 
sently available to shift the eigenvalue set {A.} to 

J 

the desired set {p^}. However there is no guarantee 
that any particular eigenvalue A k will be moved to a 
specific location p^. Without the eigenvalue location 
assignment flexibility, the resulting control in light 
of sensor and/or actuator failures could result in very 
unsatisfactory performance. 

Consider the linear, time-invariant, multivariable 
system described by 

x'=Ax+Bu (3-34) 

where x is n vector of states, u is m vector of inputs 
and m<_n. If constant gain state feedback of the form 

u = y -Kx (3-35) 

is introduced into the system of (3-34) the resulting 


3-18 



system obeys 


X = (A-BK)x '+ By (3-36) 

where y is the external control . The open-loop poles 

■(A 1 ,....X n ) of the system of (3-34) and the closed 

loop poles (P 1 »-.-P n ) of the system (3-36) are related 

by the Hsu-Chen equation (reference 3-14). 

n 

. 7T ’ (S-p . ) 

IV< S V A > r— 


IT 

J=1 




The Jordan canonical form for A is given by 

A=U 1 AU - (3-38) 

where U is a modal matrix for A: Let ( A A . . . .A ) p<m 

.be the open loop poles to be shift.ed and define the 
modal controller K, in terms of the unity rank dyadic 
product 


K=fd T V (3-39) 

■where f is an m vector d is ! a p vector and V contains 
the p rows of U. corresponding to (A^ . A ) .* With this 
definition } K Retallack (reference 3-13) reduces equation 
(3-37) to 

P 


T - -i i'V S -Pi) 

1 + d 1 (SI -A) V 

P ’ p 

TT (S-A.) 

j=l 3 


(3-40) 


where 

3 = VBf • (3-41) 

and A is the px p Jordan submatrix for (A, ,A 0 ,...A ). 

-L d p 

The remaining open-loop modes (A , _ A ) are unchanged 

>p+l n 
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under dyadic feedback. 


To overcome some of the numerical difficulties associated 
with the computation of the modal controller when clus- 
tered poles appear near the imaginary axis (reference 
3-15) the system in (3-34) is transferred to the Luenber- 
ger canonical form (reference 3-16), 

w = AW + Bu (3-42) 

where w=U ^x 

A=U“ 1 AU (3-43) 

B=U _1 B (3-44) 

An efficient algorithm to realize this transformation is 
presented in reference 3-17. Since the eigenvalues of 
U-fiK) are identical to those of( A-BK) , the pole place- 
ment process can be accomplished in the closed loop 
system of (3-42), rather than in the closed loop system 

A 

of (3-34). Once the gain matrix K has been determined 
for the canonical system, the desired modal controller 
K for the original system (3-34), is given by 

K = KU _1 (3-45) 

This pole placement method was used in this study to 
shift the open loop poles of the vehicle model for 
analysis and control system design purposes. 
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SECTION* 4 


4.0 SIMULATION OF MULTILEVEL CONTROL 

The principal objective of this section is the develop- 
ment of mathematical models of a flexible space vehicle 
that are amenable to application of multilevel control 
techniques and construction of digital computer simula- 
tions of each vehicle model and its associated multi- 
level control system. In this development effort 
mathematical models for attitude stabilization of the 
Space Construction Base (SCB) have been utilized because 
these models are deemed representative of l!he corresponding 
models of flexible space vehicles i-n general and because 
data were available at the outset for the SCB models. 

The effort expended on the simulation of models of the 
vehicle with multilevel control since the last report, 
Chichester (4-1), has been divided into three main 
areas which are discussed in .subsections 4.1, 4.2 and 
4.3, respectively. Subsection 4.1 describes the 
development of a linearized scalar state variable form 
of the three dimensional discrete mass rotational model 
for Configuration 1 of the SCB from the model originally 
presented by Cornell (4-2). This subsection, summarizes 
work presented in Lipski (4-3) and (4-4), and Chichester 
(4-5). The resulting models are suitable for direct 
application of multilevel control techniques. 

Subsection 4.2 presents the application of multilevel . 
control techniques to the models of subsection 4.1. 

This application follows the sequence previously des- 
cribed in Chichester (4-1), (4-6). The resulting model 
is of rather high dimension for all but the minimum . 
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form of the SCB, Configuration 1, and the most limited 
. application of control torques, application to the central 
body only . 

Subsection 4.3 describes simple hybrid coordinate models 
of Configuration 1 of the SCB, linear models of lower 
dimension than the corresponding models of subsection 
4.1. This subsection concludes with the application of 
multilevel control techniques to these hybrid coordinate 
models. The approach presented in this subsection is 
based upon work appearing in Likins (4-7) and Chichester, 
Kaczynski and Nachmias (4-8). 

Subsection 4.4 presents conclusions based upon the work 
presented in Subsections 4.1 through 4.3 and -recommendations 
for future work in the application of multilevel control 
to models of flexible space vehicles. 

4.1 LINEARIZED THREE-DIMENSIONAL DISCRETE MASS ROTATIONAL 

MODEL FOR CONFIGURATION 1. 

Configuration 1 of the SCB consists of two solar wings 
connected to four contiguous modules which are considered 
to be rigid- relative to the solar wings. The simulation 
model for this configuration was derived using the 
discrete coordinate approach 'and approximation by five 
rigid bodies. Each solar wing was converted to two 
rigid bodies with overall designations as shown in 
Figure 4-1 . 
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*>. (The individual bodies are connected to each other 

i 

03 via three degree of freedom spring hinge suspensions. 

Hinge numbers are circled.) 

.TOPOLOGICAL DIAGRAM OF SCB CONFIGURATION 1 
. FIGURE 4-1 





4.1.1 


Definitions 


A 


R. . 
i, J 


Skew-symmetric matrix of the Gibbsian vector A. 
0 -a„ a. 

a. 


A A 


z 


L_ a y 


Z 7 

0 -a 


x 


a 


x 


0 


Vector from C.M. of body i to hinge j in body i 
coordinates . , 


ok Inertial angular acceleration vector of body i 
expressed in body i coordinates. 


A. . 3x3 submatrix of the 15 x 15 coefficient matrix 

) J 

as shown in paragraph 2 . 

-T- 3x3 transformation matrix from body i coordinates 

J 

to body j coordinates. 


m. 


EM 


J i 


. T 


W Ri 5 j 


T . 
ei 


Mass of body i. 

Sum of masses of all five bodies. 
Inertia matrix of body i. 

Matrix transpose of A. 

Aw. x (w. x R* . ) 

1 i j j 

External torque vector on body i. 
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T . 
ai 


Actuator or control torque vector on body i . 


T . Suspension torque vector between bodies i and j . 

S X , J 

4.1.2 Rotational Equations of Motion 

The rotational equations of motion were derived using 
the method of discrete coordinates and take the gen- 
eralized form: 


— 


“ — 


r~ — 

, a l,l 

a l , 15 

lx 


c i 

i 

1 

1 

1 

1 

i 

i 

i 

i 

l 

i 

co 1 
ly- 

i 

i 

i 

. 

c 2 

i 

! 

1 

f 

1 

I 

1 

1 

a l5,l 

i 

i 

l 

i 

a i5;i5 

i 

i 

i 

i 

• 

O - 

5z 


1 

1 

C 15 





_ . 


For convenience, (4-1) can be expressed in terms of 
3x3 submatrices and Gibbsian vectors as : 


"A A A 

1,1 1,2 a 1,3 

i 

a i,4 A l,5 _ 


l 

r—t — 


V 

A A 

a 2 , 1 A 2,2 

i 

1 

1 

1 

1 

t 


i 

! 

1 

= 

! 

1 

1 

1 

i 

L A 5,1 

1 

- - - A 

5,5J 


1 

1 

CD ' 

~ 5J 


l 

in 

o 

i 
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(4-3) 


For the linearized discrete mass model, 


.T. 
i 3 


1 0 
0 1 
0 0 


0 

0 

1 



i and j 


Since the elements of the linearized coefficient matrix 
are constant, the matrix inverse need be computed only 
initially. Appropriate symmetry conditions for the 
linearized model of Configuration 1, may be stated as 
follows : 


m 3 

— 

m 2 

m 5 

= 

m 4 

*3 

= 

J 2, 

X 5 

= ' 

X 4 

R l,3 


R l£ 

R 3,3 

= 

R 2£ 

R 3,5 


R 2,4 

R 5,5 


R 4,4 


Using skew-symmetric matrices as defined and the above 
listed symmetry properties, the 3x3 sub-matrices of 
(4-2) take the following form. 
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= Ij -2(m 2+ m 4 ) [l-2(— )J R 1>2 R 1>2 (4-5) 

m 2 +rn 4 " 

= I 2 + (m 2+ m 4 ) C-^j- -1) R 2>2 R 2j2 (4-6) 

m 2 + m 4 ~ ~ 

+m 4 ^'ZM ){ R 2,2 R 2,.4 + R 2 , 4 R 2 , 2 } 

m 2 ~ ~ ' 

+m 4 C ZM _1) R 2 , 4 R 2,4 = A 3,3 


I 4 + m 4 C ZM 1) R 4,4 R 4 , 4 ~ A 5,5 

m +m 

( m 2 + m ^) ( l - 2 (- ■ )} { R R - 


(4-7) 


1,2 R 2, 2 + R l ,2 R 2,4 } 


= A, 


(4-8) 


m +m 

m 4 2 (-zm-)) %, 2 R 4,4 " '*1.5 

A T = A T 
4,1 A 5 r l 


m 2 +m 4 ~ 

(m 2 +m 4 .) (— gjp-) r 2) 2 R 2,2 


m 2 +m 4 

ZM ^ m 4 ^ R 2 2 R 


2,2 R 2,4 + R 2,4 R 2,2 } 


(4-9) 


'ZM. R 2,4 R 2 ,4 A 3,2 


m 2 +m ~ ~ n, 

m 4 ( ZM“ _1) R 2 , 2 R 4 , 4 + m 4 (1 "ZM ) R 2,4 R 4 , 4 


(4-11) 
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A 2,5 

= m 4 (- 

EM ' 

R 2,2 R 4,4 EM R 2,4 R 4,4 

(4-12) 


= A 3,4 

_ 

A T = A T 

5 , 2 R 4 , 3 



m.2 




A 4,5 

_ 4 
EM 

R 4,4 

R 4,4 = A 5,4 T 

(4-13) 


In the expansions of the 3x3 linearized coefficient 
submetrices of (4-5) through (4-13) triple vector products 
of the following form occur repeatedly. 

R lk xfij x Bjp - - [r 1]£ Rj,, ] i> } (4-14) 


These triple vector products may be transformed to the 
following form in order to derive the expressions for 
each element of the 3x3 coefficient submatrices 





" 


p — i 

e ll e 12 e 13 


ID . 

. JX 

e 21 e 22 e 23 


w . 
jy 

e 31 e 32 e 33 


^jz 



m — 


(4-15) 


The elements of the coefficient matrix in (4-15) can be 

expressed in terms of the scalar components of R^ k and 

R . „ as follows. 

3 & 


'11 


R iky 


R 


j£y 


R ikz 


R 


jHz 


e 12 R iky R j £x 
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e i3 

R ikz 

R • o 
3 Ax 

e 2 l 

_R ikx 

R • n 

j&y 

e 22 

= R.. R . 
ikx j £x 

e 23 

= -R . , R . „ 
xkz 3 Zy 

e 31 

- - R ikx 


e 32 

~ R iky 

R j£z 


(4-16) 


e o0 — R., R. n + R., R., 

33 xkx j£x iky j£y 


The linearized vector components of the right side 
of (4-2) are the following vector sums: 


C 1 T el + T al T sl,2 T s1,3 


C 2 T e 2 + T a 2 + T sl , 2 T s 2,4 


'3 = T 0 + T . 0 + T 1 0 - T n _ . 
e3 a3 si , 3 s3,5 


(4-17) 


C 4 T e4 + T a4 + T s2,4 


C 5 T e5 + T a5' + T s3,5 


4.1.3 Euler Angle Equations 

Relative Euler angle rates between bodies 1 and 2 , 
bodies 1 and 3, bodies 2 and 4 and, bodies 3 and 5 
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are used based on a ip, 6, 4> rotational sequence about 
the Z, Y and X axes respectively. Euler rates are also 
determined relating body 1 to inertial space. Assuming 
small angles, the rates become: 

A( i>l,2 = “2x _w lx 

A ®1,2 = W 2y- W ly (4-18 ) 

^1,2 = W 2z -W lz 


A< ^1 , 3 = W 3x" W lx 
A ®1,3 = u 3y" W ly 
A ^l,3 = w 3z -w lz 


(4-19) 


A< ^2, 4 = W 4x -W 2x 
A ®2,4 = W 4y -W 2y 
^2,4 = W 4z~ W 2z 

A< ^3 , 5 = W 5x -W 3x 
A S,5 = W 5y -W 3y 
A *3,5 = “Sz^Sa 


(4-20) 


(4-21) 
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(4-22) 


+1 = “lx 
®1 “ “ly 

4 

tJj _ w- 
1 lz 


4 -. 1.4 Suspension Equations 

The following torque equations are given in terms of 
rate and displacement constants (program parameters) and 
inertial rate differences and relative displacements. 


T sl , 2x “ " K sl , 2x Ai|) l , 2~ C sl , 2x (w 2x “l x ) 

T sl , 2y = _K sl,2y A0 l,2 _C sl,2y (w 2y' a) ly ) 
T sl , 2z “ -K sl,2z A ^l,2 _C sl,2z (w 2z" a5 lz ) 

T sl , 3x = " K sl,3x A(f) l,3“ C sl,3x (w 3x _ “lx ) 
T sl , 3y = ~ K si,3y A0 l,3~ C sl,3y (w 3y~ u ly) 

T sl , 3z = _ *Ssl , 3z A ^l , 3~ C sl , 3z (a) 3z -w lz ) 

T s2,4x “ _K s2 ) 4x A ' ,, 2,4" C s2,4x (a3 4x'' a) 2x) 

T s2,4y = “ K s2 3 4y Ae 2,4" C s2 J 4y (: “4y _a) 2y ) 
T s2,4z = r ' K s2 3 4z A ^2,4" C s2 J 4z^4z“ W 2z ) 


(4-23) 


(4-24) 


■ (4-25) 
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(4-26) 


T s3,5x K s3,5x A< ^3,5 C s3,5x (w 5x w 3x^ 

T s3, 5y = ~ K s3,5y A0 3,5" C s3 } 5y (w 5y" W 3y ) 
T s3,5z. = “ K s3,5z A ^3,5" C s3 J 5z (a) 5z -a) 3z ) 


4.1.5 • Rotational State Equations 

The rotational equations of motion, (4-1) and the 
Euler angle rate equations (4-18) through (4-22.), may 
he transformed to scalar state variable form as follows 


Definition of State Variables 


Y0 1 

— 

W, 

lx 

Y0 2 

r= 

w, 

ly 

Y0 3 

“ 

1 Z 

Y0 4 

= 

to . 

■ 2x 

yo k 

o 

— 

UU 

2y 

Y0 6 


w 2z 

yo ? 

— 

U 3x 

00 

O 


W 3y 

Y0 9 

= 

to„ 

3z 

c 

i— i 

o 

== 

CO 

4x 

Y0 11 

- 

W 4y 

Y0 12 

= 

W 4z 

Y013 

= 

W 5x 


(4-27) 
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Y °14 

— 

W 5y 

Y0 15 

= 

u '5z 

Y0 16 

= 

A( f ) 2,4 

Y0 17 

= 

A0 2,4 

Y0 18 

— 

A *2,4 

Y0 19 

= 

A(f) 3,5 

Y0 20 

— 

A9 sp 

Y0 21 

= 

A *3,5 

Y0 22 

= 

A< ^1,2 

Y0 23 

_ 

A9 1^ 

Y0 24 

= 

A ^l 7 2 

Y0 25 

= 

A ^l,3 

Y0 26 

= 

A0 1,3 

Y0 27 

= 

A *l,3 

Y0 28 

= 

h 

Y0 29 

= 

9 1 

Y0 30 

= 

h 


m 1 = (^o 1} 

Y0 2> “ 

-- TO 15> 

(4-28) 

yo 2 = (Y0 16 

> yo 17 , 

T 

YO ) 

(4-29) 

FT. = C. 
1 l 

i=l , 2 , - 

— ,15 

(4-30) 

FT = (FT^, 

ft 2 - 

- PT 15> T 

(4-31) 
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Using the state variable definitions in (4-27), the 
right hand sides of (4-18) through (4-22) may be 


written in the following form: 

D i = YO i-6- YO i-12> 1 * 16 ’ 17,-21 (4-32) 

D ± = Y0 i _ is - Y 0 i _2l ; i = 22,23,24 ' (4-33) 

D i = YO i-18 _YO i-24 ; 1 = 25 » 26 » 27 (4-34) 

D i = Y0 i-27 ; 1 = 28 > 29 > 30 (4-35) 


Substitution of (4-27) through (4-31) into (4-1) yields 
the following: 


Y0 1 = A 1 FT (4-36) 

where 

a l,l a l ,15 

i i 

A = | 1 (4-37) 

t i 

! I 

t 

- a 15 , 1 a 15, 15 
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With the substitution, 



b l , 15 


1 



(4^38) 


'15,1' 


15,15 


(4-36) may be written in the following form. 

. 15 

YO = E b ,FT (YO, ,Y0 9 ) ; i = 1,2, 15 (4-39) 

x j=i X > J J x ^ 

Substitution of (4-27, (4-29)- and (4-32,) through (4-35) 
into (4-18) through (4-22) yields the following form 
of the Euler rate equations. 

Y "°i = £>. (X^, Y0 2 ); i = 16, 17, 30 (4-40) 

The set of equations (4-39) and 4-40) constitutes the 
scalar rotational state equations for the discrete mass 
mathematical model of SCB Configuration 1. This set of 
equations is linear and coupled in the state variables, 
YO^. Since A is a constant matrix, A b is also and 
each YCh is a linear function of the scalar components 
of IQ X and YO^ . 
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4.2 APPLICATION OF MULTILEVEL CONTROL TECHNIQUES TO A 

LINEARIZED DISCRETE MASS MODEL 

The overall approach utilized in applying multilevel 
control techniques was outlined in a memorandum by 
Chichester (4-6) as follows: 

1. Decompose mathematical model into a set of decoupled 
equations . 

2. Construct performance index and form Hamiltonian. 

3. Develop costate equations with associated coordination 
equations . 

4. Develop control algorithm. 

5. Construct subproblem hierarchy. 

6. Discretize equations of each subproblem. 

4.2.1 Decomposition 

The first step in the application of multilevel control 
techniques to a system in coupled state variable form 
is the decomposition of the model to temporarily 
suppress the coupling between the equations of which it 
is comprised. With the omission of external and actuator 
torque terms, (4-40) may be written in the ■ following form: 

YO i = DTS i *YO i + TSC i (YO k ),k f i , i = 1, 2, 15 (4-41) 
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where : 


15 


3PT. 

= 2 
1=1 

b. . 
i, J 



3YO . 

l 

15 


27 

= .2 

b . . 

2 - 

J=1 

1,1 

k=l 



k^i 


3Y0 


k 


*Y0, 


27 


2 DTSC. . YO, 
x,k k 


k=T 

k^i 


DTSC. . 
i.k 


3TSC i 

3YO. 


15 

2 


b . . 

j=l 1,J 3YO. 


3FT. 
1 


( 4 - 42 ) 


( 4-43 ) 


(4-44) 


Prom (4-17) and (4-23) through (4-27) it is evident that 
the partial derivatives in (4-42) and (4-43) are constants 
that can be determined. 


All of the coupling terms in (4-39) are summed in (4-43 ). 
The rotational state coordination variables may thus 
be defined as follows: 

27 

= 2 DTSC ijfe . *Y0 k i - 1, 2 , 15 (4-45) 

k — 1 * 

■ k^i 


With the external and actuator torque terms restored., 
the decomposed' form of the first 15 rotational state 
equations may be written in the following form: 


YO. = DTS. 
i x 


*YO . 

l 


+ D. 


15 

2 

j=l 


ij 


(T . +T .) 
ej aj ' 


(4-46) 
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From (4-32) through (4-35) 



0 for -i = 16, 17, — ,30. 


(4-47) 


Hence, (4-40) is In decomposed form with as state 
coordination variables. The rotational scalar state 
equations are given by (4-39) and (4-40) with the 
corresponding state coordination equations (4-32) through 
(4-35) and (4-45). 


4.2.2 Performance Index For The Linearized Model of SCB 

Configuration 1 With Control Torque Applied 'Only To 
Body 1 . 

For :• 



T alx 



ta 3 


T 

alz 


P 


6 

Z 

i=l 



(4-48) 


• (4-49) 


p. = W i *(Y0.) 2 + W i+6 *(TA i ) 2 


i — 1 , 2 , 3 


(4-50) 
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p. = V (TO 1+24 -Y0* +24 ) 2 i - 4, 5, 6 


(4-51) 


where: 


W. and W . , n are scalar weighting constants to he 
1 1+6 


selected 


Y 0 v + 24 are specified optimal values of Y0 ^ + 24 for * = 


4.2.3 Hamiltonian for Linearized Mo.del of SCB Configuration 
1 With Control Applied Only to Body 1 


TE 1 = T elx 


TE 2 = T ely 


(4-52) 


TE 3 = T elz 


H = E H. 


(4-53) 


= W i (Y 0 i ) 2 +Wi + 6 (TA i ) 2 


+A . (DTS . *Y0 . +D . + £ b. .T . + S b..T .) 

l xix x,j ej xj aj 


+ B ± (YO i -YO.) 


i = 1, 2, 3 (4-54) 
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where : 


i = 1, 2, 3 


(4-55) 


D. 

l 


* 


27 

E DTSC 
k=l 
k^i 


i,k 


* YO, 


H i = W ± (TO 1+24 -YO* +24 ) 2 + \ ± D* +24 i - 4, 5, 


* * 

D i = TO i-27 


1 - 28, 29, 30 


4.2.4 Costate Equations For Linearized Model of SCB 


• = _3H 

A i 9YCL 


= -DTS i *A i -2W i *YO i ~3 i 


i = 1, 2, 3 


X i ( tf) " 0 


i = i, 2, 3 


* _3H 

A i 3YO . 


= 2 V YO i + 24- YO i + 24> 


i = 4, 5, 6 


A i (t f ) = 0 


i = 4, 5, 6 


4.2.5 Control Equations For Linearized Model of SCB 
Configuration 1 


3H 


( TA i^r+l (TA i ) r q i ( 3TA. ^ 


6 (4-56) 

(4-57) 

Configuration ] 

(4-58) 

(4-59 


(4-60) 
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(4-61) 


3E 

3TA . 
x 


■2W. . C *TA. 
1+6 1 


3 

+ Z 
3=1 



i 


1, 2., 3 


where : 


— scalar weighting coefficients 
r = subscript identifying iteration number 
X^ = ith costate variable 


4.2.6 Additional Necessary Conditions for Linearized System 


3H 

3B ± 7 

•rl 

O 

6 

* 

= YO. 

X 


i = 1, 2, 

3H 

3 

= z x . 

* 

3Dj 

6 

+ Z 

9D !+24 

3YOt 

3=1- J 
3^1 

3YO x 

3=4 

37^1 

3 9Y0* 


(4-62) 


0 


3 , 

<4 DTSC j,ib + X i+3 i = 1, 2, 3 (4-63) 

3 T 
3f± 


3H' 

3YO* 

l 


2W.*(Y0.-Y0*) 
x x x' 


0 ; 


28, 29, 30 


(4-64 ) 


YO* = YO. 
x x 


i = 28, 29, 30 (4-65) 


4.2.7 


Summary of - Coordination Equations for Linearized Model 
27 

D. = Z DTSC . *Y0. i = 1, 2; , 15 (4-46) 

k=l ljK K 

kfl 
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D i - YO i-6- YO i-12 

i = 16, 17, 

, 21 

(4-32) 

D i - YO i-18- YO i-21 

i = 22, 23, 

24 

(4-33) 

D. = YO. ,„-Y0. 
i i-18 1-24 

i = 25, 26, 

27 

(4-34) 

D i = Y0 i-27 

i = 28, 29, 

30- 

(4-35) 

B. = j WBOj.i^^i+3 

i = 1, 2, 3 


(4-63) 

Vector Representation of 

Variables for 

Linearized 

Model 


of SCB Configuration 1 With Control Applied Only to 
Body 1 


TO = 

= (Y0 1 , 

YO 

IU 2> 

— ■ Y0 30> T 


D = 

(D p D 2 

u 

i 

T 

- D 30> 


TA ; 

= < TA r 

ta 2 , 

TAs) 4 

(4-66) 

31 : 

ii 

/— N 

t-3 

H 

i — 1 

te 2 , 

— , te 15 ) t ■ 


A = 

( ^]_» ^ 2 ’ 

i 



I = 

( ^2’ 

T 




TO* = (YO| g ,YO| g ,YO| 0 ) T 
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The subproblem hierarchy lor the linearized discrete mass 
.state variable rotational model of the Space Construction 
Base is depicted in Figure 4-2. 

4.3- APPLICATION OF MULTILEVEL CONTROL TECHNIQUES TO A 

HYBRID COORDINATE MODEL 

An important objective implicit in the construction of 
a mathematical model for the rotation and flexure of 
■ Configuration 1 of the SCB is the development of technique 
that may be readily extended to successively larger 
configurations involving greater numbers of rigid bodies. 
The complexity and high dimension resulting from a 
discrete mass modeling of Configuration 1 in Section 4.1 
motivated a search for a mathematical modeling approach 
that would yield a model of less complexity and smaller 
■ dimension.- Linearization of the discrete mass model 
attained a reduction of complexity but not a- reduction 
of dimension of the model. Hybrid coordinate modeling, 
an approach developed by Likins (4-7), w.as seen to offer 
promise of both linearization and reduction of dimension 
of the resulting model. To evaluate the efficacy of 
this approach a single axis torsional space vehicle 
was modeled in terms of hybrid coordinates and subjected 
to multilevel control techniques. The results of this 
particular study were published in a paper by Chichester, 
Kaczynski and Nachmias (4-8) . 

4.3.1 Discrete Rigid Body Model of a Spacecraft 

The vehicle to be controlled is shown in Figure 4-1 . It 
consists of a central body with appendages such as solar 
panels along each direction of the Z axis. Because of 
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FIGURE 4-2 

SUBPROBLEM HIERARCHY FOR LINEARIZED 
SPACE CONSTRUCTION BASE 

ROTATIONAL MODEL WITH MULTILEVEL CONTROL 
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symmetry, the five bodies of Figure 1 may be represented 
by three masses, following Porcelli (4-9) as depicted 
in Figure 4-3. The equations of motion for the vehicle 
are the following: 


J o 9 o ■ 

- K i<V e i )+ T o 



K i< e o" e i ) - vw 

,(4-67) 

J 2 0 2 

K 2 .(6i-0 2 ) 



where: 0. = angular rotation of the ith body relative 

to an inertial frame of reference 

Tj, = disturbance or control torque on ith body 

J = moment. of inertia of ith body abopt- the Z-axis 

K. = torsional spring constant for spring 
1 between bodies i-1 and i 


4.3.2 Hybrid Coo'rdinate Model 

The standard form of the hybrid coordinate model of the 
vehicle described in (4-67) may be written as follows: 


J t 's 0 - 4 T i = T 0 


n + [ 2 ?aJ ;n + [a 2 J h.-<50 


0 


= 0 


(4-68) 


where: £ - structural damping (same for all modes) 

[ 2 «°] = 2C w , where ^a^is a diagonal 
matrix of modal frequencies. 


T1 


= q = vector of modal coordinates. 
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FIGURE 4-3 TORSIONAL MODEL OF THE VEHICLE 
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(4-69) 

(4-70) 


j = j +J n +J 0 = total vehicle moment of inertia 
t o I 2 

4.3.3 Decomposed Hybrid Coordinate Model 

Expansion of equation set (4—68) into its scalar com- 
ponents, the substitutions X^- 0 q, ^ 2 = ® 0 5 ^ 3 =ri l ’ ^ 4 -r, i’ 
X 5 =n 2 , X 0 = n 2 and the simultaneous solution of the 

resulting equations for X i , i=l, 2 , , 6 , generate a 

coupled state variable hybrid coordinate model of the 
vehicle to be controlled. State coordination equations 
that temporarily suppress coupling between the scalar 
equations of the state variable model may be written 
in the following form. 
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6 * 


(4-71) 


C i = g i CX c ); 1=1 , 2 . , . , 

where: X c = £ x i > • • > X ii:i » X i x>***> x g]] T (4-72) 

The resulting scalar equations of the decomposed 
(decoupled) state variable hybrid coordinate model are 
the following: 





(4-73) 


■ where: x i ('t Q ) = X iQ 


(4-74) 


4.3.4 Decomposed Performance Index and Hamiltonian 

When the mathematical model of the vehicle to be 
controlled is decomposed into subsystems, the corresponding 
decomposed performance index may be written in the 
form of a summation of terms each of which corresponds 
to one of the subsystems of the decomposed model . The 
decomposed Hamiltonian resulting from such a performance 
index also may be written in summation form, for the 
flexible vehicle treated here the decomposed Hamiltonian 
may be written as follows: 

6 

H = Z H. (4-75) 

i=l 1 

where: H ± = P i +X i (f i +u ]L +C i )+p i (g i -C i ) (4-76) 

♦This range for the subscript, i, will apply in the 
sequel unless a specific exception is stated. 
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(4-77) 


p. , 


W. . = constant weighting coefficients 
i j J 


= costate variable for subsystem i 

= Lagrange multiplier for the coordination 
constraint, equation (4-71). 


4.3.5 Costate and Costate Coordination Equations 

Necessary conditions for effecting optimal control are 
stated in terms of partial derivatives of the Hamiltonian. 
Two such conditions for the flexible vehicle treated 
in this section may be stated as follow's: 


-3H 

3X. 

l 



(4-78) 

(4-79) 


The first of these conditions generates the following 
costate equations . 



!£i x 

9X i 1 


6 

Z 

0=1 

jfi 


!£i r 

3X 1 p j 


where: X^(t^) = 0 


(4-80) 


(4-81) 


The second generates the costate coordination equations 
which follow: 


(4-82) 
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4.3.6 The Control Algorithm 

The iterative relationship for successive values of the 
control torque may be expressed as follows : 

(T 0 ) r+1 = (v r -q (lf 0 ) r <4 - 83) 

where r is the iteration index and q is a constant to be 
chosen on the basis of the rate of approach of the per- 
formance of the controlled vehicle to optimality. From 
(4-73), (4-76) and (4-77) it may be seen that the 
partial derivative in (4-83) is a function of 1^. 

Construction of the Subproblem Hierarchy _ 

The state equations, (4-73) and also (4-74) may be 
regarded as comprising a state subproblem. Similarly, 
the costate equations (4-80) and (4-81) may be assembled 
into a costate subproblem, and the control equations (-4-83) 
may be associated in a control subproblem. The state 
coordination equations (4-71) and the costate coordination 
equations (4-82) may be. combined into -an overall -coordina- 
tion subproblem. Then a hierarchy of these subproblems 
can be assembled as depicted in Figure 4-4. 

CONCLUSIONS AND RECOMMENDATIONS 

From the work reported in Subsections 4.1 through 4,3,. 
the following conclusions may be drawn: 

1) The .three dimensional five 'body discrete mass 
rotational model for Configuration 1 of the SCB 
has been transformed in such a way that it may readily 
be. simulated on a digital computer.. The resulting- 


4.3.7 
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FIGURE 4-4- 

SUBPROBLEM HIERARCHY FOR 


MULTILEVEL CONTROL 
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model is in scalar state variable form. 

2) Multilevel control techniques have been applied to 

the scalar state variable form of the three dimensional 
discrete mass model for Configuration 1. The resulting 
hierarchy of subproblems to be solved may be easily 
simulated on a digital computer. 

3) The multilevel hierarchical rotational form of the 
three dimensional discrete mass model for Configuration 
1 of the SCB is non] inear and of relatively high 
dimension. (30 scalar state variables ) The re- 
maining 11 Configurations would require discrete 

mass rotational models of still higher dimension. 

(Up to 120 state variables ) 

' 4) The nonlinear discrete mass rotational models of 
Configuration 1 of the SCB have been linearized 
using the assumption of small angular displacements. 

5) Hybrid coordinate modeling appears to yield models 
that are lower in dimension than the corresponding 
models generated by the discrete mass approach. 

Hybrid coordinate models also tend to be linear 

in most instances. Both of these factors render 
the model more amenable to the application of 
multilevel control techniques. 

6) In extending the combined application of hybrid 
coordinate modeling and multilevel control techniques 
from models of low dimension to models of higher 
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‘ dimension, some of the modal coordinates of the 
hybrid model should be truncated. 

The- following recommendations are offered for further 
work in the application of multilevel control to 
flexible space’ vehicles and the computer simulation of 
the resulting mathematical models. 

1) Simulation of the linearized scalar state .variable 
form of the three ^dimensional discrete mass rotation- 
al model of Configuration 1 of the SCB should be 
completed. 

2) A three- dimensional hybrid coordinate rotational 
model of Configuration 1 of the SCB with multilevel 
control should be developed for digital computer 
simulation. 

3) Candida-te .control, techniques other than multilevel 
control should be applied to either of the three 
dimensional rotational model simulations listed 
above for purposes of comparison. 

4) Since most of the candidate control techniques 
presently contemplated require access tc all of the 
state variables, techniques for generation of 
observers required to sythesize inaccessible state 
variables should be developed. 

5) Modeling, control and simulation techniques developed 
for- Configuration 1 rotational models of the SCB 
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should be extended to other Configurations involving 
more masses . 
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SECTION 5 


5.0 SIMULATION OP LQR CONTROL 

The linear quadratic regulator (LQR) approach is applied' 

• to single axis control of Configuration 1 of the Space 
Construction Base (SCB) . A single axis vehicle model 
is defined and converted to the state variable form for 
application of the LQR technique. A control performance 
index is then defined and weighting factors are selected. 

A matrix Riccati equation for the system is then solved 
as a function of time. This result is used to specify 
the optimal feedback gain coefficients. These feedback 
. gains were used in a simulation of the vehicle control 

system, and the time -responses are presented. A description 
of an alternate application of LQR is then discussed. 

5.1 VEHICLE EQUATIONS 

Configuration 1 of the SCB (References 5-1. and 5-2) is 
shown in Figure 5-1. The vehicle XYZ axes are- parallel 
to the principal axes of inertia. The vehicle is represented 
here by discrete masses with three-axis hinge springs 
between adjacent masses. The vehicle is divided into five 
masses: one central body between the solar wings, and 

each solar wing is divided into two equal masses. Because 
of the symmetry and the absence of product of inertia 
terms, the SCB can be represented as a three-mass rotational 
model as shown in Pigure 5-2. Body 1 is the central body; 
both inboard solar wing sections are represented as body 
2 and both outboard sections as body 3. Module and solar 
•wing dimensions and the inertia properties were obtained 
from Reference 5-1. Center of mass locations are given 
in Figure 5-2 relative to the reference origin for the 
XYZ coordinates. is the center of mass (CM) location 
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SPACE CONSTRUCTION BASE-CONFIGURATION 1 
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FIGURE 5-2 


INERTIA MODEL FOR CONFIGURATION 1 






of body 1 along the X axis and x^ is the CM of the total 
vehicle . 

With the vehicle modeling limited to the Z axis, the 
inertia characteristics of the central body (No. 1) are 
listed in Table 5-1. Key outputs of the table are the 
CM location and the Z axis moment of inertia about the 
CM for body 1. It is convenient to represent the single 
axis vehicle as the torsional model (as in Reference 5-3) 
of Figure 5-3, where the masses are connected by hinge 
springs and K 2 . For the Z axis , 

' ' 2 

J i = T 1 + 5T < m 2 + V R 12 (5 ' 1) 

where = effective MOI of body 1 

1^ = MOI of body 1 about it's CM 

nr = mass of the i-th body 

m = m-^ + m 2 + m^ 

R ^2 = distance from the CM of body 1 to the contact 

point with body 2 (the X axis component only in 
this case) 

The MOI of bodies 2 and 3 are obtained by assuming 
that each solar wing section is a rectangular plane. 

The MOI about the Z axis of each is then 
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TABLE 5-1 

INERTIA CHARACTERISTICS OE BODY 1 

• x = SW i x i /ZW j . = -88.74" = -7.40’ 


MODULE 

WEIGHT 

(LBS) 

POSITION 

(IN) 

RAD. OF 
GYRATION 
(IN) 

HRSH 

MOMENT OF INERTIA 
WRT CM „ 

(SLUG-FTj 


W 

X 

. 

y 

z 

mm 


*1 


16,440. 

66 

o' 

0 

57.8 

+154.7 

96,757.6 

HM 

16,487. 

-156. 

0 

0 

100. 

“67.3 

51,693.7 

TA 

650. 

• 198. 

0 

0 

43.5- 

+286.7 

1 11,795.2 

DM 

5,625. 

-377. 

0 

0 

55.1 

-288.3 

104,578.9 

TOTAL 

39,202 


264,825. 


SM - Subsystems Module 
HM - Habitability Module 
TA - Turret Assembly 


w 

g 



+ 
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FIGURE 5-3 TORSIONAL MODEL OF THE VEHICLE 
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C -'L- 



I 2 = I 3 = (26) 2 26.445/12 = 1490 slug-ft 2 and 
J 2 = I 2 and J 3 = I 3* 

Note th ( at the mass of each section is doubled since it 
represents both inboard or both outboard portions of the 
solar wings. 


The equations of motion for the torsional model (Reference 
5-4) are then defined as 


J9 — sum of constraint and disturbance torques (5-2) 


where J = 


0 

Jr 


0 

0 


2 

0 J, 


(5-3.) 




(5-4) 


Springs and K 2 and dampers D^.and D 2 cause constraining 
torques, while the disturbance is limited to a control 
torque vector u with a component on each body. The ' 
equation in matrix form is then 


J0 = D0 + KG + 


■V_ -r u 


where K = 


-K 


Kr 


0 


1 2 
K ! -( K ! +K 2 ) K 2 


0 


K, 


-K 


2 J 


(5-5) 


(5-6) 
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0 



D 



-<V d 2 ) °2 

D 2 -D 2 - 


— = [ u i 


u 


u 


J T 


(5-7) 


(5-8) 


and the 0^ and 0^ vector are defined in a manner similar 
to 0. K 0 is selected to obtain- a 1 rad/sec resonance 
with J 2 , 

K 2 = J 2 W 2 (5-9) 


is selected to be equal to Kg. Dg is set to obtain 
a damping factor of 0.005 with K 2 and J 2 and is set 
equal to Dg. All numerical values are summarized in 
Table 5-2. Each matrix formed by these parameters is 
printed out on Figure 5-4. 


To place the model in state variable form, both sides 

_1 

of the equation are premultiplied by J 

0_ = J _1 Di + J^Ki + J -1 u (5-10) 


This can be redefined as follows, while also setting 
u 2 = u 3 “ 0: 

*0* = A_0 + A„0 + Bu (5-11) 

— D— n— — 

Matrices J" 1 , A K A^ and the transpose of B are printed 
out in Figure 5-5. 


5-8 



TABLE 5-2 

SUMMARY OF PHYSICAL CONSTANTS 


9 Z AXIS OF CONFIGURATION 1 

« R = (16.50 + 7:40)FT. = 23.9 FT ALONG X AXIS 



MASS 

M.O.I. ABOUT 
EACH C.M. 

EFFECTIVE 

M.O.I. 

SPRING 

CONSTANT 

DAMPING 

CONSTANT 

i 

M. 

X 

(slugs) 

i . 

x 

(slug-ft) 2 

J. 

l 

(slug-ft) 2 

K’. 

X 

(ft-lb/rad) 


1 

1218.21. 

264 ,825. 

293,767. 

1490. 

14.9 


26.445 

1490. 

1490. 

1490. 

14.9 


26.445 

1490. 

1490. 


— 


m i 2 

J 1 = *1 + 5r <m 2 +m 3 )E X2 

m = + nig + mg 

R = DISTANCE IN BODY 1 FROM IT'S CM TO POINT OF CONTACT WITH 
BODY 2 

DAMPING SELECTED TO GIVE A DAMPING FACTOR OF 0.005. 
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FIGURE 5~4 MATRIX PARAMETERS FOR DISCRETE INERTIA, 
SPRING AND DAMPING CONSTANTS 
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•0,100000*01 


BT* 
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FIGURE 5-5 THE INVERSE OF INERTIA MATRIX, 
a k a d AND b t 


5-10 





The system can be redefined as 


• -- 

x = Ax + Bu 

where x - [© 1 0 X & 2 ®2 ®3 ®s] T 
and u = u-j^ 


(5-12) 

(5-13) 


Then A is a recasting of A^ and A^ as follows: 


A 


0 

1 

0 

0 

0 

0 

5072E-2 

-.5072E-4 

. 5072E-2 

. 5072E-4 

0 

0 

0 

0 

0 

1 

0 

0 

1 

.01 

-2 

. -.02 

1 

.01 

0 

0 

0 

0 

0 

1 

0 

0 

1 

.01 

-1 

-.01 


And -the dimension of B is increased to six: 


B 


T_ 



. 34041E-5 


0 


0 0 


0 ] (5-15) 


5.2 


OPTIMUM FEEDBACK GAIN 


The linear optimum feedback control gain can be determined 
by setting a performance index (References 5-5, 5-6 and 
5-7) as follows:' 


PI 


w 

I 


(x T Qx + 


u^Ru) dt 


(5-16) 


5-11 



where Q, R = weighting factor matrices. 


The diagonal of Q is selected to indicate the relative 
importance of each element of x. In this case, 



4 4-4 

where q is selected as = (0.1) =10 , and Wg^ is 
the control loop bandwidth. R is a 1 by 1, where r = (J^ 


) 


-2 


Once the model is placed in the state variable form 

x = Ax + Bu (5-18) 


the matrix Riccati equation which follows was solved 
until P = 0: 

P = PA + A T P - PBR _1 B T P + Q (5-19) 

The LQR feedback matrix for these conditions is defined as 

F = -R -1 B T P = -[f^ F 2 F 3 F 4 Fg F 0 ] (5-20) 

where P is the steady state solution of (5-19). 

A block diagram of the system is given in Figure 5-6. 

The input constants (Q, A-, B, R) used in the solution 
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of the Riccati equation are given in Figure 5-7, along 

with the initial conditions which were used. Note that 

2 4 —1 

q and r were normalized with q = 1 and r = ) 

The final values- of the solution, when t ^200 seconds, 
are shown in Figure 5-8. F is also printed out and 
the time solution provided the following component 
gains of F: 

2906. 

41320. 

15.0 (5-20) 
209.3 

15.1 

209.1 



It is’ interesting to note that: 


F 1 /'J 1 = Wg W « (0.1 rad/ sec) 2 

? 2 /J 1 & 2 S W bw = 2 (.707) (0.1) 


(5-21) 


Thus if the appendages were ignored, the same attitude 

* 

and attitude rate gains would be obtained for central 
body feedback if a standard non-optimal servo technique 
was applied. 


5.3 SIMULATION 

The system was simulated on the digital computer using 
the full LQR feedback matrix for an initial attitude 
error of 
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SOLUTION or MATRIX NICCATI EQUATION' >CA CONRIQ It Z AXIS 

CONSTANT OATA INPUT 


Q « 

0,190090 01 

0,0 


0.0 

0,0 


0,0 

0.0 


0,0 

0,0 


0.0 

0,0 


0,0 

0,0 

Am 

6,0 

0,100000 01 


■ 0 , 907200 » Q 3 

* 0 , 907200*01 


0,0 

0,0 
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0,0 

0,0 


0,0 

0,0 

11 

0,0 
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nm 
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Ts 

0,0 SEC ?• 

0,0 

o.o 

IT *. 

1 7 * 

0,0 

0,0 
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0,0 



0,0 

0,0 



0,0 

0,0 



0,0 

0,9 



0,0 

0,0 


ROOT * 
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0,0 



0,0 

9.0 



0,0 

0,0 



0,0 

0,0 



0,0 

0,0 



0,0 

0.0 


0,0 

o.o 

0,0 


0.0 

0.0 

0,0 

0.0 


0.0 

0.0 

0,0 

0,0 


0,0 

0,0 

0,0 

0.0 


0,0 

0,0 

0,0 

O.o 


o.o 

0,0 

0,0 

0,0 


0,0 

0,0 

0,0 

0,0 


0,0 

0 , 507200 * 02 . 

0 , 507200 * 0 * 

. 0,0 


0,0 

0,0 
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0,0 


0,0 

• 0,200000 01 
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01 
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0,0 

0,0 
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01 
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0,0 

0,0 

0,0 

” 

“ 0.0 


0,0 

0,0 

0 , 0 . 

0,0 

0.0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

• 0,0 

o.o 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

0.0 

0,0 

0.0 

0.0 

0,0 

0,0 

0.0 

0,0 

o.o 

0,0 

0,0 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

0,0 

0,0 

0,0 


FIGURE 5-7 SYSTEM MATRICES A AND B, WEIGHTING FACTORS Q AND R, 
AND INITIAL CONDITIONS FOR SOLUTION OF THE MATRIX 
RICCATI EQUATION. 
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0,511900 00 
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0,141290*02 
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0,711190 91 
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0,259440*02 
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*9,191270*94 

0,4537*0.0* 
0. *44510.03 
•4,4215*0.0* 
■0,211*20*05 
■0,4*1460*97 
.0,1976*0*01 

0,751090*05 

•0,371040*04 

•0, 231*20.05 
0,15*010*0* 
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0,251*10*04 
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0,1144)0.04 
•0,591870*04 
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0,251610*04 
•0,711*10*03 
0,409**0.04 


FIGURE 5-8 "FINAL" CONDITIONS OF MATRIX RICCATI EQUATION 
SOLUTION AT TIME = 200 SECONDS 
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r 


9(0) 


1 degree 

0 

0 


(5-22) 


The 9 1? u, r\^ and n 3 variables are shown on print-plots 
of Figure 5-9 and 5-10, where: 


= relative motion of inboard appendage = 9 9 - 9 1 

2 4 ± ( 5-23) 

n 3 = relative motion of outboard appendage = ®3 - 

Since the central actuator is located at body 1, the 9^ 
transient is well behaved. n 2 an <* n 3 , however, oscillate 
at about 0.1 Hz and appear to be lightly damped. The 
control u also appears as a smooth transient - even 
though it is also dependent upon oscillatory states from 
bodies 2 and 3. All six states have been applied to 
minimize the performance index based upon the weighted 
integration of 9^ and u. 

To examine this further, the simulation was repeated with 

Fg = F^_ = Fg = Fg = 0, since these gains were' already 

very low in the optimal case. The result of this "partial ■ 

LQR" feedback was basically the same as that obtained 

previously for 9^., and ri 3 (Figures 5-9 and 5-10). But 

the control torque u had a 0.1 Hz oscillatory component 

superimposed on the basic response obtained in Figure 5-9. 

These results are shown in Figures 5-11 and 5-12. 

Obviously, the performance index here is not minimal 

T 

since the integral of the u Ru term would approach 
infinity . 
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FIGURE 5-9 FULL STATE LQR FEEDBACK, CENTRAL BODY ATTITUDE 
AND CONTROL TORQUE INITIAL CONDITION: 

0-^0) = 1 DEGREE 

5-17 
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-0,616610 00 , 
•0,676610 00 ■— — 
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•0, 1 0966D 01 , 
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-0,121660 01 t 
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; RELATIVE MOTION OF pODY -2 


: RELATIVE- MOTION OF BODY 


>,100D 01 


50 

(TIME SECONDS) 


FIGURE 5-10 FULL STATE LQR FEEDBACK, RELATIVE MOTION 
OE INBOARD AND OUTBOARD APPENDAGES. 
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5.4 DISCUSSION 

In the previous sub-sections the LQR control technique 
described is a straight forward approach which i-s relatively 
easy to apply once a system is described in the state 
variable form. Some skill is involved in the selection 
of weighting factors (Q and R) for the performance 
index, otherwise the technique is relatively simple for 
linear, non-stochastic systems. 

For the discrete mass and spring model defined here, only 
the full state variable feedback obtained completely 
satisfactory results. When solar wing appendage motion 
gains were reduced to zero, a lightly damped oscillation 
was obtained for the .control torque u-^ although the 
response of the other variables appeared to be unchanged. 
This raises the question as to' how the system would react 
for a multi— dimensional finite element representation 
for the solar wings. Feeding back hundreds of signals 
to a single actuator set seems to be impractical . Re- 
presenting flexible modes of the solar wings by hybrid 
coordinates (Reference 5-8) would be a more reasonable 
approach, where truncation of high frequency or low 
mass solar wing modes may be accomplished in a selective 
manner . 

The LQR control approach could be expanded in several 
different directions: 

a. More elements or dimensions as discussed above by 
modeling appendages from finite element data, which 
may then be represented in hybrid coordinates (Ref- 
erence 5-9). 
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b. Three axis vehicle models including cross coupling 
terms . 

c. Including feedback sensor (and possibly actuator) 
noise sources and generating a feedback filter to 
generate a better estimate of x, usually referred to 
as the linear quadratic Gaussian (LQR) approach. 

d. Attitude stabilize two or more locations on a large 
flexible spacecraft. 

The application of LQR to d. can be demonstrated by a 
single axis representation of a long flexible vehicle as 
in Figure 5-13. If an intermediate mass and angle are 
represented as and 9^, respectively, locations I, k and n 
on the vehicle can be stabilized by selecting the rel- 
ative magnitudes of weighting factors in Q 1 , Q k and 



5-22 




FIGURE 5-13 SINGLE AXIS N-DI MEN S I ONAL REPRESENTATION OF A 
LARGE FLEXIBLE VEHICLE 
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01 


(5-25) 


where Q. 

1 


~ q 2i-l 

0 


0J 


for i = 1, k, n. 


The relative magnitudes of q 1 ^k-l and q 2n-l may be 
selected according to the pointing performance require- 
ments of 9 9, and 9 . If a pointing stability or 

I K II 

"jitter" requirement is specified, at say location m, 

9 m can be reduced by a high weighting of q 2m relative 

to q n t and others . 

H 2m-1 
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6.0 SIMULATION OP MNA CONTROL 

Some- findings regarding the mathematical analysis, 
control system design and simulation of a multivariable 
Nyquist array (MNA) control method are presented here, 
as they apply to control and stabilization of a large 
flexible vehicle. 

The mathematical model for the space vehicle (SCB 
Configuration 1) was obtained using the discrete 
coordinate modeling method described in Section 2. 

The model was then simulated on a Systems Engineering 
Lab's SEL/32 digital computer for dynamic evaluation 
and control system synthesis. 

The model was reconfigured into state space equation 
form for eigenvalue calculations and control system 
design. The full-order three-axis model was reduced 
to three five-body single-axis models to determine 
the effect of cross-axis coupling on the eigenvalues 
of the three-axis model. It was determined through 
eigenvalue assignment, that for control system design 
purposes? each axis could be independently controlled. 
Further analysis revealed that symmetry in the torsional 
axis could simplify the control design effort. This 
condition does not prevail, however, in either the 
normal or lateral axis. 
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Using the symmetric single axis model for the torsional 
mode, a control system design study was performed using 
the MNA method and the retallack pole placement 
procedure. The MNA design method reflects a need for 
single loop bridged-T compensation when a limited 
number of measurements are available to the control 
design. This is due primarily to the location of the 
lightly damped modes within the control system band- 
width. As the bridged-T compensator is a form of 
closed loop pole assignment, the remainder of the 
design study for the torsional axis was performed 
using the Pole Placement method. 

Simulation results for the single axis reduced model 
and the full order three axis model were identical 
when the pole shifts where made parallel to the real 
axis. This type of pole shift is a form of "active 
damping" and was used successfully to achieve 
acceptable levels of closed loop system performances. 

The full scale test results confirm the negligible 
effects of cross-axis coupling from the torsional axis 
to the normal and lateral axis. 

6.1 CONTROL DESIGN EOR THE Z-AXIS 

It was decided at this point in the analysis to 
concentrate the control design effort on the Z axis of 
Configuration 1 for two reasons. First, the symmetry 
assumption detailed above simplified the mathematical 
model for design and analysis purposes. Since the 
design methods to be described later are new to the 
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flexible spacecraft application, a complete under- 
standing of the advantages, disadvantages and limita- 
tions of the proposed methods is required before more 
involved models are used. Secondly, the control system 
design using the reduced Z axis model can be tested on 
the full five-body three-axis simulation and compared 
with the time response of the single axis reduced 
model simulation. A comparison of these time responses 
will provide the experience and knowledge necessary to 
pursue the more difficult X and Y axis designs for 
this and subsequent SCB configurations. 

Using the five-body three-axis model, the transfer 
matrix G^(S) was computed using Danielevsky ' s method 
(reference 6-1). Several of the ninety transfer 
functions were plotted using Bode methods to reveal 
that all flexible modes are located within the control 
system bandwidth. Design problems related to this 
condition are detailed in References 6-2 and 6-3. This 
analysis was repeated for the symmetric Z-axis model 
with the rigid body dynamics removed. A typical magni- 
tude Bode plot is indicated in Figure 6-1, which shows 
the modes as they appear in the system bandwidth. The 
complete set of transfer functions is indicated in 
Appendix B. 

The eigenvalues for the reduced order model (ROM) are: 

a. 0.0 

b. 0.0 

c, d. -.00500564 + jl. 00112 

e, f . -.00160836 + jO. 32045 
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CXE 1) 



DIRECT NYQUIST PLANE 

FREQUENCY RANGE t 0.000 > OMEGA > 3.000 

FREQUENCY INCREMENT » 0.03297 NPTS * 90 

Figure 6-1. A Typical Bode Diagram for the Z Axis 


This set corresponds to a denominator form of: 

S 2 (S 2 + 26 1 w 1 S+w 1 2 ) (S 2 +26 2 w 2 S+w 2 2 ) 

where the roots of each quadratic term are: 

S 1 ,S 2 '= - Su + 

For 6<<1, S^,S 2 = <$w+ jw and 

co = Im( X) 

6 = Re( X)/Im( X) 

where X is a complex eigenvalue. With the eigenvalue 
set above 

(0 = 1.00112 6 = .00500004 

co 2 = .32045 S 2 = .00501907 
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The control system design for the ROM Z-axis model 
will proceed using the MNA method and the Retailack 
pole placement method. The MNA method will be con- 
sidered using only and as the controlled outputs 
with and X 2 g as the system inputs. The pole place- 

ment method will shift specific eigenvalues to 
prescribed locations using full state feedback. 

6.2 MNA DESIGN FOR ROM Z-AXIS 

Using the appropriate transfer functions from Appendix B, 
a direct Nyquist array (DNA) design was initiated over 
the frequency range 0<w<3.0 using the technique 
described In Subsection 3.3. Recent results described 
in Reference 6-4 for finite frequency .Nyquist array 
designs establish this range as acceptable for design 
and stability purposes. 

For this DNA design the post compensator matrix was 
prespecified as the identity matrix with the precom- 
pensator matrix selected in accordance with the 
algorithim in Reference 6-5. The closed loop diagram 
for this axis with 

Q(S) - L G(S) K (6-1) 

representing the forward loop transfer matrix is 
indicated in Figure 6-2, 

where: G(S) = spacecraft vehicle transfer matrix 

L, K = Post- and pre-compensator matrices 
F(S) = Feedback transfer matrix. 
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Figure 6-2. Closed Loop Spacecraft 
Control Diagram 

In Reference 6-6, a procedure is described whereby 
Multivariable Bode diagrams can be obtained from MNA 
results using Reference 6-7. This procedure has been 
incorporated into an interactive design suite in 
Reference 6-8. Results for the ROM Z-axis are presented 
in Figures 6-3 and 6-4. The dominance levels obtained 
using equation (3-32) are very small, thus reflecting a 
decoupled condition as indicated by the virtual 
coincidence of the envelope curves in the Bode diagrams. 

From Figure 6-3, the effect of solar wing interaction 
on the rigid body dynamics is seen to be negligible. 
Hence the control in loop 1 may proceed using a 
rigid body assumption. The dynamics- indicated in 
Figure 6-4 reflect the lightly damped poles at to = .32 
and to = 1.0 and also indicate the presence of a 
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multivariable .zero near w = .32. Since this anti- 
resonance is sharper than that of the pole at w = .32, 
the damping ratio for the zero is less than 0.005. 

One possibility for the control of the mode at w = .32 
is to decrease the damping ratio of the pole via a 
bridged-T compensator and shift the pole location to that 
of the multivariable zero. This approach would 
correspond to a pole-zero cancellation in the traditional 
single-input single-output case. Uncertainties in the 
dynamic model could, however, have a deleterious effect 
on the overall system behavior. 


Alternatively the poles could be shifted outside the 
control bandwidth (active stiffening) and/or reduce 
the pole resonance peaks (active damping). Figure 6-5 
indicates the time response of in the closed loop 

system without dynamic compensation. In Figure 6-6, 
the Bode diagram for loop 2 is presented following 
a bridged-T compensator insertion to increase the 
damping ratio at w = 1.0 to 5 = 0.1. 


P 22< S ) " 


s 2 + .01 S + 1.01 

S a + .2 s + 1.01 


( 6 - 2 ) 


Figure 6-7 indicates the closed loop time response for 

i^ 2 i(t) using (6-2). Retaining the pole at w = 1.0 has 

the effect of reducing the impact of the resonance at 

* 

co = 0.32 as it will appear principally as an envelope 
modulation in the time response. 
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FIGURE 6-7 CLOSED LOOP ij> 2 TIME RESPONSE FOR 
COMPENSATED SYSTEM 
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To examine the shift in frequency of the pole at 
to = 1.0, a bridged-T compensator was placed in loop 
2 to increase the damping ratio to 5 = 0.1 and simultan- 
eously shift the resonance to w = 2.0. This is equival- 
ent to placing a complex zero pair at w = 1.0 and a 
complex pole pair at to = 2.0. The Bode diagram for 
this case is indicated in Figure 6-8 with the correspond 
ing closed loop time response presented in Figure 6-9. 
Note that the zero placement removes the impact of the 
pole at to = 1.0 but increases the effect of the 
resonance at oi = .32. This latter effect is due to the 
pole relocation at a higher frequency; i.e., system gain 
is .effectively reduced by 50%. To recover the effect 
of Figure 6-7 the loop gain should be doubled. 
Alternatively, pole shifts in the resonance at si = .32 
could be used. However, this would ultimately require 
a fourth order compensator for fgg(S) in the closed 
loop control. 

The above analysis and design synthesis for the ROM 
Z-axis model demonstrates the utility of the MNA 
method in the design of controls for flexible spacecraft 
At this point, it was decided to postpone further 
design efforts using 'MNA procedures to examine the 
potential of the Retallack pole placement method. 

6.3 RETALLACK POLE PLACEMENT FOR ROM Z-AXIS 

The dyadic pole placement technique developed by Simon 
and Mitter (.Reference 6-9) and further refined by 
Retallack and MacFarlane (Reference 6-10) utilizes the 
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model structure of a linear system to give the designer 
direct control over the movement of given system modes 
to desired closed loop locations. Unlike other 
methods, the Retallack algorithm is a state feedback 
procedure which allows for the direct assignment of 
specific open loop poles to specific closed loop 
locations. This method is particularly attractive 
from a design viewpoint whenever the open loop poles 
can be directly associated with specific parts of the 
physical system. In this way only those poles directly 
associated with the subsystem will be shifted. 

Let (.j = , n ) represent the open loop poles 

of the system and let p ± (i = 1/ • • • , n) be the desired 
closed loop set. Several algorithms are presently 
available to shift the eigenvalue set { X . } to the 

«J 

desired set {p^}. However, there is no guarantee 
( that any particular eigenvalue Xj will be moved to a 
specific location Without this eigenvalue location 

assignment flexibility, the resulting control in light 
of sensor and/or actuator failures could result in very 
unsatisfactory performance. 

Given the dynamic system 

x = Ax + Bu (6-3) 

and the control law of the form 

u(t) = r - Kx (6-4) 
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where r is the set point control, the eigenvalues 
{p i > of the resulting closed loop system 

x = (A-BK) x + Br (6-5) 

are related to the open loop eigenvalue ( A ^ ) by the 
Hsu-Chen equation (Reference 6-11) 

w (S-p ± ) 

I I - (SI - A)” 1 BK I = (6-6) 

Inn l n (S-A.) 

ir — 3 

j=l 

The Jordan Cononical form for A is given by: 

A = IT 1 AU (6-7) 

where U is the modal matrix for A. Let ( A^, , , ^p) 

p<n be the open loop poles to be shifted and define 
the modal controller K using the unity rank dyadic 
product as 

K = f d T V (6-8) 

where f is an m vector, d is a p vector and V contains 

the p rows of IT 1 corresponding to (A^, , A p ). 

Using (6-8), Retallack reduces (6-6) to 

m (S-P X ) 

l+d T (SI - A) 1 8 = Azl (6_9) 

P 

j=l 

where 

8 = VBF (6-10) 


6-16 



The remaining open loop poles ( \p + -j_>* ’ * •» remain 

unchanged under dyadic feedback. 

To overcome some of the numerical- difficulties associated 
with the computation of the modal controller when 
clustered poles appear near the imaginary axis (Reference 6-12 ) y 
the system in (6-3) is transformed to the Luenberger 
canonical Form (Reference 6-13), 

w = Aw + Bu ; (6-11) 

where w = U - ' 1 ' x and 

A = u -1 AU • (6-12) 

B - U -1 B (6-13) 

An efficient algorithm to realize this transformation 
is presented in Reference 6-l4 . 

Since the eigenvalues of (A-BK) are identical to those 
of A-BK, the pole placement process can be accomplished 
in the closed-loop system of (6-11) rather than in the 
closed-loop system in (6-3). Once the gain matrix 
K has been determined for the canonical system, the 
desired modal controller K for the original system 
(6-3) is given by 

K = KIT 1 (6-14) 

For the ROM Z-axis model provided in Appendix B 

the system equations can be partitioned into two 
distinct sets: one for the rigid body dynamics and one 
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for the flexible appendage. The open loop eigenvalue 
sets are, respectively, 

Set 1: 0.0, 0.0 

Set 2: -.00500564 + j 1.. 00112 

-.00160836 + j 0.320495 

Here was, added to the state vector to produce the 
second pole at the origin. The state variables are 

then 1*21’ w 42’ ^21’ ^42 ’ anc * ^1 wa th controls 
and For the five-body single-axis control 

implementation is set equal to since symmetry 

was used to generate the Z-axis ROM. ■ 

The rigid body eigenvalues were shifted to the following 
positions 

al . r~ . 1 , — . 1 
bl. -.2, -.2 
cl . - . 3 , - . 3 
dl. -.7, -.7 

while the modes corresponding to the solar wing were 
shifted to 

a2 . -.1+j 1.00112 
-.1+j 0.320495 

b2 . -.2+j 1.00112 
-.2+j 0.320495 

c2 . -.7+j 1.00112 
-. 7+3 0.320495 
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The resulting control laws were implemented in a time 
domain simulation of the Z-axis reduced order model with 
^4 (t Q ) = .02 radians. Since all velocities, (t ) 

and i {)^2 (t Q ) are zero, this condition corresponds to an 
equilibrium point with all masses displayed .02 
radians from the desired equilibrium point. In each 
case, the control restored the vehicle to the desired 
location in the manner indicated by the closed loop 
eigenvalue characteristics. .Table 6-1 compares' 
several of the important characteristics of the sim- 
ulator responses. 

The main purpose of this study was to demonstrate that 
the Retallack Pole Placement method can be effectively 
used in the control of flexible space vehicles of the SCB 
configuration. Several observations resulting from 
this effort are indicated below: 

1. As the eigenvalues of the rigid body are moved 
further into the left ha'lf plane with fixed solar 

wing eigenvalues, the peak to peak excursions of 

and increase appreciably with little variation in the 

time to steady state. 

2. Tor fixed rigid body eigenvalues and variable 
solar wing eigenvalues, the time to steady state 
decreases in proportion to the rate of increase in the 
damping ratio with small variations in the peak to 
peak excursions of and 
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TABLE 6-1 


CLOSED LOOP RESPONSE 


if^CO) = 0.02 


EIGENVALUE 

FIGURE 

MAX , 

,MAX , 

TIME TO 

STEADY STATE AND PEAK VALUE 



SET 

NUMBERS 

l T izl 

| T 2zl 


1 

J 21 : 


12 






SECONDS 

RADIANS 

SECONDS 

RADIANS 

SECONDS 

RADIANS 

a l» 

a 2 

— 

53.7 

0.320 

50.0 

0.02 

25.0 

0.00050 

25.0 

0.0016 

b l> 

a 2 

— 

214.82 

1.112 

29.0 

0.02 

30.0 

0.0021 

24.0 

0.0047 

c l 5 

a 2 

— 

483.35 

2.025 

20.0 

0.02 

36.0 

0.0041 

30.0 

0.0079 

a l> 

b 2 

— 

53.7 

0.461 

50.0 

0.02 

21.0 

0.00085 

25.0 

0.0017 

b l’ 

b 2 

— 

214.82 

1.4518 

29.0 

0.02 

25.0 

0.0026 

22.0 

0.0051 

c r 

b 2 

6-10 to 
6-14 

483.4 

2.717 

17.0 

0.02 

25.0 

0.0051 

17.0 

0.0088 . 

d r 

b 2 

6-15 to 
6-19 

8631.6 

9.071 

8.0 

0.02 

18.5 

0.0171 

18.5 

0.0218 

d i- 

°2 

6-20 to 
6-24 

2631.6 

30.637 

8.0 

0.02 

13.0 

0.0207 

12.0 

0.0303 
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a-oooot-oo... 
9.5000E 00 
l.OOOOE 01 
1.0300E 01 
1.100QE 01 
1.1500E 01 
1.2000E 01 
1.2500E 01 
l_3000E_0l — 
1.3500E 01 
1.4Q00E 01 
1.4500E 01 
1.5000E 01 
1.5500E 01 
1.6000E 01 
1.6500E 01 
1.700Q£_0l_ 
1.7500E 01 
1.8000E Ol 
1.05OOE, Ol 
1.9000E Ol 
1.9500E Ol 
2.0D00E 01 
2.0500E 01 
2.1000E 01.. 
2.1500E Ol 
2.2000E 01 
2.2500E 01 
2.3000E 01 
2.3500E Ol 
2.4000E 01 
2.4500E Ol 
. 2.5000E. Ql 


HINIMUH X5 VERSUS TIRE „ 

•<*,7?57E-03 S.8173E-03 

X5 I 1 

0.0 

1.2993E-04 ........... + 

. 8.6814C-04 • + 

2.3718C-03 * — 

4.3998C-03 — ------ — ♦ 

6.4783E-03 — ♦ 

8.Q868E-03 — * •_ — 

8.6173E-03 — — ZZZZZZ.Z. 

6.4753E-03 — — --«— — 

7.J.092E-03 * * 

. -4.974SE-03-. ■ -- - 

2.4534E-03 — — ........ * 

-4.3427E-05 

-2.1676E-03 
-3.6949E-03 •••«♦ 

-4.5503E-03 ♦ 

-4.7957E-03 ♦ 

-4.5904E-03 * 

...4.13751*03 ... — *■ • - 
-3.6310E-03 — — ♦ 

-3.2154E-03 '+ 

-2.9644E-03 

-2.aa08E.03 

-2.9117E-03 * 

-2.9749E-03 
-2.9864E-03 — — ♦ 

_ «2.8834E«03 ■*•»»**♦ — - 

-2.6378E-03 — — + 

-2.2588E.03 

-1..7862E-03 

-1 ,2777E*03 

-7.9307E-04 

-3.8123E-04 — + 

-7.3370E-05 ♦ 

1.2475E-04; + 

2.2530E-04 ——————— 

2.5451E-04 ' 

2.4390E-04 

2.2253C-04 + 

2.1126E-04 — — 

2 . 20 02E-04 

2.4814E-04 --------- — - — ♦ 

2 # 87 09E-04 — ------ — 

3.2449E-04 — 

3.4828E«04 — — ♦ 

3.S002E-04 ----------- — ----+ 

> -3.2665E-04 

2.8073E-04 — + 

2.1919E-04 
1.5129E-04 
— 8-6321E.05 


TIME 

0.0 

5'. DOOOE-01 

l.OOOOE 00 — 
1.5000E 00 
2.0000E 00 
2.5000E 00 
3.0000E 00 
3.5000E 00 
4.0000E 00 
4.5000E 00 
3.0000E-00 — 
5.5000E 00 
6.0000E 00 
6.5000E OO 
7.0000E 00 
7.5000E 00 
8.0000E 00 
B.5000E 00 
.9.0000E -00 
9.S000E 00 
l.OOOOE 01 
1.0500E Ol 
1.1000E Ol 
1.1500E Ol 
1.2000E Ol 
1.2500E Ol 
-1.3Q00E .01 — 
1.3300E 01 
1.4000E 01 
1.4500E 01 
l.SOOOE Ol 
1.5500E Ol 
1.6000E Ol 
1.6500E 01 
-1»7000E-U1- 
1.7500E 01 
1.8000E Ol 
1.8500E OJ 
1.9000E Ol 
1.9S00E Ol 
2.0000E 01 
2.0500E Ol 
2.10QOE. 01- 
2.1500E 01 
2.2000E Ol 
2.2S00E 01 
2.3000E 01 
2.3500E, Ol 
2.4000E 01 
2.4500E Ol 
jt.soooe-jai.. 


RINIMLK U1 VERSUS time 

-4.0335E 02 6.S334E-01 

U1 1 * 

-4.8335E 02 + 

-3.5362E 02 * 

..-2 .5065E 02 ... ---------- — . 

-1.69511 02 

-l.o&iiE 02 — — — i— — I— * 

-5.7080E 01 — — — — — « — ” — t— + 

-1.9652E oi — — — — — — 

8.4568E 00 — — ----- 

2.9116E 01 IIIIII.- - -------- — + 

4.388&E oi — : i — — — .1 

6.0338E oi — ..... .... 

6.3918E 01 ZZZZZZZ Z-Z 

6.5330E 01 — — - .* 

6.5108E Ol '‘"I" .. -------- 

6.3681E 01 

6.1388E 01 Z .... — .+ 

. ..5.5223E- oi- - ~zzzzzz.z~z.z~z- 

5.1724E oi — — — ..... 

4.4532E 01 I - ....... / 

4.1003E 01 HU. 1. 1-.- ....... I-.....-. I 

3.4338E 01 — TTTTI ~~~ZZZZZZZZZ.Z — — ♦ 

— 2.8373E- oi — ««»»»*>• II-.— -—-H— I-———; 
2 . 3194 E oi — — : 

1.6869E 01 .......... . . . 

i,5ii6E oi ZZZZZZZZZZZZZZZZZZZZZZZZZZZ-Z— - 

iTo78ol oi 

9.6055F 00 - — ' T 

8.5494E 00 — -• — * ’ 

6.7512E 00 — — — — — ’ . Z--ZZZZ-ZZ I-.I-.t 

5.3108E oo — — — ZZZ Z-ZZ.Z.-ZZ 

4.1635E 00 —————— — ' ’ 

3. f 8221 00 - T 

3 . 2542 E oo ;-r 

2.8740E 00 

2.53651 00 * I 

2.2373E 00 — 

. -2.237.7e 00 • ” * 



FIGURE 6-14 CONTROL TORQUE BODY 2 


SS-9 


sit. -i 


iC" 


e? 

c£ 


' t S o 

i&9 




TIME 

0>0 

5.0000E-01 
_l„OOOOE-gO~ 
1.5000E OO 
2.0000E DO 

itsoooe oo 

3.00D0E 00 
3.5000E 00 
4.0000E 00 
4.5000E 00 
5 • 0 0 0 0 E_0Q — . 
5.5000E 00 
6.0000E 00 
6.5000E 00 
7.0000E 00 
7.S000E 00 
S.QOOOE 00 
O.5000E 00 
5.0000E-00_ 
9.5000E 00 
1.0000E 01 
i.osooe 01 
1.1000E 01 
1.1500E 01 
1.2000E 01 
1.2500E 01 
_L, 3000 E_ 01 ^ 
1.3500E 01 
1.H000E OJ 
1.450GE 01 
1.5000E 01 
1.5500E 01 
1.6000E 01 
1.6500E 01 
_1*7000E-01- 
1.7500E 01 
1.8000E 01 
l.BSOOE 01 
1.9000E 01 
1.9500E 01 
2.0000E 01 
2.0500E 01 
S.1000E Cl 
1.1500E 01 
2.2000E 01 
2.2500E 01 
2.300QE 01 
2.3500E 01 
2.4000E "I 
2.4500E 01 
2.5000E 


MINIMUM 02 VERSUS TI " E 

-2*6960E 00 

U2 l . 

0. 0 T 

-1.9498E 00 * 

-2-.663SE 00 — ‘ 

-2.6969E 00 ♦ 

-2.3757E 00 

-1.9744E 00 *!!!!*.. 

-1.6190E 00 

-1.0743E 00 ~ ZZZZZ.—---Z~-+ 

- 7 .a 066 t -01 1' I , ♦ 

*3.995SE.01 Z...Z-Z---...Z 

6.6272E-01 ......... — — ' 

1.2704C 00 , 

1 . e 4 o 5 r oo — — — — 

hinii 

^^8479e -oo 

2.4268E 00 """" ............. 

2.1019E 00 """".I — 

i,7270e oo 

i.oiooe oo — — — + 

7.2370E-01 I" + 

__3.1734E*01 

1I7339E-01 

4.7594E-02 

-7.2564E-02 

-1.9324E-01 — ■*■**. 

_j|:!23lf:oi— - 

-6.4724E-01 — 7ZZZZ.—.-.. -♦ 

-8.5959E-01 - '"I . 

-tl^isbE-oi ZZZlllllll 

-5.3979E-01 ~ ZZ.ZZZ..... 

•4*o528E-gi 

.--•3.4291E-01 — ...... 

-2.6850E-01 ZZ.Z.ZZZ 

-2,425iE-ni 

-2.0364E-01 — — — + 

-1.6967E-01 — .... 

-1.383!iE-01 ~ ZZZ-. 

-1.0804E-01 ZZZZZZZ 

-7 ■ 8O58E-02 ZZZ. 

«Jt-a7tJ£E-02-.. ^ ' 


MAXIMUM 
^. 717 lE- 00 




FIGURE 6-15 ANGULAR POSITION BODY 1 FIGURE 6-16 RELATIVE ANGULAR POSITION BODY 




HiNimm 
4.5403E-16 ■ 


VERSUS tihe 


X6 I 

2.0000E-02 

x*9027e”02 — 

-4.6884E.02 — - — ..... 

1.4347E-02 — -- — 

X.1837E-02 

9.5575E-03 

- 7.5922E-03 

5,?54gE»93 

•»IS2l5r»03 

3J5567E-03 —•••••♦ 

..2.7178E»03 — »-«•».+- - - 
2;0641E-03 

- g:?iiS g :88 

4,0812E»O4 «♦ 

3.6222E-04 ♦ 

2.6810E»04_. 

1.9799E-04 ♦ 

1.4590E-04 ♦ 

1.0732E-04 ♦ 

7.8794E*05_ *- 

5I7759E-05 + 

4.2276E»05 ♦ 

Itlls^l.ot *- 

1*6447E»05 ♦ 

1 J J.978E*05 ♦ 

8,7143E«06 ♦ 

-6.3337E-06 
4,5991E“06 ♦ 

3.3367E-06 ♦ 

2.4188E«06 ♦ 

--l:22i2g:8i-r 

9.172SE-07 ♦ 

6,6301E»07 + 

4.7094E-07 - * 

3,4578E«07 + 

2.4948E-07 ♦ 

i;7991E«07 ♦ 

. . -1.2967E.07 *■ 

9.3413E-08 + 

6;7263E.CP ♦ 

_j:S 8 yesi. : 

2.5045E-08 ♦ 

1.8B03E.08 ♦ 

1.2936E-08 ♦ 

9.2916E.99 — * 


MAX I HUM 
2.0000E-02 


IRE 

o 

OOOOE-Ol 
OOOOE-OO 
5000E 00 
OOOOE 00 
5000E 00 
0000C 00 
5000E 00 
OOOOE 00 
5000E 00 
OOOOE-OO 
5000E 00 
OOOOE 00 
5000E 00 
OOOOE- 00 
5000E 00 
OOOOE 00 
5000E 00 
OOOOE 00. 
SOOOE 00 
OOOOE 01 
0500E «1 
1000F 01 . 
1500E Oi 
2000E Ol 
2500E Ol 
3000E-01- 

33881 » 

4500E Ol 
SOOOE Ol 
S500E 01 
6000E Ol 

6500E Ol 
7000E- Ol.. 

28881 8* 
8500E «1 
9000E- Ol 
9500E Ol 
OOOOE Ol 
0500E Ol 
1000E Ol 

i§88 8} 

2S00E 01 
,3000E.0l 
3500E 01 
4OO0E Ol 
4500E Ol 
SOOOE -Ol 


4 
0 

314UE-04 

8993E-04 - 

0811E.03 

2966E-03 

1590E-02 

2467E.02 • 

0293E-02 

3975E-03 

1339E-03 


HINIHLK 

1.2467E-02 


VERSUS TIHE 


■•»+ 


8445E-03 — * 

3327E-02 — 


6b6SE-02 -• 

7140E-02 
5228E.02.- r? 
1499E-02 

8072E-O3 
3504E«03 — 

3416E»03. — 
7372E-03 
7419E-03 
5697E-03 
6348E*03 4* 

4194E-03 
3506F-03 
1377E-04 — 

1535E-04 

9402E»04 *• 

6689E-03 

4065E"03 

9879E-03 

2b0lE“03 

m 3E-03 

2E-03 

0546E.03 4* 

2865E.03 
6049E-04 
563SE-0S 
8417E.04 
3930E-04 -• 

2033E»04 *■ 

9352E.04 *■ 

3217E.04 “.<• 

8740E-OS 
2779E.05 -» 

6472E.0S' -• 

2733E-04-- — 
4613E-04 
6728E-04 
3873E-04 




****** — 






FIGURE 6-17 RELATIVE ANGULAR POSITION FIGURE 6-18 CONTROL TORQUE BODY 

BODY 4 


SS-9 


TIME 

0.0 


OOOOE-Ol 

ooooe ue 

SOOOE OD 
OOOOE 00 
SOOOE 00 
OOOOE 00 
SOOOE 00 
OOOOE 00 
5000E 00 
, OOOOE 00... 
,3000c 00 

S88BI 88 

OOOOE.OO. 
SOOOE 00 
.OOOOE 00 
,3000e 00 

OOOOE-OO . 

5000E oo 
OOOOE 01 
.OSOOE Ol 
.1000E Ol 

JS88I 81 

2500E Ol 
SOOOE Ol.. 
,3500E Ol 
,4000£ Ol 
45Q0E Ol 
SOOOE Ol 
SSOOE Ol 
,6000E Ol 
, SSOOE Ol 
.7000E Ol 

n 

S500E Ol 
9QGQE Ol 
9S00E Ol 
OOOOE Ol 
0500c Ol 
100GE Ol 
1500£ Ol 
2000E Ol 
2500E Ol 
SOOOE Ol 
SSOOE 01 
OOOOE Ol 
,1!500 e Ol 
,5000E-0l 


MINIMUM X5 VERSUS TIME MAXlHim 

-2.0169E-02 S*1702e*O2 

XS I 1 

0,0 

6.4039E-04 . 

3.fl735E-03 — 

1.S(j12C-02 — — — 

2.0460E-02 * : — * . 

2.1702E-O2 

1 .5i01E»02 - .... 

l!2S54E*02 

5.3474E-03 

_-4.672‘(e-03 . t-4 

-1.2502E-02 
• -1,7670E»02 ••+ " 

-2.0169E-02 ♦ 

«.1,9495E«02 - ♦ 

-1.6520E'02 «■— *♦ 

-ll2247E-02 - 

-7.7429E-03 * 

- -S.895SE-03 -- 

-1.2564E-01 : +_ 

1.246SE-05 — •■*+ 

9 .4952E-0S 

. *5.9257e-04 T--T---* 

.1.54S9E-03 — 

-2.309SE-03 

,2;5765E-03 

_-.2,2330E»03 T-T-4** 1 *— 

.1 .3509E-03 •* 

-1.3676E-04 + 

1.1411E-03 + 

2.2277E-03 — — — + 

2.9394E-03 

3.19S0E-03 T+ 

3 .O224E-03 

2,5269E*0S- - 

6.2954E-04 + 

2.6052E-04 
9.6769E-05 
1 . 0363E-O4 

2.1453E-04 + 

3,S302E*04 --------- — 

4.5284E-04 

4.7180E-04 

3.9744E-04 

2,'44S6E*04 

4.6920E-05 — ... — + 

-1.5417£*04 

T3.2094E-04 — ♦ 

-4.2748E-Q4 — 


TIME 

0.0 

5. OOOOE-Ol 

i , OOOOE 00 

_:8oo°o°i 88 


3, OOOOE 00 
3, S000£ 00 
4,5000? 00 
4, SOOOE 00 
5, OOOOE 00 
5. SOOOf 00 
s,5oooi oo 
6, SOOOe UO 
7, OOOOE oo 
7, SOOOe 00 
0.0000? 00 
fl.SOOOg 00 

sloooog 00 
9, SOOOE 00 
l.OOOOr Ol 
i;o50oe oi 
i.iooot: gi .. 
i.isooe oi 
i: 2 ooo£ gi 

1.2S00E Ol 
.1.3000E Ol— 

l.ssggE gi 

1 ,4000e Ol 

I»4500E Ol 

1.5000E Ol 

llssggc 2 1 

l.figgOE 5l 

1.6500E 
1.7000E 
i;7500|: 
1.8000E Ol 
1, SSOOE Ol 
1.9000c Ol 
1.9S00E Ol 

!:8s88E 8i 
!:!jP « 
i:iS88i 

1:18881 8} 

2 »4oggE gi 

2.4S00E gi 

2.5Q00E .01 


Ol 

01 . 

01 


Si 

i 


MINIMUM U1 VERSUS TIME MAXIMUM 

-2.6316E 03 3,544|e 02 

-2.6316E 03 ♦ 

_ _ ........ — ......... — ......................... + 

,4295E 02 — — 

3.5446E 02 — — - — .... — . — — — 

2 ,Sao 4 E 02 »-f 

2.4244E 02 ....... .............. . + 

- 1.9866E 02 — 

1.5959E 02 ...... — ....... — ........ 4 . 

9.8718E Ol — 

7.6424E 01 - .... — ... ..... . — ..... 

5, 6089c Ol — • — — * 

4.4763E Ol - 

3.3944E Ql — . ^ 

-2.5612E_0l- — 

1.924DE Ol 

1.4398E Ol 

i;o738E oi 
7,9842t 00 - 

S ,9203c 00 — — — * — ....... 

4.3791E 00 ............................... 

3.2319E 00 .— 4 . 

- -2.3803E-00 — s-sf*— 

1.7499E 00 ...... .... — — ♦ 

1.2842E 00 

9;4095e-01 

. 6,8844E»0l 

5, 0301E-01 .......... ...... ............ ..o.. + 

3.07O6E-OI 4 

2.6754E-01 t— — 

i;4167E-Sl ... ....... 

ll0294C-0l ...... .... -♦ 

7.4731EO02 — ......... — ....... 4 

5.4204E-02 . 

3,9284E»02 

f ,6449E-02 

.S5B8E-02 .... .... 

1,4886E*02 — ........ — ,.... — ............ 

1.0760E-02 ..... . 

sisimn 

4, 0468E-Q3 ...... ............... + 

2.9177E-03 

2.1026E-03 ......... .... ... 

1.5145E-03 ...................... 

- 1.0904E-03 - 


FIGURE 6-19 CONTROL TORQUE BODY 


9S-9 


to 


00 

00 

00 

00 . 

00 

00 


time 

o,o- 

5.0000E-O1 
■*1«0000C DO - 
1.S000E oo 
a.ooooE oo 

а. sooaE oo 

3.0000E 00 

5:58881 “ 

4,sooor 

5.0000L 

1:58881 „ 

б. 5000c 00 

,7«ooooe oo * 
7,5000c 00 
e.ooooc oo 
gisoooE oo 

9,OOOOE -00- 

?;5888i |S 

lI0500C «1 

i:i888i oi • 

i:«OOE Ol 
l,3000r-0l- 
1.3500E 01 

ijnoooc 01 

lI<*500E 01 

l-.5G00C.0l 
1.5500E 01 
1.6O00E Ol 
iU^OOF 01 
l,7000C-Ol 

1:18881 1 

..1:18881, !i 

1;!888E 81 

2.0500E Ol 
2.1000E °1 
2* 1500C Ol 
2,2000E Ol 
II2500E 01 
2,3000C Ol 
2*3500E 01 
2.4000C Ol 
2.450QC Ol 
2.3000C Ol 


minimum ' 

■9,071SE 00 


U2 


VErSUS time 


U2 

0,0 

-0 ,5630C 00 «♦ 

— 9.0715E 00- — * 
-6.4126C 00 


maximum 
6,4731c 00 


---- 


.. 

- 3 , 3066 c 00 

-9.7164E-01 ............. 

3 , 5371c * 01 . 

9 , 4438 c - oi . 

1 . 2534 E 00 

1 . 6740 E 00 ! .!. 

....2,413»E -00 *--*-*** 

3.4842F. 00 — 

4,64&7C 00 
5.7O20E 00 
. 6,3692E .00 — 

6 J 4731 E 00 - ..... , 

3.9639E 00 — 

4 ,sieaE go — :::::::::::::::: 

-3.5142E- 00 -*’*-**” --! ... 

-6 , 4975 E * 0 l — — — — 

:?#!!! 

— »l,6767E—00 





00 

00 


—1 ,4467! 

.1,2806_ 
t1,2041E 00 

--1.2103E S2 ♦ 

-i,2644e 00 . ♦ 

-1.3199E 00 : - 

-1 ;275bc-oo — ■ 

item. ? 

-6,6296E-0l 
-4,2413E*01 

-J:§?S?g:8J 

. 2*0972C*01 
2I3472E-S1 

2.2693E-01 ; ......... 

2.0380E-01 - ' 

i,ao24E-oi 

i,6596e-oi ■ — _ + 

i:739or-oi — — "’I";,..;*..* 

.. 1,6848c* 01 - * 











FIGURE 6-22 RELATIVE ANGULAR POSITION FIGURE 6-23 CONTROL TORQUE BODY 

BODY 4 . 




FIGURE 2-24 CONTROL TORQUE BODY 2 


62-9 


MNimm 

-3.Q637E 01 


VERSUS TIflE 


.3.0637E 01 
-1.1177E- 01 




1.1177E- U1 — *“+ 

4.7994E 00 """Ill” 

3i4633E 00 ; 


-1.5327E 00 
•3.4673E 00 
-1.918SE 00 




5.5000E 00 
6.0000E 00- 
6.3000E 00 


♦ 

• -+ 


B;4034E 00 


1:S§SIE 


OOOOE 00 * .7 'S066£ OO 


5.3841E 22 
5;2aa2E oo 
U3043E 00 


— ♦ 


3843E 00 — — 

4718E-0a 


1:18811 55 

l*7500E*”Ol 
1 ,S0QQE gl 
1.B500E 01 
1.9000E gl 
1.9500E gl 
2.0000E gl 
2.050GE 0l 
2^10001 01 
2.1500E gl 
2.2000E gl 
2,2500E 01 
2.3000E gl 
2.3500E 01 
2.4000E gl 
2.4500E 01 
2 ( S00QE 01 


-i;6282E-0l 
-1,3206E 00 
-i:4696£ 00 
_-l,4222£- gg 
.1.2709E 00 
-i;O0O7e gO 
-Bj9l05E*0l 
-•*7,21416-01- 

SZSttHl 

-3, 65446-01 
-.2.8813E-01 
-2.2530E-01 

:1:1SS8|:81 

-wins i- 

-5.4464E-C2 

t3J9450e*02 

-2,a326E”«2 

-2,02226-02 

-1.4391E-02 

zi-.nniqi 

-2.6149E-03 

♦1.8499E-03 

.i;3028E-03 

-9;i387E-04 










2*a326t-02 

2,8149E"03 1 

1.8499E-03 — •'" ’’''“"I"..;.:... -I 

i,3028e*o3 ■! 

6.3733E-04 + 

4.43876-04 - * 


3. The maximum absolute value of t 1z increases 
exponentially as the rigid body dynamics are moved 
farther into the left half plane. See Figure 6-25. 

4. The maximum absolute value of t 2Z is proportional 
to the peak to peak excursions in the state variables 
associated with the solar wings. From Table 6-2 the 
feedback gains increase significantly as the real part 

of the eigenvalues increase. Also the rate of increase 
in the gains associated with \p 2 ^ and is greater than 

the rate of increase in the gains associated with 
and as Re( X) increases. The significance of this 

will appear when sensor failure studies are undertaken. 

The above analysis relates to the time domain simulation 
of the Z axis reduced order model. Application of the 
control laws developed above were applied to the 
full scale five-body three-axis simulation. 

These results are provided in Figures 6-26 through 6-30 
for eigenvalue sets cl and b2 (corresponding to 
Figures 6-10 to 6-14 of the ROM). A comparison of the 
simulations for the full order model (FOM) with the 
simulations of the Z-axis reduced order model (ROM) 
demonstrates the effectiveness of the design method 
and the validity of the model reduction procedure. 

The simulation results for the FOM for the other eigenvalue 
sets obtained using the ROM control also showed excellent, 
comparison with the ROM simulation results. 
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6.4 CONCLUSIONS AND RECOMMENDATIONS 

An eigenvalue analysis of the five-body three-axis Space 
Construction Base/Configuration I vehicle shows that 
. the eigenvalue shifts from a double pole configuration 
are due to the interaxis coupling of the solar wings. 
Cross-axis coupling for this vehicle has been shown to 
be negligible thus leading to the conclusion that from 
a control design viewpoint, each axis can be designed 
separately. 

Using a symmetric reduced order model for the torsional 
axis, two new design procedures were studied. The multi- 
variable Nyquist array •' method, a frequency domain design 
technique, was shown to be effective for the difficult 
case of the flexible modes within the control system 
bandwidth. Utilization of the multivariable Bode dia- 
gram facilitated the design process since standard 
single-input single-output methods could be used. The 
Retallack pole placement technique when coupled with a 
computationally efficient numerical algorithm to obtain 
the Luenberger canonical form was effective in placing 
specific open loop poles in desired positions. Time 
domain simulations demonstrated the potential of this 
method. 

It is apparent from the results obtained that the methods 
used in this study should be investigated further. 
Specifically., several areas of consideration are identi- 
fied: 

a) Design and evaluation associated with pole shifts 
in frequency as well as in damping ratio. 
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b) Repeat design for the Y axis with and without the 
symmetry assumption. Compare designs using both 
MNA and pole placement procedures. 

c) Repeat the design for the X axis using MNA and pole 
placement . 

d) Sensitivity analysis of the control configuration 
under model parameter uncertainties. 

e) Sensor failure analysis in all axes and a simpli- 
fication in the system gain space description for 
failure accommodation. 

f) Develope a test for robustness of the control design 
• as a means of simplifying the design for future 

configuration studies. 

g) Incorporate design constraints on control magnitude 
and state measurement in the selection of control 
weighting vectors used in the pole placement procedure 

h) Incorporate gain magnitude constraints and develop a 
criteria for gain selection based upon time domain 
waveforms. 

i) Consider a Kalman pole placement algorithm to sequent- 
ially place poles "in desired locations using the 
criteria suggested above. 

j) Compare results using standard pole placement (Simon 
and Mitter) techniques and evaluate subsequent designs 
for robustness . 
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7.0 ALTERNATE CONTROL APPROACHES 

Among alternate approaches to the application of control 
to the Space Construction Base one of the most promising 
appears to be the general one presented in this section. 
In the sequel this approach is developed for a plant 
(system to be controlled) that may be represented by a 
■ broad class of linear, constant parameter multivariable 
mathematical models. 

i 

The general development' of, this approach begins with a 
description of the salient characteristics of the class 
of linear constant parameter multivariable plants for 
which the approach is applicable. The matrix properties 
of the plants are then utilized, in the development of 
techniques for the generation of state regulator control 
and estimation .gain matrices, that optimize plant per- 
formance . 

Following the development of the general approach to 
the design of state, controllers and estimators for 
this class of linear constant parameter multivariable 
plants., the considerations involved in applying this 
general approach to the rotational mathematical model 
of the Space Construction Base are discussed. These 
considerations include treatment 'of cyclical aero- 
dynamic disturbance torques, selection of weighting 
coefficients for state variable estimation error in 
the computation of optimal estimator gain and decom- 
position of the known portion of input disturbances in- 
to terms corresponding to gravity gradient, actuator 
reaction and actuator desaturation torques. 
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LINEAR, CONSTANT PARAMETER, MULTIVARIABLE CONTROL 
APPROACH 


7.1.1 


Plant 

Consider' the following plant 

x = A x +C( u+ v) +D w 
— p P—P — 

v=B x 

^ p-p 

where : 

x = plant state vector of dimension n 
“P 


(7-1) 

(7-2) 


u = plant control vector. of dimension c 

'v - plant control disturbance vector of dimension c 

w = plant disturbance vector of dimension, d 

y = plant measurement vector of dimension m 

A = uncontrolled plant 'matrix of dimension nxn 
P 

C = plant control matrix of dimension nxc 
D = plant disturbance matrix of dimension nxd 

Bp = plant measurement matrix of dimension mxn 


Suppose that the jth column of the plant control matrix 
is linearly dependent. Then we can write: 


C (u + v) = C j(Uj+Vj)‘+Cj(Uj+Vj) 

= C.(u,+v.)+C.c(u +v.) (7-3) 

«J' J ' J J , J J 

■ c j[VV ( VV£] - c '<*j + 2p 


where : 



C . 
3 

= C 

minus 

-j 

— u 

minus 

-j 

“ V 

.minus 


the jth column 
the jth component, 
the jth component 
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c. = jth column of C 

J 

u. = jth component of u 

«J 

v. = jth component of v 

J 

c = constant vector of dimension c-1 


Thus we can assume, without loss of generality, that C 
has rank c. If M is a real .symmetric matrix, then M and 
M^M- have the "same rank. 'Hence,' defining: 

M *=(M T M‘)~^M T ( 7-4) 

we can write: 

C ( u+v ) +Dw=C ( u+ v+C*Dw ) + ( I • - CC* ) Dw 

- C(u+v")+D''w (7-5) 


T ^ 

Since CD ;=0, we can assume without loss of generality 
that the columns of C are orthogonal to the columns of 
the plant disturbance matrix. Now suppose that the ith 
row of the plant measurement matrix is linearly depen- 
dent. Then we can write: 



(7-6) 



(7-7) 


where : 

B' = B minus the ith row 
p P 

y-' = y minus the ith component 

y^ = The ith component of y. 

c = a constant vector of dimension m-1 
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Thus we can assume without loss of generality that 
has rank m. Furthermore, it is assumed that: 


1) c < n and m ■< n 

2) All matrices are constant 

3) The unknown part of v and w can be written as 
B d x d where x d satisfies the following matrix 
differential equation. 


A <#d 


(7-8) 


That is , 


v = 


v. 


+ B x 
V— V 


(7-9) 


where: 


,x = A x 
—v v— v 


and 


w = w, + B x 
— — k w— w 


(7-10) 


(7-11.) . 


where : 


x 

— w 


A x 
w— w 


' (7-12) 


7.1.2 Design of State Feedback Controller 

We wish to design a state feedback control system for 
this plant so that the plant state vector will remain 
near some equilibrium plant state vector (Xp) d - 

Since , A p (x p ) d =0, (7-13) 

we can write - : 

x" = Ax" + C(u+v)+Dw (7-14) 

P P P 
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(7-15) 


y" = B x" 
“ P- P 


where : 


-'p = V% 5 d 


(7-16) 


and 


v' ~ y-B (x ) , 
^ - p -p'd 


(7-17) 


Hence, we can assume without loss of generality that: 


( Vd " 0 


'7-18) 


Since 

C T D = 0/ ^7-19 ) 

none of the plant disturbance vector can be cancelled 
by the plant control vector. However, all of the plant 
control disturbance vector can be cancelled. Therefore, 
we let : 

u = g(x p )-v (7-20) 

To simplify the design, we will limit the feedback 
structure to a linear, time-invariant function so that 

u = ~G5£ p “v' (7-21) 

where G is a controller gain matrix of dimension c x n. 
Substituting equation (7-21) into equation (7-1) yields: 

x = A x +Dw' (7-22) 

— p c p — 

where the controlled plant matrix may be expressed: 

A = A - CG (7-23) 

c p v ' 
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Note that if we could choose G such that 


A c “ 


I, 


(7-24 


where I = identity matrix of the same dimensions as A , 
then the present plant state vector would he 


x (t) = 
~P 


e T x(0)+D 


“P 


r 


T-t 

e T w(x)dT 


(7-25) 


where x p (0) is the initial plant state vector. This 
would be a desirable choice for G since, under the 
assumption that T is small, only recent, large plant 
disturbances could cause the present plant state not to 
be small. However, since c is less than n, there 
generally is no matrix G that will satisfy equation (7-24). 
Therefore we will choose G to render a measure of the difference 
between the two matrices in equation (7-24) as small as 
possible . 


We could define the measure of the difference between 

these two matrices in the same way as for the measure of 

T i 

the difference between two vectors as (trE E ) where: 

E c = A c + ^1 (7-26) 

This yields the following solution. 

G = C*(A + Jl) (7-27) 

(The derivation and analysis of this solution appear in 
Appendix C) 


However, we wish to make E small in the sense that for 

c 

any x p , E c — p is sma -*-l compared with x . Furthermore, 
for various physical reasons, some plant state variables 
or combinations of plant state variables may be more 
important to control than others. Consequently we will 
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choose G to minimize the following quantity: 


max 


W. E.z 
c c— 


zfo | W^z 


(7-28) 


where the columns of W are linearly independent. 

c 


T 

Since W W is nonsingu-lar , we can write: 
c o 


max 

z^o 


l W c E c^l 

W z 
c— 


max 

z^o 



where 


E ■*= W E VT 
c c c c 


z = W z 
— c— 


(7-29) 


(7-30) 

(7-31) 


Dividing the numerator and denominator of the right hand 
side by |z"j and then squaring, we see that the above 
problem is equivalent to choosing G to minimize the 
quantity : 

max z T (E ") T E "z ^7-32) 

— c c — ^ 

|zi=l 


Introducing the constraint, 

\z\ - 1, ' (7-33) 

by the Lagrange multiplier approach, a necessary con- 
dition for a critical value of z is.: 

9l[- T(E c )T V- -MhV 1 )]^ (7-34) 
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where 


3_ 

dz 


indicates a gradient with respect to z. 


Completing the dif ferentiation indicated on the left 
hand side and dividing by 2 yields: 


f(E ") T E ^ -Allz=0 (7-35) 

i. C — C ^ jji 

Premultiplying both sides by' z and rearranging terms 
yields: 


X = z T (E ") T E" z (7-36) 

— c c — 

From equations (7-35) and (7-36), we see that: 

T T T 

max z (E •*) E 'z = largest eigenvalue of (E ') E ' 
c c c c 

121=1 (7-37) 


Hence, we wish to choose G such that A, the largest 


eigenvalue of (E ') E ' 

c c 

v- w cg + 


is a minimum where: 
* 


(7-38) 


Since G is unconstrained, a necessary condition for the 
minimum value of A is 


3 A 
3G 


= 0 


(7-39) 


where — indicates the gradient with respect to G. 

9G 

T 

Furthermore, since (E ") E ' is symmetric, A can be 
approximated by the fact that: A = lim | — k+i | ( 7.40) 


k->°° r 




where: 

e 


T =(E ') E "e, 
k+1 c c — k 


(7-4 T) 


(Appendix D gives a numerical algorithm for generating 
the optimum controller gain matrix) . 
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7.1.3 


Design of State Estimator 

The plant control given by equation (7-21) requires com- 
plete knowledge of the present plant state and control 
disturbance. Generally, these quantities are not com- 
pletely known. However, they can be estimated from 
the known part of the past plant disturbances and the 
past plant measurements. Let x g be defined as the solu- 
tion to the following vector differential equation: 

x = (A-HB)x +[d T !o!o1 T w +Hy (7-42) 

— e — e i i t J — k 

where : 


A= 


A 1 0 ! DB 
c | ^ w 

0 I 0 I A 
! I W . 


(7-43) 


H = constant estimator gain matrix 
B = [ B p | 0 j ° ] (7-44) 

x e (0) = B T (B p B p T )" 1 Z (0) (7-45) 

Combining the assumptions about the system disturbances 

with equations (7-2) and (7-22) yields: 

• T T 1 i IT 

X = Ax + [D i 0 i 0 j w k (7-46) 

y = Bx (7-47) 

where : 

K|Sw] T < 7 - 48 > 

Subtracting equation (7-46) from equation (7-42) and 

substituting equation (7-47) for y, it is evident that: 
Ax = (A-HB) Ax (7-49) 


7-9 



where : 


Ax = x 0 -x 

Therefore, if we choose H to satisfy: 


(7-50) 


A-HB 



(7-51) 


then x would be a good estimate of x since the solution 

— e 

to equation (7-49)_miuld be 

-20t (7-52) 

Ax(t)=Ax(0)e T 

However, since m£n, there exists no H such that equa- 
tion (7-51) is satisfied. Furthermore, the estimator 

T >4 

gain matrix that minimizes (trE E )® 


where: E = A-HB + isr I . (7-53) 

e I 

yields estimates for the unknown part of the control and 
plant disturbances of zero at all times. (See Appendix 
E for a derivation and analysis of this solution). 


Consequently, we will choose H to minimize the largest 

T 

eigenvalue of the matrix, (E ') E " 


where: 


E W E W 
e e e e 


(7-54) 


(See Appendix D for a numerical algorithm for generating 
the optimum’ estimator gain matrix.) Figure 7-1 depicts 
the relationships between the plant, state controller and 
state estimator in the complete system. 
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DIGITAL I 
COMPUTER' ) 


FIGURE 7-1 

BLOCK DIAGRAM OF COMPLETE' SYSTEM 
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7.1.4 


System Design Parameters 

Since the eigenvalues of equation (7-42) are approximately 

- the sampling frequency of the digital computer 

■L 400 

should be approximately — . The matrices C and yield 
trade offs between control system performance and cost . 

The number of linearly independent actuators and sensors 
and, consequently, the system performance, increases as 
the ranks of these matrices increase. However, this 
better control system performance may not be worth the 
increased control system cost. The scalar T also affords 
a trade off between control system performance and cost . 
This trade off is more complex than those for C and B p 
since smaller values of T not only increase control 
system cost but also may decrease control system per- 
formance. A control system simulation and trial and 
error will be necessary to choose the best values for 
C , Bp and T . 


If a satisfactory control system cannot be found, it 
may be possible to significantly increase the performance 
and/or decrease the cost of the control system by making 
some relatively minor changes in the design of the un- 
controlled plant. For instance, assume that 



where f . . (rr) is the ith row and jth column element of 
— J 

Ap and i is a vector of plant parameters which can be 
chosen to make control easier . 


If we could render: 



(7-56) 
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7.2 


then G=0 and the control system would be optimal . How- 
ever, since generally there will he no tt such that these 
matrices are equal, we will render the difference between 
■the two matrices of equation .(7-56) as small as possible 
by choosing tt to minimize the largest eigenvalue of 

(E ") T E ** where: 

P P 



(7-57) 


APPLICATION TO SPACE BASE 

A major problem in designing control systems, for flexible 
spacecraft is the development of an accurate yet simple 
plant model. Because of the inherent difficulties of 
solving partial differential equations a finite element 
approach must be used. Space Base can be modeled as a 
number of rigid bodies connected in' a tree configura- 
tion by either dissipative spring-hinges .or electro- 
magnetic torque actuators. The greater the number of 
rigid bodies used in the model , the better the accuracy 
of the model and hence the better the performance of 
the control system. However, if more rigid bodies are 
used than are necessary, then the order of the model 
and, consequently, the cost of the control system will 
be greater than necessary. It is conjectured that, for 

the control scheme of subsection 7.1, if the lowest 

20 

resonant frequency of a physical structure is or 
higher, then that structure can be considered rigid. 


Assuming that the control system is holding or very 
slowly changing the orientation of Space Base, the 
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equations of motion of the rigid body finite element 
model can be linearized and written in the form of 
equations (7-1) and (7-2) with the angular rate and 
position error of the most precisely pointed body and the 
relative angular rates and position of adjacent bodies as 
state variables. 


7.2.1 Weighting of State Variables 

There are two obvious reasons to weight the individual 
plant state variables when computing the optimum con- 
troller gain matrix. First, it makes no sense to 
add quantities of different physical dimensions such as 
angular rate and position. Second, it may be more 
important to control one section of Space Base than 
another. Consequently, if 


x 


= [e T ; e T ' 


(7-58) 


where the vector 9 contains the angular rate plant 
state variables and the vector 9^ contains the angular 
position plant state variables, then a logical choice 


for W would be: 
c 



TW { 0 

0 ] W 


(7-59) 


where the scalar, T, nondimensionalizes the angular 
rate plant state variables and W is a positive definite 
diagonal matrix which weights the relative importance 
of controlling the individual angular position plant 
state variables. For this choice of W , we have: 
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* T — 1 T 

W = (W A W ) X W 
c c c c 


— W ^ 


0 


I 

I w 


(7-60) 


Furthermore, a better choice for the starting value of 
G in the algorithm of Appendix D might be: 


K - (W c C)*.W c (A p+ il) 


(7-61) 


7-. 2.2 Aerodynamic Disturbance Torque 

Assuming the earth's gravity gradient disturbance torque 
to be computed and thus known, the primary unknown dis- 
turbance torque acting on Space Base is aerodynamic. 
Since this torque is approximately sinusoidal with orbit 
frequency, we can write: 


v u = a. cos w Q t + b sin w t (7-62) 

where v u is the unknown part of the control disturbance 
torque and a and b are constant but unknown vectors.. 
Differentiating v^ twice with respect to time yields: 


■ 


* 

1 


- 

v 

— u 


0 

i 2 

1 “o 1 


• 

• ^u 

• 

V 

— u 


I 

1 

1 o 


v 

— u 


(7-63) 


Accordingly, 

■ [i 


T ' 

i 

l l 


B 

v 


[°h 


TIT 


T] 
V - 

-u J 


(7-64) 
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Similarly , 


r. T 1 T]T 

X = W I w 
— w L— u I — u 

b w = [° N 


(7-67) 

(7-68) 

(7-69) 


7.2.3 Weighting Coefficients For Disturbance State Variable 
Estimation Error 

It is reasonable in computing the optimum estimator gain 
matrix to weight various combinations of the disturbance 
state variable estimation errors according to the magni- 
tude and location of the jerk, acceleration and velocity 
they cause. Thus a logical choice for W is: 
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w = 

e 


W 


0 


0 


0 


0 


r T s 

T 
I 

-r 

l 

-l. 


I T W c C 


0 


I 0 

-I 

! o 
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0 


0 


TW c C 


0 




1 0 | 0 

I 0 I 0 

_± 4 

_l 0 * 0 

i ~ ~2 i 

T w c D 0 
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i 


TW D 

C ' 
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For this choice of W , we have 

e 


w: 


* 

0 

J_ 0 

i -C)* 

j T Z C 

! 

! 

o j o 
! 

i 

1 

! 

o | o 
1 

1 0 

r°~~ 

0 

1 0 
1 

J ^(W c C)* J o 

1 

1 o 
1 

•1 

0 

1 0 

i 

! 0 

?-< W c D > 

*| 0 

0 

1 0 

1 0 

. 0 



(7-71) 
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7.2.4 


Known Input and Plant Disturbances 
Assuming the use of rigid body momentum exchange 
actuators , the known part of the input disturbance can 
be divided into three components: gravity gradient, 
actuator reaction and actuator dosaturation as follows. 

v = v + v + v. (7-72) 

-k -g -r -d . 

The known part of the plant disturbance is just gravity 
gradient . 
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8.-0 STRUCTURAL ANALYSIS STUDY 

The determination of flexibility body characteristics 
for a series of connected substructures is an important 
problem in the study of spacecraft dynamics. Assuming 
that the flexibility characteristics of each of the 
substructures comprising a vehicle system are known 
in terms of the free-free modes which are obtained 
from a discrete coordinate formulation of the individual 
substructures, a method (Reference 8-1) has been 
developed recently for combining the substructures by 
deriving the forces and torques - that exist at the 
connection points between contiguous bodies while 
satisfying the appropriate boundary conditions. It is 
the objective of this section to determine the validity 
of the above proposed method by considering a homogeneous 
beam as an example. 

8.1 SUMMARY OF APPROACH 

A homogeneous beam made of steel is considered as the 
example to investigate the validity of the method 
mentioned above. Initially, the- beam is idealized to 
be composed of elastically interconnected discrete 
rigid bodies with no damping mechanism involved. The 
springs connecting the bodies are assumed to be capable 
of having both 1 , -translational an d rotational motion. 

The free-free modes of the beam are determined from the 
equations of motion of the discrete rigid bodies 
constituting the beam. Next, the beam is halved and 
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each half of the beam is modelled in a similar' manner 
as the total beam. Hence, the free-free modes of each 
half of the beam can be determined. Then the two 
half-beams are connected by two connection points; and, 
by writing the forces and torques that exist at the 
connection points, the equations of motion of the 
connected beam are developed to obtain the free-free 
modes of the connected beam. This result is compared 
to the modes of the beam obtained previously. 

8.2 TECHNICAL DISCUSSION 

The beam is subdivided into 2n rigid bodies. The 
adjacent rigid bodies are connected by springs which 
are capable of both translational and rotational motion. 
Thus the idealized model of the beam, illustrated in 
Figure 8-1, is composed of elastically interconnected 
discrete rigid bodies with no damping mechanism . involved. 
Figure 8-2 and Table 8-1 explain some of the nomenclature 
used in this Section. 

It is assumed that the discrete bodies composing the 
beam undergo deformat ions . so small that only terms 
involving the first-degree deformations need be retained 
in the analysis. Further, the beam is allowed unrestri- 
cted .motion in response to any applied forces and 
torques. Before proceeding with the analyses, it is 
necessary to define the reference frame with respect 
to which all measurements are made. If an inert ially 
fixed reference frame is chosen, the deformations 
relative to the frame may grow large if the beam 
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FIGURE 8-1 IDEALIZED MODEL OF THE BEAM 



FIGURE 8-2 ILLUSTRATION OF SOME NOMENCLATURE 
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TABLE 8-1 


NOMENCLATURE 


vector from inertial point to CM of the beam 

vector from inertial point to instantaneous 
(inst.) CM of body i 

vector from CM of body i to instantaneous 
(inst.) CM of body i 

vector from inst. CM of body i to an elemental 
mass in body i 

vector from CM of body i to the point of appli- 
cation of external force 

inertial rate of the rigid body frame of the 
total beam 

inertial rate of the rigid body frame of the 
half beam #£ 

angular rate of body i relative to the beam 

vector from the connection point to the CM of 
bodjr i 

vector from the CM of body i to the interface 

point that lies on the line joining the CM’s 

of bodies i and j . 

interface force between adjoining bodies i and 

constraint force between adjoining bodies i 
and j 

external force applied to body i 

angular momentum of body i about the CM of the 
beam 



TABLE 8-1 (CONT'D) 


J i 


inertia of body i about its CM 


K, K„ 


K e, K e 1 




m . 

i 

M. . 


M. 


C id 


M n 


M n 


H- 


n 


Q o,V 


£ r i 


R i 


T 

e 


i 


translational stiffness matrix between adjoin- 
ing bodies 

rotational stiffness matrix between adjoining 
bodies 

refers to half-beam #1 for Z=1 and to half 
beam #2 for &=2 

mass of body i 

interface torque at the adjoining bodies of i 
and j 

constraint torque at the adjoining bodies of 
i and j 

total mass of the beam 
mass of the half beam #£ 


the number of -bodies constituting a half beam, 
or one-half of the number of bodies composing 
the beam 

vector of the deformation coordinates 

vector from the CM of the beam to the CM of 
body i 

vector from the CM of half beam #£ to the CM 
of body i in the half beam 

vector from the CM of the beam to inst. CM 
of body i 

external torque applied to body i 
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undergoes an appreciable motion so that the analysis, 
using the assumption that’ the deformations are small, is 
not valid. However, if the reference frame is chosen 
such that it moves with the beam in a well-defined 
manner, then the deformations relative to this "floating” 
frame may be assumed to be small and hence the '!first 
order deformations only" analysis can be carried out. 

The Tisserand frame (reference 8-2), with attractive 
properties for the purposes of this report ,• presents 
itself as such a floating reference frame. It is defined 
by the set of axes that moves in such a way as to set 
the internal angular momentum and the internal linear 
momentum always to zero. The latter requirement 
implies that the origin of the frame be located at- the 
(rigid body) center of mass. 

8.2.1 Preliminaries 

The translational and rotational equation of motion 
'for any discrete body i are given by: 

m i a ± (8-1) 

H. (8-2) 

where F i and T i are, respectively, the force and the 
torque exerted on body i,- m^ is the mass of body r, 
a^ the inertial acceleration of the center of mass (,CM.) 
of body i, and the angular momentum of body i 
referred to its CM. The small circle (o) indicates time 
differentiation with respect to inertial space. 
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The following equations are of use in the ensuing dev- 
elopment of the equations of motion for the bodies of the 
beam. By the definition of the CM of the beam 


2n 

£ ' m i r i = 0 (8-3) 

i=l 


and by the translational equation of motion of the beam 


oo 


p 


o 



(8-4) 


where the total 


external force is given by 


F = 
e 


2n 

E 

i=l 



(8-5) 


and the total mass of the beam is given by 


m t 


2n 

Z 

i-1 


m. 


( 8 - 6 ) 
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The rotational equation of motion of the beam is given 
by: 


* 

J 


o 

« 

0 



(8-7) 


where the total external force is given by 
2n 

T = £ T (8-8) 

e i-i e t 


and J is the inertia matrix about the CM of the beam. 

* 

Since the beam is deformable, J is not a constant. 

Equations (8-4) and (8-7) define the linear and angular 
accelerations of the Tisserand frame with respect to 
which all measurements are made. 

8.2.2 Constraint Forces and Constraint Torques 

Since the beam is modeled as elastically interconnected 
discrete rigid bodies, there exist constraint forces 
and constraint torques that are impressed on the discrete 
bodies due to the deformations of interconnected bodies.. 
10 derive the expressions for these forces and torques, 
the interconnection of a typical discrete body i with 
bodies i+2 and i+1 is illustrated respectively in 
Figures 8-3 and 8-4. 


8 



K/4 

FIGURE 8-3 CONNECTION OF BODIES i AND i+2 


Z 



FIGURE 8-4 CONNECTION OF BODIES i AND i+1 
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Body i is connected to body i+2 (or i-2) via four 

springs, each with a translational spring constant jK 

1 ^ 

and a rotational spring constant ^-K^. Similarly, it is 

connected to body i+1 via four springs each with 

translational and rotational spring constants 4 k and 
1 ^ x 
4 K 0' n ’ 


Let denote the vector from the CM of body i to the 
vertex j' of body i. Then the constraint force that 
is exerted on body i due to body i+1 can be written as 

,1 T 


F 


c . 

i, 


i+1 


. E ‘ "A^i+l^i^ 
J=1 ' „ 

4 


+ ill jKj <B 1+1 x dj +1 - Si x dj) (8-9)- 


Since 


4 i 
2 d. 

1=1 J 


4d i,.i+l’ 


2 d^ +1 = 4d. ,, . 

j=l J 1+1 


( 8 - 10 ) 


d . . , , — d . , , . 

i,i+l i+l,i 
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and with the skew symmetric matrix A defined by 


A 




o 




( 8 - 11 ) 


where A x , A , A z are the x, y, z components of the vector 
A, equation (8-9) simplifies to 


i+1 


= Ki(y 1+ i-y 1 ) + K 1 d ljl+1 (6 I+1+ e i ) (8-12) 


Similarly, the constraint force impressed on body i+1 
due to body i is given by (8-13) 

* c - K 1 (u r w l+l )+K l d iH-l,i (B i +S i+l) 

i+1, i 
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In a similar manner, the constraint forces between bodies 
i and i+2 can be obtained; to" be 

F =K(y -Vi )+Kd (0 +0 ) (8-14) 

i , i+2 i+2 i i;if2 i+2 . i 


? c. , . = K u i+2 ) + K d 1+2 i <e ± + e 1+2 ) 

1+2 , i 


(8-15) 


It can be seen from the above equations that 


r p 

c . . — c . . 

Ji 


The constraint torque that is experienced by body i due 
to body i+1 can be written as 


M. 


'i , i+1 


E d 1 x | K ± (u ±+1 - W± ) + 

j=l 3 


..(8-16) 


4 

E 

j=l 


dj *1% (B 1+1 x df 1 - B ± X dj) 
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With d\ as the vector in the xy plane of the body i 
from the center of elasticity of the four springs conn- 
ecting bodies i and i+1 (the center is the point of 
meeting of the vectors d i i+1 and d ±+1 ± ) to the j th 
spring, equation (8-16) using (8-10) can be written as 


- d i,i + l y C V i - y + 

■X 

2 

j=l 

1 | 

( d i >i+1 

+ d 1 .) x 
oj' 

K l + ~>i - 

4 ■ 
2 

j = l 


^ d i,i+l 

+ d 1 • ) x 
oj ' 


. K 1 «i + l,i + O'] B i + 1 < 8 - 17 > 


The coefficient of 6 ± in equation (8-17) the term in 

brackets, simplifies to (noting that d .= d . ) 

oj oj 


d. .... K- d. 
i,i+l 1 i,i+l 


4 ~ 
2 

j=l 


2 d 


l 

oj 


ri i+1 K 
d oJ K ] 
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The second term in the above expression may be recognized 
as the rotational spring constant between bodies i and 
i+1 and hence is replaced by •' Thus equation (8-17) 

can be written as 


M 


c . 


i , i+1 


" d i,i+l K 1 (y i+l ' w i ) + 

d i , i+1 K 1 d i, i+1 ^i+l + + 

K 0, < 6 i«l - 


(8-18) 


Similarly, the constraint torque on body i+1 due to 
body i is given by 


M 


c i+l, i 


d i+l,i K 1 <lI i y i+l ) + 
d i+l,i' K l d i+l,i (B l + B i+1 ) + 

K n <»i " W 

.1 1 * 


(8-19) 
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In a similar way, the constraint torques between bodies 
i and i+2 may be obtained to be 


M 

c i , i+2 


d i,i+2 K Cv i+2 v i^ + 
d i,i+2 K d i , i+2 ^ 6 i+2 + 6 i^ + 


^0 ^i+2 ~ ^i^ 


( 8 - 20 ) 


M 

c i+2 , i 


d . K (y. - y.' ) + 

r+2 ,x l 1+2 J 

d. 0 . K d. 0 . (8. + 8 ■ , o ) + 

i+2,i i+2,i v l i+2 ' 


K 0 ^i ^i+2 ) 


( 8 - 21 ) 


8.2.3 Translational Equations of Motion for the Full Beam 

The translational equation for a typical body i can be 
written, using equation (8-1) and recognizing that the 
constraint force between the adjacent bodies contributes 
to the total force impressed on body i, as 
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+ F„ + F + F 

6 i °i, i-2 c i,i+l c i,i+2 


oo 
m . p . 
I X 


i = 3,5, ..... . ,2n-3 


( 8 - 22 ) 


From Figure 8-2, it follows that 


P = P + Rj 
i o 


R. = r. + v . 
x x K x 


(8-23) 

(8-24) 


i = 1 , 2, . . . , 2n 


Combining (8 r 23) and (8-24); and differentiating twice 
in the inertial frame results in 


OO oo OO OO 

p i = p o + r i + Vi 


(8-25) 


The terms appearing in the right-hand side of the above 

equation can be written in terms of body i rigid body 

coordinates : 

' **• «• 

oo 

r -i = % x r -i + x x ~ w x r. 
x o x o o x o x. - 
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00 ’* * * *• 

y i =li i +a) o xu i +2u o xu i +w o Xw o xy :L ^y i 

where the dot («) indicates time differentiation in a 
reference frame located at the CM of body i. The 
approximation results due to the assumption that the 
deformations are small, which entails the deletion of 
nonlinear terms. Thus, with the above expressions and 
the equation (8-4), (8-25) becomes (writing the results 
in body i rigid body coordinates): 

oo f • , »• 

P ± = ® + % x r ± + (8-26) 

M T 

Substituting for °p\ from (8-26) and for the constraint 
forces from section 8.2.2, the translational equation 
of motion for body i can be obtained as 


m 


i U 1 +(2K + K 1 ) Vi -K Wl _ 2 -K 1 P 1+1 -Kp i+2 -K 1 d i;1+1 e 1 +K4 l _ 2 _ i 8 1 . 2 + 

K 1 d i+l,± K d i+2- i ®i+2 = - • 

m . ~ 

i * 

P — x-, — P p m . X- <0 

e i e 1 1 ° 


i=3 , 5 , ..... ,2n-3 


(8-27) 


Similarly, the translational equations of motion for 
the remaining bodies can be obtained. A summary of 
the translational equations is given in Appendix F.l. 
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8.2.4 Rotational Equations of Motion for, the Full Beam 


Using equation (8-2) and noting that the constraint 
torques,, which include the effect of constraint forces 
between the adjacent bodies, contribute to the total 
torque applied to body i (again, a typical body in the 
top row of the beam model), the rotational equation 
of motion for body i can be written as 


T +Y . xF +M +M 

e i 1 e i c i , i-2 c 


+M- 


'± , i+1 




(8-28) 


i=3 j 5, . . . ,2n-3 


The rate of change of the angular momentum of body i 
is given by 


o • 

H = T..H, 
10 11 


(8-29) 


with 


H± = J-, 0 T* (S 0 + B 1 )+J 1 [ 0 T|(i 0 +B 1 ) Xo Tj«; 0i +5 1 )] (8-30) 


where T. -is a fixed transformation from the rigid 
o i 

body coordinate frame of body i to that of the inertial 
(Tisserand) frame. Neglecting nonlinear terms in (8-30) 
and noting that the transformation matrix c'an be approx- 
imated by the identity matrix, equations (8-29) and 
(8-30) combine to yield 


8-18 - 



(8-31) 


o 

H ± ^((0 


o 


*» 




Substituting for the constraint torques from section 3.2, 
equations (8-28) and (8-31) with the following defini- 
tions . 


n 


K *- d i,x±2 Kd i,i ±2 


(8-32) 


K. • , -,Kd. 

^ ■ 1,1+1 i,i±l 


(8-33) 


for all i such that the subscripts are positive, yield 
J i ® i + ( 2K e +K e i ^ e i“ K e g i-2 _K e 1 B i+l _K 0 e i+2 


+(2 V% > 1 » + % B i~2 + % g i+ 1 + S e i+2 
11 ! 


+d i , i+l K l y i d i , i-2 Ku i-2~ d i , i+l K l vl i+l _d i , i+2 Ky i+2 


T + y, p - J. M 
e^ 1 1 e^ i o 


i=3 ,5, , 2n-3 


(8-34) 
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The rotational equations for bodies i=4,6, . . . . . ,2n-4 
can be written in a similar manner. To illustrate the 
nature of equations for the end bodies (i=l , 2 , 2n-l , 2n ) , 
the rotational equation of motion for body i is given: 


j 1 e 1 +(K 0 +K 0i )3 1 k 0 _^b 2 k q b 3 


+( W 6 1 + V2 + V 3 


(8-35) 


+ (dL L2 K L +d 13 K ) yi- 4^ ^ - d~ 13 K y 3 


= T + y 


1 F e, 


- J, 


to 


o 


The rotational equations for i=2, 2n-l and 2n can be 
written in a similar manner. A summary of the rotational 
equations is given in Appendix F.2. 
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8.2.5 


Matrix Equation of Motion for the Full Beam 


In the equations of motion developed in the previous 

two sections, the right-hand sides (rhs) are functions 

of only external forces and torques, except seemingly 

for a) . The rate of change of the inertial rate of the 
0 

rigid body frame of the beam, w o _, may be defined as the ’ 
product of the inverse of the inertia dyadic of the 
discrete bodies about the CM of the beam and the total 
external torque impressed on the beam; i.e., 


tu. 



(8-36) 


where : 


J = 
o 


2n 

2 

i=l 


(J.-m.r .r . ) 
v x l x x' 


(8-37) 


and 


T = 
o 


I 


T_ + y J 


e i 


i e. 


(8-38) 


With the above definition (8-36), the rhs of the equations 
of motion is now only a function of external forces 
and torques. 


^--3 
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The translational and rotational equations of motion of 
the discrete bodies constituting the beam may be combined 
into a single matrix equation to facilitate computations: 


M C Q + K c Q = L q 


(8-39) 


where M and K Q are, respectively, mass and stiffness 
matrices of the beam and are 12n x 12n, L 0 is a 12n-column 
vector of external forces and torques, and Q is a 12n-column 
vector that characterizes the flexible beam deformations 
and is defined by 



fy t pt y t gt 

L 1 1 2 2 


yt pt *1 
2n 2n J 


where y^ ' s and 8^ ' s are 3-row vectors and the superscript 
t denotes the transpose of the matrix or vector. 


The beam model that is used in the computer simulation 
consists of 20(n=10) identical bodies so that 

= m 
J. = J 

for i=l,2, . . . ,2n. Using the equations of motion developed 
in sections 8.2.3 and 8.2.4, the mass and stiffness 
matrices M Q and K 0 can be defined as follows. 
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Mass Matrix 


M 

o L ij J 
m. . : 3x3 matrix 

0. : 3x3 null matrix 

r m i=j=l,3, , 4n-l 

=< J ±=j-2 , 4 , 4n 

Lo 

Stiffness Matrix K 0 = 

K^-j : 3x3 matrix 

O' : 3 x 3 null matrix 




K+K x 

i=j=l, 3, 4n-3,4n-l 

VVVS 

i = 3 "2 , 4 j 4n-2 , 4n 

2K+K 1 

i=j = 5, 7 , 4n-5 

2K e+K +2VV 

i=j-6,8 ,4n- 


K. .=K*. 
ij Ji 


~K 

- K 1 

%" K 0 

Sr K e- 


(i,j)“(l,5),(3,7), . ... , (4n-5,4n-l) 
(i,j)-(l,3), (5,7), . . . . , (4n-3 , 4n-l) 

(i,j)-(2,6) , (4,8), ,(4n-4,4n) 

( i, j )=(2, 4) , (6, 8) , ,(4n-2,4n) 


/ 


K. -=rf. 
1J 31 


= 


-K 

-K 


-K 


-K 


-K, 


+K 


(i,j)=(l,6), (3,8), , (4n-5 3 4n) 

(i, j)=(l,4),(5,8), , (4n-3,4n) 

(i,j)=(2,5),(4,7), . ...,(4n-4,4n-l) 

(i,j)=(2,3),(6,7), , (4n-2,4n~l) 

(i, j)=(5,6),(9,10), . . . , (4n-7,4n-6) 

(i, j)=(7,8),(ll,12), . . , (4n-5,4n-4) 
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K ir K ji 


- 


, -K„ -K a 

( 0 

" K P +K 0 1 

K 0 “ K 0 1 
V K 0 +K 0 X 


(i } j )“(3,4) 

(i, j )=(4n-3,4n-2) 
(i, j)=(4n-l,4n) 


K. . = 0 otherwise. 


The following definitions have been used in the above 
characterization of the matrix K 0 : 


V, i d i,i + i K i 


K 0 - d i,i + 2 K 


K 0 X d i,l+l 

% - K 0 d i , 1+2 


Also note that 


and 


d. . = -d . . 
ij J 1 


d i,i+2 d i,i-2 
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8 . 2.6 


Split Beam - Interface Forces and Torques 
The beam is divided into two halves so that each half 
of the beam may be modeled to be composed of n(=10) 
elastically interconnected discrete rigid bodies. The 
two halves of the beam are put together, as shown in 
Figure 8-5, by connecting the discrete bodies n-1 and 
n+1 and the bodies n and n+2. The beam thus connected 
is referred to as the 'split beam'. Since the two 
halves are modeled in a manner similar to that of the 
full beam, the free-free modes of the half-beams can be 
determined in a similar way as the full beam. Then the 
free-free modes of the split-beam can be obtained from 
the free-free modes of the half-beams and the forces 
and torques acting at the connection points of the half- 
beams. For computational simplicity, the above two 
steps can be combined to accomplish the same goal; 
namely, obtaining the free-free modes of the split-beam 
from the equations of motion of the 2n discrete bodies 
that compose the split-beam, keeping in mind that it 
is the interface forces and torques, not the constraint 
forces and torques, that act at connection points. 

Half-Beam #1 Half-Beam #2 



FIGURE 8-5 CONNECTION OF THE HALF-BEAMS 
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The equations of motion for the discrete bodies not 
involved in connections can be derived in an identical 
manner as for those 'of the fu-1-1 -beam“ In order to 
derive the equations of motion for the bodies involved 
in the connections of the half-beams, it is first necessary 
to determine the interface forces and torques acting at 
each connection point. 


The interface force at the connection point between 
bodies n and n+2 is determined by writing the transla- 
tional equations of the bodies n and n+2: 



+ F , +F 

c n ,n~l °n, n-2 


+ F n,n+2 


oo 

mp 


n 


(8-40) 


F 

e n+2 


+ F + F 

c n+2 , n+1 °n+2,n+4 


+ F n+2,n 


oo 

mp n+2 


(8-41) 


where F. . denotes the interface force at the connection 

point between bodies i and j . The bodies n and n+2 are 

connected such that they move as a rigid body so that 

F = p 0 Hence, the translational equation for 

n, n+2 n+2 , n 

this composite body can be written as 


F +F +F +F +F 

e_ e_ . „ c 


n ^n+2 v "n, n-l +F<3 n, n-2 ±c n+2,n+l - c n+2,n+4 2m p n-,n+2 

(8-42) 

where p n n+2 is the vector from the inertial point 
to the connection point. Further, 


oo 


p = p 

K n ^ 


+ e 


(8-43) 


n , n+2 


n 


n+2 


p n+2 p n,n+2 + e 


(8-44) 



From equations (8—40), (8—42) and (8—43), and noting 
that (neglecting nonlinear terms) 


oo • 


= (“oi + S n ) x e n 

where w 0l is the in ertial rate of the rigid "body frame of 
the half-beam #1, the interface force can be solved 
to be j- 

F , „= -me -me w -+£ 1 F +F C 

n,n+2 n n n ol \ c n+2 ,n+4 c n+2,n+l 


'n,n-l n , n-2 


+ £ (F 


e n+2 


^ ) 


(8-45) 


Using equations (8.41), (8-42) and (8-44), 

F n+2 , n = - m5 n+2 6 n+2 +4 | I 'c n _ n . 2 + F c n>n _ 1 + \ +2 , n+1 +1?0 n+2 , n+4 


~ me n+2“o2 + 2 (F e n F e n+2 > 


(8-46) 


where w q2 is the inertial rate of the rigid body frame 
of half beam #2 . 


Since the two halves of the beam have the same inertial 


rate; - i .e . , 

“ol + " B n = “o2 + B n+2 


(8-47) 


e = -e , 0 it follows from (8-45) and (8-46) that 
n n+z, 
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(8-48) 


F n,n+2 -F n+2,n 

which is as expected and says that there is no net 
force acting at the connection point when the split— beam 
is viewed as a whole. 

Similarly, from the translation equations of bodies n-1 
and n+1, the interface force at the connection point of 
these bodies can be obtained to be 


~ i 

P 1 ,4.1 = - me r, V 1 + 9 F - 

n-l,n+l n-1 n-1 ^ 

m 


'n+1 , n+2 n+1 , n+3 


-F - F 1 

c n-l , n c n-l , n-3 I 

-^n-l Sol + fctVfVj 


(8-49) 


F n+1 , n-1 me n+l 3 n+1 + 2 


c n+l , n+2 c n+l ,n+3 


F . +F 

L c n-l , n-3 c n-l,n 

n+3] + ^ [ ?e n-l e n+ll 


-me n +i “ 


(8-50) 


and again noting that £ n _i “ "^n+l and 


10 - + 8 t = 0) o + 8 n 

ol p n-l o2 K n+1 


(8-5i; 
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it is seen that 


F n-l,n+l F n+l,n-l 
which is as expected. 


(8-52) 


The interlace torque at the connection point between 
bodies n and n+2 is determined from the rotational 
equations of the bodies n and n+2: 


T +_y F +M +M +d . 9 F n+9 +M n+ o = H (8-53) 

e 1 T n e c . 0 c n,n+2 n,n+2 n,n+2 n 

n n n 7 n j n -*■ 


T + v F +M +M' +d F - + 

e n+2 2 11+2 e n+2 C n+2.n44 C n+2 . n+1 n+2 > n n+2 ’ n 


M n+2 , n H n+2 


(8-54) 


Since the net torque at the connection point should be 


zero.it follows that 
M n , n+2 - ~ M n+2 , n 


(8-55) 


Using (8-55) in (8-53) and (8-54) results in the interface 
torque 


M 


n, n+2 


1 " 

2 

- 


M +M 

c n+2 , n+4 c n+2,n+l 


-M 


1 1 

' ~ U o -U 

'n,n-l C n,n-2J 


T 

e n ,n+2 


(8-56) 


where T denotes the total external torque on bodies 

e n , n+2 

n and n+2. Use is made of the facts that 


F n,n+2 “ - F n + 2.n (equation (8 -4 8) ) 

*^n , n+2 ^n+2 , n 
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that the inertial rates of the bodies n and n+2 are the 


same 


n 


+ w 


ol 


«• 

= 6 


n+2 


+ « 


o2 


and that the inertia about the CM of bodies n and n+2 
is the same; i.e. , 



J 


n+2 


Similarly, the interface torque at the connection point 
of bodies n-1 and. n+1 may be determined to be 


M 


n-1 , n+1 



M +M -M -M 

L c n+l , n+3 c n+l , n+2 c n-l,n n-l,n-3 


1 

2 


T e 

e n-l ,n+l 


(8-57) 


where T e n "i n +l is the total external torque on bodies 
n-1 and n+1. Again, the interface torques are such that 


M n~l,n+1 M n+l,n-l 


(8-58) 


8.2.7 Translational Equations For Bodies of the Split-Beam 
The translational equation for a body i which is not 
involved in connections may be written, as before, as 


F +F +F +F 

e i c i , i-2 c i, i+1 c i,i+2 


oo 

m i p i 


l £ I 1} 1 3 


(8-59) 


8-30 



(8-60) 


F +F +F +F 

e. c - ■ t-, c. , c. . , 
x 1 , 1-2 i,l-l i,l+2 


O 0 

= m. Pi 




where the indicator functions are defined as 


1^ = {3,5,... ,n-3> 

I 2 = {4,6, . . . , n-2} 

Ig = {n+3,n+5, . . . ,2n-3} 

1^ = {n+4 , n+6 , . . . , 2n-2} 

From Figure 8-2, it follows that 


p i = 


„p + n r . +y . 

I o H 1 1 


(8-61) 


where the preceding subscript i can be either 1 or 2 
which, respectively, refer to half-beams #1 or 2; and 
'r. is the vector from the CM of half-beam #1 to the 

% l 

CM of body i in that half-beam. Differentiating (8-61) 
twice and retaining only the first-order terms in 
deformations, results in 



oo 

& p o 


ol 


X 


£ r i 


V. 


(8-62) 


For half-beam #£ 


F 

oo _% e 
£ P o M 

L l 


(8-63) 


where 0 F and M are, respectively, the total external 
force applied on the half-beam and the total mass 
of the half-beam #£ . Substituting equations (8-62) and 
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(8-63') in (8-59) and (8-60) results in the translational 
equations of motion for bodies not involved in connections. 
For example, the translational equation for bodies which 
are in the top row of the split-beam and are not involved 
in connection of the half-beams is given by 


m i^i + ( 2K+K i )^i ~ K ^i-2 “ K l u i+1 Ku i+2 ” K l d i } i+l 6 i 
+ Kd i-2 , i^i-2 + K l d i+l,i P i+l + ^ d i+2, i 6 i+2 


m i - • 

= F - F + m. .r- u „ 
e. e i l i o£ 

1 A JS, Z 


(8-64) 


£ = 1 for i £ 1^ 
H = 2 for i£ I 3 


Comparing (8-64) with (8—27) it can be seen that in 
general, the left-hand sides of the translational 
equations of the bodies (not involved in connections) of 
the split-beam .are identical to that of the corresponding 
bodies of the full-beam. Further, it may be noted that 
the right-hand side (rhs) of (8-64) is an appropriately 
modified form of the rhs of (8-27). 


The translational equations of motion for the four 
bodies involved in the connection of the half-beams are 
as follows: 


m 


mil -t — F -F — F .. - = F 

ln_T n n —l n+1 g 


n-1 


n — 1 c -i „ c -i „ q u 1 , n+1 e -> 1 o 

n-1 , n n— l,n— o n— l 


F + 


m n r , w , 
n-1 1 n-1 ol 
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m 

_p = F - F 

'n “c n c „ n , n +2 e M„ 1 e 
n,n-l n,n-<i n i-. 


-¥ 


+ ra n I r n u 


ol 


( 8 - 65 ) 


m 


m ^n+l ‘ F c n+ltn+ 2 c n +x t n +3 n+1 > n 1 e 


n+1 


n+1 M T, 


2 F e 


+ m n+l 2 r n+l (i) o2 
ra 


m ^ +2 ^n + 2,„ t r r °n + 2,n + 4' F 


= F 


n+2 


n+2 , n e n +2 M t, 


F 

2 e 


+ m n+2 2 r n+2 w o2 


Substituting into ( 8 - 65 ) for the constraint forces from 
section 8.2.2 and for the interface forces from section 
8.2.6 results in the translational equations of motion 
for the bodies involved in connections. These equations 
are given in Appendix F. 3 . 


8 . 2.8 Rotational Equations For Bodies of the Split-Beam 

• The rotational equation of motion for the bodies not 
involved in connections of the half— beams may be written, 
as in Section 8 . 2 . 4 , as 


T + • xF + M 

e i * 1 e i c i , i-2 


+ M 


c i, i+1 


o 

+ M = H. 

c i , i+ 2 


( 8 - 66 ) 


Z = 1 for i 6 I-l 
l - 2 for i ^ ^3 
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(8r-67) 


e i 


Vi xF ei 


. +M 


C i, i-2 


+M 


c i, i-1 


+M 


c i , i+2 


o 

= H. 


£ = 1 for i £ I 2 
£ = 2 for i €. 


As in Section 8.2.4, the rate of change of the angular 
momentum of body i may be obtained to be 


o 

H. 



1 ^ I I, 

J i ( “o2 + *i ) 

1 ^ I 3, 


( 8 - 68 ) 


Substituting for constraint torques from Section 8.2.2, 
equations (8-66) and (8-67) with the aid of (8-68) yield 
the rotational equations of motion for bodies not involved 
in connection. For example, for bodies in the top row 
of the beam, J ± 8 i + (2K 0 +K 0 _^) 3 1 ” K eBi-2" K e 1 ^i+l _K 0‘ e i+2 


+(2k. +K. ) 6-+K. 3. 0 +K. e.,,+K,$., 0 
ip i|; 1-2 l+l t|j 1+2 


+K, y.+K.y. _-K, y.,T-K,u.^ n 
(j)^ l (J) 1-2 l+l 4> K l+2 


T + n y.F -J.w 
e . 1 ’ 1 e . 10 

1 1 


(8-69) 


£ = 1 for i (* 

£ = 2 for i € I 2 
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Comparing (8-69) with (8-34), it is again seen that the 
left-hand sides of the rotational equations of the bodies 
(not involved in connections) of the split-beam are 
identical to that of the ■ corresponding bodies of the 
full-beam. Further, the right-hand side of (8-69) is 
an appropriately modified form of the rhs of (8-34). 


The rotationa.l equations for bodies involved in the 
connection of half-beams may be written as 


T +.y .F +M +M +d n , , F , , 

e 1 n-1 e t c t c - 0 n-1 , n+1 n-1 , n+1 

n-1 n-1 n-l,n n-l,n-3 ’ 


+M •, , — J -j(w t +3 -i ) 

n-1, n+1 n-l v ol n-1' 


(8-70) 


T + 0 Y +M +M +d 

e n+2 2 n 2 e n +2 c n+2,n+l c n+2,n+4 n 2,n 


+M „ = J , q ( CjO 0 +g , 0 ) 

n+2,n n+2 v o2 M n+2' 


F 

n+2, n 
(8-71) 


The rotational equations for bodies n+1 and n would be 
identical, respectively, to (8-70) and (8-71) because 
(i) The expressions for the interface torques have been 
derived using the rotational equations of bodies n-1 
and n+1 and bodies n and n+2, (ii) The inertial rate 
of the two halves of the beam is the same and hence for 
the bodies in connection. 


“ol 

** 

+ *n-l 

= “o2 

* »n + l 

(8-72) 

• 

*1 

• 

M 

(8-73) 

“ol 

+ *n 

<M 

0 

3 

1! 

+ 6 n+2 
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If w , f (o o, then (8-72) and (8-73) may be utilized 
as the rotational equations of bodies n + 1 and n, 
respectively. Substituting for the constraint and 
interface forces and torques from sections 8.2.2 and 
8.2.6 into (8-70) and (8-71) results in the rotational 
equations of motion for bodies n-1 and n+2. These equations 
are given in Appendix F.4. 

8. 2. .9 Matrix Equation of Motion for the Split -Beam 

If the inertial rates of the rigid body frames of the 
two half-beams are assumed to be equal, then these inertial 
rates must equal the inertial rate of the rigid body 
frame of the original beam. Hence, it is seen again that 
the rhs of the equations of motion developed in the pre- 
vious two sections are functions of only external forces 
and torques. 


For the beam under consideration, and the 

inertial rates of the half r beams are equal; viz,, - 

u . In such a case (8-72) and '(8-73) cannot be used 
as the rotational equations of bodies n and n+1 . 
Moreover, these equations imply that 


6 

/ 

3 


n+1 

n 


’n-1” 


n+2 


(8-74) 

(8-75) 


Thus, £ n and 3 n+ ^ ma y be eliminated from the equations 
of motion. The resultant translational and rotational 
equations of motion of the bodies composing the split- 
beam may be combined into the matrix equation 
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(8-76) 


M 2 Q 2 + K 2 Q 2 " L 2 

where Mg and K g are the mass and stiffness matrices, 
respectively, of the split-beam and are (12n-6) x(12n-6). 
L 2 is a (12n-6) column vector of external forces and 
torques; and Qg is u (12n-6) column vector 

that characterizes the flexible beam deformations 
(excluding the angular deformation of bodies n and n+1). 
and is defined by 



Let Kg be the 12n x 12n stiffness matrix of the split-beam 
that can be obtained from the equations of motion of the 
bodies of the split-beam. 


Since the left-hand sides of the equations of motion of 
the bodies which are not involved in the connection of 
the half-beams are identical to that of the corresponding 
bodies of the full beam, the rows of the Kg matrix 
corresponding to the bodies not involved in connections 
would be identical to the corresponding rows of the 
stiffness matrix of the full beam, K . So, let K' be 
the 12n x 12n matrix K q , but with rows 2n-3, 2n-2 , . . . , 2n+4 
which correspond to the equations of motion of the bodies 
in connection, set to zero. These rows are defined 
by the elements of a 12n x 12n matrix, K g , which has all 
zero entries for rows 1,2, . . . ,2n-4,2n+5,2n+6, . . . ,4n which 
correspond to the equations of motion of bodies not 
involved in connections arid for which the remaining rows 
are defined by the equations of motion of the bodies in 
connections. However, the entries in rows 2n and 2n+2 
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are all equal to zero due to the rotational equations 
((8-72) and (8-73)) of the bodies 2n and 2n+2. Once the 12nxl2n 
matrix "kj, = ~Kq +K q Is formecf , ' the" rows and columns 
corresponding to g~ and g n+1 have to be eliminated 

keeping in mind that g n = 3 n+2 and B n+1 = ^ n -l ^ for 
the case w ^ = w q2 )- This is accomplished by adding 
column 2n to column 2n+4 and column 2n+2 to column 
2n-2 and then deleting .the columns and rows 2n and 2n+2, 
thus resulting in a (12n-6) x (T2n-6) matrix which 
defines K 2 - 


The (12n-6) x (12n-6) mass matrix is obtained from the 
equations of motion developed in sections 8.2.7 
8.2.8 and noting the fact that 8 n = 3 n+2 and P n+1 = 
g which helps eliminate g n and 8 n+ -^ and results 
in the following: 
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8 . 2.10 


Mass Matrix M 2 = [ m ij] 


( 12 n- 6 )x( 12 n- 6 ) matrix 


m. . 
ij 

0 


3x3 matrix 
3x3 null matrix 




m. . 
ij 


= < 


m 

i=j=l, 3, . . . . ,4n-3 • 

J 

i=j=2,4, . . . . , 2n-4, J 

J+J* 

i=j= 2 n- 2 , 2 n +2 

m 

i=j= 2 n 

rai n-l 

(i, j )=(2n-3, 2n-2) 

mi n+l 

(i, 3 )=( 2 n, 2 n- 2 ) 

mS n 

(i,j)=( 2 n~l, 2 n+ 2 ) • 

me n +2 

(i, j )=( 2 n+l, 2 n+ 2 ) 

0 

otherwise 


An-2 


Determination of tlie Vibration Frequencies and Mode Shapes 
The equations of motion for the beam being considered 
can be written, as pointed out earlier, as: 

MQ + KQ = L ' (8-77) 

To determine the free-free modes of the system represented 
by (8-77), the external forces and torques on the rhs 
of (8-77) is set to zero. Thus, (8-77) may be written 
as : 

Q + RQ = 0 

where 

R = M X K 


(8-78) 

(8-79) 
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The eigenvalues, which will be real and positive if M 
and K are mass and stiffness matrices of a stable 
structural system, of the matrix R represent the square 
of the (natural) frequencies (in radians) of the modes 
of vibration which are possible in the system. The 
eigenvectors of the matrix R corresponding to the eigen 
values computed, characterize the mode shape of the 
vibration. Thus, if A* is an eigenvalue of R and X i is the 
eigenvector corresponding to A _j , then "\ZX~/2Tr is a natural 
frequency of the system and the components of the 
vector Xi describe the mode shape of the vibration. 


8.2.11 Stiffness Matrices et al 

For simulation purposes, a 20'x2 , x4' homogeneous beam 
made of steel is considered. The beam is modelled by 
the twenty 2 ft. cubes. For the material being 
considered, the modulus of elasticity, the modulus 
of rigidity and the density are, respectively, given 

by 

E = 29x10 psi 
G = 11 . 5x10® psi 
p = 0.283 lbs/in 3 

With the coordinate system as defined in Figure 8-6, 
the following may be defined: 


the axial spring constant 

the shear spring constant 

** 

the bending spring , 

constant 


K. 


y 



K x = K z = 

K = ^2. 
.^0 X 12a 


_5 be 
6 a 

E 


G 



bc^ 

12a 


E 
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the torsional spring constant 


K 0y a G 


where J is the polar moment of inertia and is given by 

P _ 

r , ,4 

J D = b 3 c 1/3 -0.21 £(1~ 
p c 12c^ 


With the above definitions, the stiffness matrices for 
body i in relation to bodies i-2, i+1, .and i+2 (see 
Figure 8-7) are as follows. 

(i) the translational stiffness matrix for bodies i and i+2. 


5 be r 

6 a~ G 


5 be 


(ii) the translational stiffness matrix for bodies i and i+1 


5 be 


| — G 0 
6 a 


be E 
a 
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(iii) the rotational stiffness 


matrix for bodies i and i+2 : 


(iv) 


K 


0 


X 


12a 

0 

0 


E 


0 


-z G 


0 


0 

0 

be 3 
12 a 


E J 


the rotational stiffness matrix for bodies i and i+1 


K 


6n 


r 3 

a c 

12b 
0 

0 


E 


0 

ac 3 

12b 

0 


E 


0 

0 

ip 

a 


G J 


where J 


PP 


a 3 c 


1/3 - 0.21 f (1 ^4 

c 12b* 


m 


The factor of one-half arises in the above . definitions 
due to the fact that the spring rates of the adjacent 
bodies are in series. 


The mass of each discrete body composing the beam can be 
computed from 


m = p v/ g 

where v is the volume of each body and g is the gravi- 
tational acceleration. 
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The moment of inertia of a body about its CM is given 
by: 

j = / (<$ T SI-<SS T )dm. 


where 6 is the vector from the CM of the body to an 
elemental mass dm and is given by 

S T = (ixjy.kz) 

and the elemental mass is given by 
dm = p dx dy dz/g 

Eor a body of a x b x c 


J = 


b/2 a/2 c 42 

(5 T 6I-6S T ) 

-Id/ 2 -a/2 -c/2 g 

Carrying out the integration results in the inertia dyadic 


III (<5 T SI-6S T )£ dx dy dz 


J 


m 

12 


a 2 +b 2 0 0 

0 b 2 +c 2 0 

0 0 c 2 +a 2 


The various vectors that are required for the simulation 
of the equations of motion of the beam are as follows: 


(i) the vectors from the CM of the beam to the CM of the 
body i: 


r ± = (0, - 




+ 1 


-a, (-1) 


i+1 b 


2 ' 
i=l,2, 


,n 


r 2n-i+l = " r i 
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(ii) 


( iii) 


(iv) 


where 



denotes the integer part of 


i+1 
2 4 


the vectors from the CM of the half-beam #1 to the CM of 
the body i in the half-beam #1: 


r . = (0 , - 

1 l v ’ 


1 - 
L 2 


, - .i+1 b. 
a, ( -1) f ) ) 


2 

i=l,2. 


n 

2 


l r n-i+l l r i 

l r j = 2 r n+j J=1,2, • • * > n 

the vectors from the CM of body i to the interface 
point that lies on the line joining the CM of body i 


the 

CM of 

the 

adj 

oining 

body : 



d. 

a, 

, i+2 

(0, 

a 

2’ 

0) 

i=l, 

2,.. 

. . , 2n-2 

d i; 

, i-2 

(0, 

a 

” 2’ 

0) 

i=3, 

4 , . . 

. . ,2n 

d i, 

, i+1 

(0, 

o, 

_ b) 
2 } 

i=l> 

3, . . 

. . ,2n-l 

d. 

i. 

, i-1 

(0, 

o, 

2 ' 

i=2, 

4,.. 

. . , 2n 


the vectors from the connection point to the CM of body i 
£ n-l = e n = ~ £ n+l = ~ e n+2 = (0 ' " 2’ 0) 

Finally, it may be pointed out that, for the beam 
example being considered, 

a = b = c = 2 ft. 
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8 . 2.12 


Simulation Results 

The eigenvalues and the corresponding eigenvectors are 
obtained, as discussed in Section 8.2.10, for 
both the full beam and the split-beam. All the eigen- 
values are, as expected, non-negative in both cases. 

For the full beam, there are six zero eigenvalues which 
are associated with the six rigid body modes of the 
beam. The remaining 114 eigenvalues yield the beam's 
natural frequencies ranging from 58 . 2Hz to 2066.8Hz. 

The eigenvectors corresponding to these 114 eigenvalues, 
the imaginary parts being zero, are real and hence 
describe the mode shapes of the vibrations. Due to the 
truncation and roundoff errors ' in computations, four 
of the eigenvalues corresponding to the rigid body modes 
(instead of being zero) are complex conjugates with 
the real and imaginary parts being less than 10 and 
the other two eigenvalues are real and are less than 
10~ 5 . Thus, the computed eigenvectors corresponding to 
these rigid body modes do not portray the shapes of the 
rigid body modes. 

For the split-beam, there are twelve zero eigenvalues, 
six of which are associated with the rigid body modes 
of the beam and the rest of which result due to the 
way the half-beams are put together. The remaining 102 
positive eigenvalues yield natural frequencies ranging 
from 66.8Hz to 2014.5Hz. Thus, there are twelve 
frequencies of vibration of the beam which are missing 
for the split-beam case; otherwise, the frequencies 
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of vibration of the split-beam compare very favorably 
with the natural frequencies of the full beam, the 
maximum error being 14.8%. A comparison chart of the 
frequencies of vibration obtained for the full beam and 
the split-beam is shown in Table 8-2. 

The mode shapes of the non-zero natural frequencies of 
the split-beam compare, in general, quite favorably 
with the corresponding mode shapes obtained for the 
original beam. Since the angular deformations of two 
of the four bodies involved in the connection of the 
half-beams have been eliminated from the equations of 
motion of the split— beam, the mode shapes corresponding 
to these deformation coordinates are absent in the 
eigenvectors of the split-beam. In some cases, the 
relative amplitudes of vibration of the bodies involved 
in connections, given by the components (of the eigen- 
vectors) corresponding to the deformation coordinates 
of the bodies in question, are close to zero. However, 
it may be pointed out- that the corresponding components 
of the eigenvectors- of the original beam are small, 
generally less than 20% of the maximum component of the 
eigenvector being considered. 

In summary, the lower-order modes (that is, modes 
corresponding to the frequency range from the lowest 
to a frequency at least ten times the lowest frequency) 
of the split-beam compare very well with those of the 
original beam. The natural frequencies of the beam 
that are missing from those obtained for the split-beam 
are all in the high frequency range. Further, there 
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TABLE 8-2 


COMPARISON OF THE FREQUENCIES OF VIBRATION 
FOR THE ORIGINAL BEAM AND THE SPLIT-BEAM 

ORIGINAL BEAM SPLIT -BEAM 


58.2 

66.8 

102 _4 

115.1 

125.5 

140.6 

151.4 

153.5 

225.6 

233.1 

249.9 

252.7 

275.6 

305.3 

291.9 

324.1 

363.7 

384.6 

372.8 

409.3 

418 

429.4 

495.4 

504.3 

495.6 

509.6 

568.7 

622.0 

576.7 

576.7 

617.7 

652.7 

618.7 

667.5 

623.2 

657.9 

719.3 

746.5 

735.4 

756.8 

741.7 

767.2 

780.3 

783.2 

847.2 

* 

859.8 

917.4 

862.3 

916.2 
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TABLE 8-2 
CONT ' D 

ORIGINAL BEAM 

873.4 

933.1 

934.5 

948.3 
965 

989,9 
991 

1024.8 
1037 

1044.1 

1054.1 

1096.9 

1128.6 

1131.3 

1145.4 

1184.7 

1192.3 

1192.5 

1199.4 

1200.4 

1201.5 

1205.9 • 

1209 

1229.8 

1263.7 

1278.4 
1289.3 

1290.8 


SPLIT-BEAM 

929.1 

933.1 
974.9 
951.8" 
984 
1045 
991 

1071.6 
* 

* 

1074.1 

1096.9 

1128.6 

1139.5 

1157.2 
1128 

1192.5 

1192.6 

1199.4 
1202.8 

1201.9 

1247.2 

1181.4 

1231.5 

1219.7 
1278.4 

1291.3 

1294.7 
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TABLE 8-2 
CONT'D 


ORIGINAL 

1307.9 

1309.9 

1316.3 
1319.6 
1319.6 
1319.6 

1324.9 

1340.5 

1341.1 

1358.8 
1359 

1373.1 

1376.5 

1387.9 

1401.3 

1402.3 

1409.9 

1411.9 

1417.4 

1431.3 

1433.4 

1441.1 

1461.6 

1466.3 
1503.3 ' 
1509.8 

1513.2 

1517.4 


BEAM SPLIT-BEAM 

1317.8 

1311.3 

1362.4 

1319.6 
* 

* 

1344.4 
* 

1346.3 
* 

1307.3 

1374.2 
* 

1387.9 

1399.6 

1419.8 

1410.4 

1427.3 

1378.8 

1394.9 

1408.6 

1485.7 

1413.8 

1492.7 
* 

1509.8 
1515.7 
1565.1 


8-50 



TABLE 8- 
CONT ' D 


ORIGINAL BEAM 

1529.2 
1530.5 
1556 

1579.3 

1579.8 

1600.8 

1601.3 

1601.9 

1611.3 
1627 

1635.7 

1635.9 

1656.4 

1656.5 

1662.8 

1663.4 

1686.9 

1690.1 

1713.5 
1721 
1723 

1751.2 

1759.7 

1761.6 

1774.8 

1786.3 

1786.9 
1800.2 


SPLIT-BEAM 

* 

1546.6 
1554 

* 

1562.9 

1564.7 

1601.4 
1601.2 

* 

1656.3 

1624.8 
1636 

1634.1 

1623.5 
* 

1669.2 
* 

1691.5 

1719.6 
* 

1657.8 

1724.5 

1759.7 

1761.6 

1774.8 

1756 .4 

1752.9 

1829.2 
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TABLE 8-2 
CONT ’ D 


ORIGINAL BEAM 


SPLIT-BEAM 


1843.2 

1866 

1944.8 
2017.5 

2066.8 


1753.6 

1850 

2002.2 

2014.5 

* 


♦Denotes frequencies for which there are no corresponding 
mode shapes. 
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are four split-beam's natural frequencies for which 
there are no corresponding mode shapes in the original 
beam. However, these are again in the high frequency 
range. The fact that the higher-order modes of the 
split-beam do not compare as well as the lower-order 
modes do with the corresponding modes of .the original 
beam, poses no problem since the control systems 
associated with spacecraft are, in general, low frequency 
systems . 

The computer outputs showing the natural frequencies 
and the mode shapes for the original beam and for the 
split-beam are given in Appendices G and H, respectively. 

8.3 CONCLUSIONS 

The computer simulation results indicate that the 
frequencies of vibration and the corresponding mode 
shapes obtained from the free-free modes of the half- 
beams and using the method developed for combining the 
half-beams compare very favorably with those obtained 
for the original beam using the standard discrete 
coordinate formulation. The distortions that are 
encountered in the modes of the split-beam are mostly 
in the higher— order modes. These distortions are 
mainly due to the fact that the dynamic model of the 
split-beam includes, in effect, only 18 bodies whereas 
the original beam is modeled with 20 bodies. Two of 
the bodies are "lost" because, for the split— beam, there 
are only 18 rotational equations and two of the 20 
translational equations give rise to six more rigid 
body modes that are not in the original beam. Thus, 
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ten percent of the components of the system are, in 
effect, "lost" when the new modeling technique is 
utilized for a beam modeled with a small number of 
bodies. In practice, however, the number of bodies that 
constitute a spacecraft is large and the number of 
components involved in connections would be very small 
thus resulting in a very small percentage of the bodies 
being "lost" when the modeling technique being 
considered is used. Thus;, in practice, it is expected 
that the difference between the flexible body characteris- 
tics obtained via the new method and the [actual flexible 
body characteristics would be insignificant. 


Hence, it can be concluded from the above remarks 
and the simulation results obtained that the method 
(reference 8-1) of combining substructures of a system 
achieves a good degree of fidelity in the determination 
of .flexible body characteristics for -the overall system 
from the characteristics of the substructures constituting 
the system. 
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SECTION 9 


9.0 CONCLUSIONS AND RECOMMENDATIONS 

Equations have been generated for three different 
approaches to controlling a large flexible spacecraft. 

Two of the three (LQR and MNA) were simulated con- 
trolling a single axis of the Space Construction Base 
(SCB), Configuration 1. The third approach, multilevel 
control, was developed for a three axis SCB simulation; 
but the computer program was too complex for the actual 
simulation to be completed before the end of this re- 
porting period. 

One control approach utilized some conventional optimi- 
zation techniques based upon design in the time domain. 
This approach (LQR) can be easily expanded in dimen- 
sions (more axes, actuator locations, flexible modes, 
etc) , readily extended to include noise sources at 
sensors and actuators, and the filtering of state mea- 
surements . 

The second control method (MNA), is esentially based 
upon design in the frequency domain. It utilizes system 
matrix manipulation to obtain diagonal dominance, there- 
by reducing the problem to a set of independent single 
loop designs . 

The third control method utilizes a multilevel hierarchal 
technique which is amenable to implementation by a digital 
processor. Optimization of the control system is obtained 
through a performance index using weighting factors on 
components of the state variable and control vectors. 
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The other portion of- the study was related to mathe- 
matical modeling of flexible spacecraft dynamics. The 
previously generated SCB vehicle equations were derived 
for large angular displacements of the central body and 
the appendages (solar wing elements). The model was 
simplified considerably by restricting motion to small 
angular displacements. This linearized model was then 
utilized in single-axis control simulations and the 
multilevel, three-axis, simulation equations. The 
three axis linearized equations were also examined with 
respect to MNA control while including the cross-coupling 
between axes (the full order mode, FOM) . The reduced 
order model (ROM) was generated to study control of the 
SCB without cross -coup ling, thereby reducing the problem 
from one larger three-axis analysis to three much sim- 
pler single-axis analyses. 

A structural analysis study (Section 8) illustrated a 
technique whereby the structural characteristics of two 
or more submodules may be combined when the submodules 
are firmly connected, as through docking or berthing ports. 
Prediction of the resulting dynamics- may be based upon 
the theory presented. This was verified by means of an 
example of' connecting two shorter beams to obtain one 
which is longer. 

9.1 CONCLUSIONS 

The following conclusions are based upon the formulation 
of control equations and the simulations of Configuration l 
of the Spcae Construction Base: 
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a. Excellent simulation results were obtained when 
feeding back the full state variable vector (for 
LQR, MNA and Retallack pole placement). 

b. The results for partial feedback of the state 
variable vector were satisfactory, but power 
expended at the central body actuator system was 
excessive when appendage deflection states were 
omitted (LQR) . 

c. The MNA approach also produced good performance for 
the difficult case of vehicle flexible modes being 
within the closed loop bandwidth of the attitude 
control system. 

d. The Retallack pole placement method was very success- 
ful in placing specific open loop poles in desired posi- 
tions when the control system loop is closed. 

e. The three-dimensional, five-body, linearized 1 vehicle 
model was transformed into a scalar state variable 
equation set which was readily programmed for a 
digital computer. 

f. The three-axis, five-body multilevel control method 
was similarly transformed into an equation set 
amenable to programming. 

g. Cross-axis coupling of the subject vehicle was vir- 
tually negligible; that is, differences between 
the reduced order model (ROM) and the full order 
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model (POM) were negligible. 

In addition, the structural analysis technique developed 
in Section 8 was fairly successful in predicting the 
frequencies of vibration and mode shapes of the "split- 
beam” structure where the two submodules are fastened 
end-to-end. This is particularly fortunate since the 
number of finite elements used was not very high; 20 
bodies were used in defining the original full beam 
and two sets of 10 bodies were used for each of the two 
connected submodules. Finite element modules for space- 
craft in later stages of design usually are defined to 
a much higher degree; the fidelity of the method des- 
cribed, therefore, is actually expected to be much better 
when applied to clusters of "real" space vehicle 
modules . 

9.2 RECOMMENDATIONS 

The following items are recommended for future studies 
regarding control of large flexible space vehicles; 

a. All three control methods should be expanded to 
stabilize space vehicles modeled with more discrete 
elements and. dimensions. 

b. Similarly, the control should be expanded to stabilize 
vehicles modeled in three axes, including cross- 
coupling. 
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c. More attention must be devoted to quick, accurate 
generation of spacecraft mathematical models. For 
vehicles such as Configuration 1 of the SCB, the 
hybrid coordinate technique is an excellent candi- 
date for this type of study because of the flexible 
solar wings. Decisions for truncating modes may 
then be more confidently made on the basis of modal 
dominance and frequencies . 

d. Noise on sensor measurements of states and at actu- 
ators may be included in future studies; these are 
customarily extensions to the LQR method. 

e. The MNA and pole placement method should be expanded 
to placement of closed loop poles with regard to 
the magnitude of frequency (in addition to damping 
factor) . 

f. Control simulations of the LQR and MNA with pole 
placement methods should be repeated for the X and 
Y axes of the same vehicle configuration. 

g. Because of uncertainties of vehicle model parameters, 
sensitivity analyses must be included in control 
system evaluation to determine variations in per- 
formance . 

h. It may be desirable to try other pole placement tech- 
niques (the standard Simon and Mitter technique and 
the Kalman pole placement algorithm), and evaluate 
the subsequent designs and performance. 
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i. A test should be developed for robustness of the 
control design as a means of simplifying design for 
future configurations. 

j . The three-axis multilevel control simulation should 
be completed (the SCB Configuration 1 model equations 
have previously been programmed). 

k. When access to measurement of essential state vari- 
ables is impractical, filtering techniques must be 
developed for generation of observers required for 
the inaccessible states . 

l. The LQR method should be expanded for stabilization 
of two or more locations on a large flexible space 
vehicle; for example, when multiple Instrument 
Pointing Systems (IPS) are used. There’ may be a 
need for stabilization of all six degrees of freedom 
at the mounting location, particularly when the 
center of mass of the IPS inner gimbal and payload 
is not located at each gimbal pivot axis . 


• It is recommended also, that the structural analysis 
technique that was developed and evaluated on beams with 
a limited number of elements, be verified for a high 
degree of fidelity on more complex structural models 
with many more finite elements. Further, the problem 
of substructure model truncation prior to determining 
the overall system’s characteristics should be examined. 
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APPENDIX A 


NUMERICAL VALUES OF 
MATRICES USED IN 
ANALYSIS AND DESIGN 


Space Construction Base, Configuration 



Values of matrices P, N, KD, KS, A, and B used for 
analysis and control system design purposes are shown 
on the following pages. A list of constants used in 
the model is also shown in this Appendix. 

Most of the parameter values used in the simulation of 
five-body three-dimension model of Configuration I are 
based upon the request issued by Marshall Space Flight 
Center . 
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0000000E 00 0. 0008000E 00 

0000000E 80 0. 0000000E 00 
0000000E 08 0. 0000000E 00 


000000 0E 
3148030E 

0800000E 

8000000E 


00 0. 0000000E 00 
04 0. 0000000E 00 


00 0. 0000000E 00 
00-0. 3203908E 01 


0. 0000000E 88 
0. 2677729E 03 
0. 0000000E 00 
0. 0000000E 00 
0..0000000E 00 
0. 0000000E 00 
0. 0000800E 00 
-0. 2677729E 03 
8. 0008000E 00 
0. 0000000E 00 
0. 0000000E 00 
0. 0000000E 00 
0. 0000000E 00 
0. 0000000E 00 
0. 0000000E 00 


KD ARRAY 



V P o 
^po 

J} 

p ? 

o ^ 

A- f 


BEND IX RESERRCH LABORATORIES 


*** KS STIFFNESS MATRIX +** 

< TWO LINES PER R ROW > 

0 000PI080E 00 0. 0000000E 00 0. 0000000E 00 0. 1966330E 05 0. 00O0000E 00 0. O000000E 00 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 000O0Q0E 00 0. 0Q0O00OE 00 0. 0000000E 00 

0 000000 0E 00 0. 0000000E 00 0. 0000000E 80 0. 0000000E 00 0. 5998410E 07 0 8OO0000E 00 

Q. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 000000OE 00 O. 00O000OE 00 

0 000B000E 00 0. 0000000E 00 0. 00000O0 E 00 0 0000000E 00 0. 0000000E 00 0. 65572O0E 03 

0* 6557200E 03 0. 0008080E 00 0. 0000080E 80 0. 0000000E 00 0. O00O000E 00 0. 0000O00E 00 

0. 0000000E 00 0. 0080000E 00 0. 0000000E 00-0. 1966330E 05 0. 090O000E 00 0. 0000000E 00 

0! 0000000E 00 0. 1841408E 04 0. 008000OE 00 0. 0000OO0E 00 0. 000000OE 00 0. 000O000E 00 

0. 0000000E 00 0. 0008000E 00 0. 0080000E 00 0. 0000000E 00-0. 5998410E 07 0. 0O000O0E 00 

0. 0000000E 08 0. 0000000E 00 0. 7024030E 06 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00-0. 6557200E 03 

0. 0000000E 00 0. 0000 000 E 00 0. 0000000E 00 0. 89549O0E, 02 0. 0000Q08E 00 0. 0000000E 00 

* 

0 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0800000E 00 0. 0000000E 00- 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000O00E 00 0. 1841403E 04 0. 0000080E 00 

0 0000000E 00 0. 0000000E 00 0. 0000800E 00 0. 0000000E 00 0. 0000080E 00 0. 0000000E 00 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 00Q0000E 00 0. 7024030E 06 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 

-0. 6557208E 03 0. 0000000E 00 0. 0000000E 00 0. 0000008E 00 0. 0000000E 00 0. 0030003E 00 

0 00 0000 0E 00 0. 0003000E 00 0. 0000000E 00 a 0000000E 00 a 0000000E 00 a 0000000E 00 

a 0000000E 00-0. 1841400E 04 0. 0080000E 00 0. 0000OO0E 00 0- 00000Q0E 00 0. 00O0000E 00 

0 0000030 E 00 a 0800000E 00 a 0008000 E 00 a 0000000E 00 0. 0000000E 00 a 0000000E 00 

a 0080000 E 00 a 0000000E 00-a 7024030E 06 a 0000000E 00 a 0000000E 00 a 0000000E 00 

a 0000000E 00 a 0000000E 00 a 0000000E 00 0. 0000000E 00 a 0000000E 00 a 00 

0. 0000000E 00 0. 0000080E 00 0. 0000000E 00-0. 8954900E 02 0. 0000O00E 00 0. 0O0O000E 00 

0 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 a 0000000E 00 a 0000080E 00 

a 0000000E 00 0. 0080000E 00 0. 0000030E 00 0. 00000Q0E 00-0. 1841400E 04 0. 0000000E 00 

a 0000080E 00 0. 00000O0E 00 a 0008000E 00 a 0000000E 00 a 0000000E 00 a 00 

0 . 0000000E 00 a 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00-0. 70248- CE 06 

a 000B000E 00 0. 0080000E 00 a 8030000E 00 a 0000000E 00 0. 0000000E 00 a 0000000E 00 

0 . 0000000E 00 0. 0000000E 00 a 0000000E 80 0. 0000000E 00 0. 00O0000E 00 0. 0000000E 00 


1966830E 

0000000E 

0000000E 

0000000E 

0000000E 

0000000E 

0000000E 

0000000E 

0000000E 

0000000E 

0000080E 

0000008E 

. 1966830E 
. 0000000E 

. 0000000E 
. 0000000E 

. 0000000E 
. 895490 8E 

. 8008000E 
. 0000000E 

. 0000800E 
. 0000000E 

. 0000000E 
. 0000000E 

. 0000000E 
. 0Q0000QE 

I. 0000000E 
I. 0000000E 

I. 0000000E 
I. 8954900E 


05 0. 0000000E 00 
00 

80 a 5998410E 07 
00 

00 0. 0000000E 00 
00 

00 0. 0000000E 00 
00 

00 0. 0000000E 00 
00 

00 a 0000080E 00 
00 

05 a 0000000E 00 
00 

00-0. 5998410E 07 
00 

00 0. 0000000E 00 
02 

00 0. 0000000E 00 
00 

00 0. 0000000E 00 
00 

00 0. 0000000E 00 
00 

00 8. 0000008E 00 
00 

00 0. 0000000E 00 


00 0. 0000000E 00 
02 


KS ARRAY 



BEND IX RESERRCH LABORATORIES 


*** All PART OF R -THE PLANT MATRIX <ROMS=l, 13 ; COLS=l, 13) +** 

< TWO LINES FER A ROW ) 

“ 0. 0008000E 08 0. 0000000E 00 0. 0000000E 00 0. 1288899E-02 0. 0000000E 00-0. 5771S43E-12 0. 1280399E-02 0. 000000OE 00 

-0. 5771543E-12-0. 6926712E-03 0. 000000OE 00 0. 0000Q00E 00-8. 6926707E-03 0. 0000800E 80 0. 0000000E 00 

0. 0000800E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 1000111E-02 0. 00080S0E 00 0, 0000000E 00 0. 1000111E-82 

0. 0000000E 00 0. 0000000E 00-0. 1348562E-02 0. 0000000E 00 0. 00O0000E 00-0. 1348562E-02 0. 0080000E 00 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00-0. 7501572E-04 8. 0008800E 00 0. 1336089E-04 0. 7501571E-04 0. 6000080E 00 

0. 1336089E-04 0. 4367963E-03 0. 0000000E 00 0. 000090OE 00-0. 4367961E-03 0. 0008O00E 00 0. 8O00000E 00 

0. 0000080E 00 0. 0000000E 00 0. 0008000E 00-0. 2962485E-02 0. 0000808E 00 0. 7431766E-05-0. 1354121E-02 0. 0000000E 00 
0. 7431766E-65 0. 1767937E-01 0. 00000OOE 00 0. 0000000E 00 0. 42634S7E-03 0. 0000800E 00 0. 0000000E 00 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000800E 00-0. 1118445E-01 0. 0080000E 00 0. 0000000E 00-0. 9886064E-83 

0. 0000O00E 00 0, 0000000E 00 0: 29S3305E 00 0. 0000O00E 00 0. 00Q0000E 08-0. 1739376E-03 0. 00000O0E 00 

0. 0000000E 08 0. 0000000E 00 0. 0000000E 00 0. 7501571E-04 0. 0000000E 00-0. 4751936E-02-0. 7501572E-04 0. 8O00000E 00 

-0. 1336089E-04-0. 4367963E-03 0. 0008000E 80 0. 4231419Er02 0. 4367961E-03 0. 00O000OE 00 0. 0000080E 00 

0. 0000800E 00 0. 0000000E 80 0. 0800008E 00-0. 1354121E-02 0. 0000O00E 00-0. 7421767E-05-0. 2962486E-02 0. 0000000E 00 
-0. 7431767E-05 .0. 4263471E-03 0. 0000000E 00 0. 0000000E 00 0. 1767937E-01 0. 0000000E 00 0. 0000000E 00 

0. 0008600E 00 0. 0000000E 00 0. 0000000E 00 0. 0000800E 00 ~ 0 . 9886064E-03 0. 8000080E 00 0. 0000900E 00-0. 1118445E-01 

0. 0000000E 08 0. 0000000E 80-0. 1739395E-03 0. 0080000E 00 0. 0000O0QE 00 0. 2963274E 00 0. 0O0O0O0E 00 

0. 0000000E 00 0. 0000008E 00 0. 0800000E 80 0. 7501571E-04 0. 0000008E 00-0. 1336089E-04-0. 7501572E-04 0. 8000080E 00 

-0. 4731936E-02-0. 4367963E-03 0. 00000OOE 00 0. 0000000E 00 0. 4367961E-03 0. 0800000E 00 0. 4231419E-02 

> 0. 0000000E 00 0. 0000000 E 00 0. 0008000E 00 0. 4130814E-02 0. 0000000E 00-0. 1011885E-04 0. 9961669E-04 0. 00O0000E 00 

^ -0. 1011085E-84-0. 5660854E-01 0. 0000000E 00 0. 08000O0E 00 0. 3623310E-03 0. 0000000E 00 0. 0000000E 00 

0. 0000080E 00 0. 0000000E 00 0. 0080000E 80 0. 0000000E 00 0. 2528601E-81 0. 0000000E 00 0. 0000B00E 00-0. 1479540E-04 

0. 0000000E 00 0. 0000000E 00-0. 9839008E 00 0. 0008Q00E 00 0. 0000000E 00 0. 1957956E-02 0. 0000000E 00 

0. 0008000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000800E 00 0. 0000008E 00 0. 4738573E-02 0. 0000000E 00 0. 000000QE 00 

0. 0003000E 00 0. 0000000E 00 8. 0000000E 00-0. 8462839E-02 0. 0000000E 00 0. 0000000E 00 0. 000QQ00E 00 

0. 0000000E 00 0. 0000000E 00 0. 0008000E 00 0. 9961659E-04 0. 080O000E 00 0. 1011084E-04 0. 413081 4E-02 0. 0000000E 00 
0. 1011084E-04 0. 3 623313 E-03 0. 000O008E 00 0. 000O0O0E 00-0. 566OS54E-01 0.- 0QO0000E 00 0. 000OO00E 00 • 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00-0. 1479524E-04 0. 0000000E 00 0. 0000000E 00 0. 2520575E-01 

0. 00O0000E 00 0. 0000000E 00 0. 1957956E-82 0. 00008Q0E 00 0. 000O000E 00-0. 9828948E 00 0. 000O080E 08 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 00O00O0E 80 0. 000000OE 00 0. 00000O0E 00 0. 0000000E 00 0. 0000000E 00 

0. 4738373E-02 0. 0000080E 00 0. 00000OOE 00 0. O0O0000E OO 0 . 0O0O000E 00 0. 0OO0000E 00-0. 8462S39E-02 


A = A 11 A 12 , All PART OP A MATRIX 


I 


O 



BENDIX RESEARCH LABORATORIES 


> 

I 

~4 


*** A12 PART OF A -THE PLANT MATRIX <ROWS=l> 15 jCOLS=16.30) *** 

< TWO LINES PER A ROW > 

0. 0000000E 00 0. 0000008E 00 0. 0000000E 00 0. 6956348E 00 0. 0000000E 00-0. 1054800E-09 0. 635634SE 00 0. 0000000E 00 

-0. 1054800E-09-0. 8228935E-02 0. 000O000E 00 0. 0008000E 00-O. 8228928E-02 0. 00000O9E 0O 0. Q000O00E 00 

0. 0000000E 00 0. 0000008E 00 0. 0000000E 00 O. 0000000E 00 0 2240359E 02 0. 00008O0E 00 0. 0000008E 00 0. 2240359E 02 

0! 0000O00E 00 0. 0000000E 00-0. 2008974E 00 6. 000OQ08E 00 0. 00O0000E 00-0. 30O8974E 00 0. O0O0000E 00 

0. 0000000E 00 0. 0800000E 00 0. 0000000E 00-0. 4073977E-01 0. 0600080E 00 0. 2441819E-02 0. 4073976E-01 0. 0000000E 00 

0 ! 2441819E-02 0. 5189139E-02 0. 000OS00E 00 0. 0000000E 00-0. 5189139E-02 0. O0000O0E 00 0. 8000000E 00 

0. 0000000E 00 0. 0000800E 00 8. 0000000E 00-0. 1608875E 01 0. 0000800E 00 0. 135S220E-02-0. 7354001E 00 0. 0000000E 00 

0. 1358220E-02 0. 2100309E 00 8. 000O00OE 00 0. 0000000E 00 0. 5064998E-02 0. 0O00008E OO 0. 0000000E 00 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0 . 0000000E 00-0. 2505439E 03 0. 0000 OQ 0 E 00 0. 0000000E 00-0. 2214587E 02 

0. 0000000E 00 0. 0000000E 00 0: 6611862E 02 8. 0000000E 00 0. 0000000E 00-0. 3880975E-O1 0. 0006O00E 00 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 4073977E-01 8. 0000000E 00-0. 8684571E 00-0. 4073976E-01 0. 0000000E 00 

-0. 2441819E-02-0. 5189139E-02 0. 00O0000E 00 0. 1182682E .00 0. 5189129E-02 0. 0000000E 00 0. 0000000E 00 

0 0000000E 00 0. 0000008E 00 0. 0000000E 00-0. 7354001E 00 0. 0000000E 00-0. 1358220E-02-0. 1608376E 01 0. 0000000E 00 
-0. 1358220E-02 .0. 5065005E-02 0. 0000000E 00 0. 0000080E 00 0. 2100309E 00 0. 0000800E 00 0. 0000000E 00 

0. 0000000E 00 0. 0000000E 00 0. 0800000E 00 0. 0000000E 00-0. 2214587E 02 0. 0000880E 00 0. 0000000E 00-0. 2505439E 03 

0. 0000000E 00 0. 0000000E 00-0. 3881019E-01 0. 0000008E 00 0. 0G00000E 00 0. 6S11795E 02 0. G008O00E 60 

0. 0000000E 00 0. 0000000E 00 0. 0000800E 00 0. 4073977E-01 0. 0000000E 03-0. 2441819E-02-0. 4073976E-01 0, 0000O0OE 00 

-0. 8884571E 00-0. 5189129E-02 0. 0000000E 00 0. 0000000E 00 0. 5189139E-02 0. 0000O00E 00 0. 1182682E 00 

0. 0000030E 00 0. 0000030E 03 0. 0000000E 00 0. 2243375E 01 0. 0000800E 00-0. 1847846E-02 0. 5410015E-01 0. 0000000E 00 

-0. 1847846E-02-0. 6725094E 00 0. 0000000E 00 0. 0000000E 00 0. 4304495E-02 0. 0000O88E 00 0. 0000000E 69 


0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 5646426E 03 0. 0000000E 00 0. 0000000E 00-0. 3314323E 00 

0, 0000000E 00 0. 0080000E 00-0. 2195325E 03 0. 0000000E 00 0. 0000080E 00 0. 4368683E 00 0. 0000000E 80 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0008000E 00 0. 8660153E 00 0. 0009000E 80 0. 0000000E 00 

0. 0000000E 00 0. 0000000E 00 0. 0000000E 00-0. 2365264E 00 0. 0080000E 00 0. 0000000E 80 0. 0000000E 00 

0. 0000080E 00 0. 0000000E 00 0. 0000080E 00 0. 5410010E-01 0. 0000800E 00 0. 1847S45E-02 0. 2243375E 01 0. 000000OE 00 
0. 1847845E-02 0. 4304495E-02 0. 0000080E 00 8. 0000000E 00-0. 6725094E 00 0. 0000000E 08 0. 0908000E 00 


0 . 

0 . 


0000000E 00 0. 0000000E 00 0. 0000000E 00 0. 0000000E 00-0. 3314285E 00 0.‘0080000E 00 0. 0000000E 00 0. 5646367E 03 
0000000E 08 8. 0000000E 00 8. 4368683E 00 0. 000O0O0E 00 0. 000QO00E 80-0. 2195212E 03 0. 00000O0E 00 


0. 0000080E 00 0. 0000000E 00 0. 0000000E 00 8. 0000000E 00 0. 08O0000E 00 0. 0000000E 00 0. 8000000E 00 0. 0800000E 00 
0* 86S0153E 00 0. 0090000E 00 0. 0000000E 00 0. 000O000E 00 0. 00Q0000E 00 0. 0000000E 00-0. 2365364E 00 


A12 PART OF A MATRIX 



BEND IX RESEARCH LABORATORIES 
*** B ARRAY CCONTROL MATRIX) FIRST 15 RONS *** 

0. 3208360E-04 0. 0000000E 00-0. 1608613E-12-0. 3284731E-05 0. 0000000E 00 0. 08000O0E 00-0. 3284727E-05 0. 000O000E 00 0. 0000080E 00 

0. 0000000E .00 0. 3401813E-05 8. 0000000E 00 0. 00O0O00E 00-0. 3331886E-06 0. 0000000E 00 0. 080O000E 00-0. 3331086E-06 0. 0000800E 00 

-0. 1608613E-12 0. 0000000E 00 0. 3723875E-05 0. 2071341E-05. 0. 0080000E 00 0. 000000OE 00-0. 2071341E-05 0. 0000000E 00 0. 0000000E 00 

-0. 3536833E-04 0. 0000000E 00 0. 2071341E-05 0. 4643211E-04 0. 0000000E 80 0. 00O0800E 00 8. 2021786E-05 0. 0000000E 00 0. 0000080E 00 

0. 0008000E 00-0. 3734922E-05 0. 0800000E 00 0. 0080000E 00 0. 3803347E-84 0. 0000000E 00 0. 0009900E 00-0. 4296493E-07 0. 0000800E 00 
0. 1608613E-12 0. 0000000E 00-0. 3723875E-05-0. 2071341E-05 0. 0000000E 00 0. 1320709E-02 0. 2071341E-05 0, 000080OE 00 0. 0000000E 00 
- 0 . 3536833E-04 0. 0000000E 80-0. 2071341E-05 0. 2021789E-05 8. 0000000E 00 0. 8000800E 89 0. 4643212E-04 0. 0O0O000E 00 0 . 0800080E 00 

0. 8000000E 00-0. 3734922E-05 0; 0800000E 00 0. 0000000E 08-0. 4296493E-07 0. 8Q00000E 00 0. 8000000E 80 0. 3883347E-04 0. 0000080E 08 

0. 1608613E-12 0. 0080000E 00-0. 2723875E-05-0. 2071341E-05 8. 0000008E 00 0. 0O08000E 89 0. 2071341E-05 8. 008000OE 80 0 . 1320789E-02 

0 . 4468847E-05 0. Q000000E 00-8. 2818841E-05-0. 1095916E-83 0. 080O080E 00 8. 8000880E 89 8. 1718220E-05 8. 8800009E 00 0 . 0000800E 00 

0 . 0000000E 00 .0. 4283829E-06 0. 0080008E 00 0. 0088800E 08-0. 9370365E-04 0. 0808080E 83 0. 0008008E 00 0. 4836362E-06 0. 8000088E 00 

0 . 0000000E 00 0. 0000008E 00 0. 0009000E 00 0. 0080000E 00 0. 0080000E 89-0. 1320789E-02 0. 8080008E 00 0. 0008000E 00 0 . 0000000E 00 

0. 4468843E-05 8. 8800008E 00 0. 2818039E-85 0. 1718219E-05 0. 000808OE 00 0. 8880008E 00-0. 1095916E-83 8. 0008880E 00 0. 0800888E 00 
0 . 0000000E 00 0. 4283829E-8S 0. 8000008E 00 0. 0088000E 00 0. 4836357E-06 0. 0O08008E 80 8. 0008800E 80-0. 9370269E-04 8. 0080008E 00 

0 . 0000000E 80 0. 0000080E 00 0. 0000008E 08 0. 0000808E 00 8. 0000000E 00 0. 0O00000E 00 0. 0000008E 00 0. 0080009E 00-0. 1320709E-02 

*** RONS 16 THRU 30 ARE ALL ZEROS *** 


B ARRAY 



BEND IX RESEARCH LABORATORIES 


***** SPACE BASE, CONFIGURATION U5B0DY> ***** 


CD 

MSN1 

KXSM1 

KXSMS 

i 

• 0. 5103602E 03 

a 5125000E 01 

0 2991782E 02 

CD 

MHM1 

KVSM1 

KVSWB 

2 

0. 5124204E ©3 

0. 4S16702E 01 

0. 3091702E 82 

CD 

MDM1 

KZSH1 

KZSWB 

3 

0. 1748300E 03 

0. 4S16702E 01 

0. 75667S2E 01 

CD 

, . MTfll 

KXSHfl 

KXTfil 

4 

0. 2020201E 02 

8. 2991702E 02 

0. 3091702E 01 

CD 

MSWR 

KVSHfl 

KVTOl 

5 

0. 2644302E 02 

0. 30.91782E 02 

0. 3625000E 01 

CD 

MSWB 

KZSHfl 

K2TA1 

6 

8. 2644982E 02 r 

0. 7566702E 01 

0. 3625000E 01 

CD 

' DTA1 

LDM1 

' LHM1 

7 

0. 8750090E 01 

0. 1088O00E 02 

a 2600000E 02 

CD 

LStli 

LTfll 

IX0M1 

8 

a 1100000E 02. 

8. 11000OOE 02 

, 0. 3972300E 84 

CD 

IXHM1 

IVDM1 

IVHM1 

9 

0. 1345900E 05 

0. 3686088E 04 

0. 3558500E 05 

CD 

I2DH1 

I2HM1 

MSWA 

10 

0. 3686000E 04 

a 3558500E 05 

0: 2600000E 02 

CD 

WSHB 

XDMi 

VDM1 

11‘ 

a 2600090E 02 

-0. 3141698E ©2 

0. 0000000E 00 

CD 

2DM1 

XSWA 

Y5WA 

12 

a000Q000E 00! 

0. 1650080E 02 

0. 0000000E 08 

CD 

ZSWH 

XHM1 

YHM1 

13 

a 595S301E 82 

-8. 1300080E 02 

,0. 0000000E 00 

CD 

ZHtll - 

XSMB 

VSWB 

14 

0. 0000000E 00 

0. 165800OE 02 

0. 0000000E 00 

CD 

ZSWB 

XSM1 

VSM1 

15 

-a 5959299E 82 

a 5500008E 81 

9; 0008880E 88 

CD 

* ZSW1 

XTA1 

YTA1 

16 

0. 00OQ000E 08 

0. 1658008E 02 

8. 00O0000E 00 

CD 

ZTfll 

KS12X - 

KS12V 

17 

0. 0000000E 00 

0. 1966830E 05 

0. 5998410E 07 

CD 

KS12Z 

CS12X 

CS12Y 

18 

0. 6557202E 03 

0. 36216O0E 82 

0. 2677732E 03 


MODEL PARAMETERS 



BEND IX RESEARCH’ LABORATORIES 



CD 

CS122 

KS13X 

KS13V 


19 

0. 3587900E 01 

0. 1966830E 05 

0. 5S98410E 07 


CD 

KS13Z 

CS13X 

CS13V 


20 

0. 6557202E 03 

0. 3621600E 02 

6. 2677732E 03 


CD 

CS13Z 

KS24X 

KS24V 


21 

0. 3587900E 01 

0. 1841400E 04 

8. 7824030E 06 


CD 

KS24Z 

CS24X 

CS24V' 


22 

0. 8954901E 02 

0. 1550Q08E 03 

0. 3148030E 04 


CD 

CS24Z 

KS35X 

KS35V 


23 

0. 32O3980E 01 

0. 1841400E 84 

0. 7024830E 06 


CD 

KS35Z 

CS35X 

CS35V 


24 

0. 8954901E 02 

0. 155O000E 03 

0. 3148030E 04 


CD 

, CS35Z 

XLSWA 

XLSWB 


25 

0. 3203900E 01 

0. 1040800E 03 

0. 104000OE 03 


CD 

RE IX 

RE1V 

RE1Z 


26 

8. 0000000E 00 

0. 0000000E 00 

0. 00O0000E 00 


CD 

RR1X 

RA1V 

RA1Z 


27 

0. 00O0000E 00 

0. 0000000E 00 

0. 0000000E 00 

It* 

CD 

RE2X 

RE2V ' 

RE2Z 

r 

28 

0, 0000000E 00 

8. 0000000E 00 

0. 000O0O0E 00 

H- 





o 

CD 

RR2X 

RA2V 

RA2Z 

■O o 

29 

0. 00O0000E 00 

0. 0008000E 00 

0. O000000E 80 

■*» § 

CD 

RE3X 

RE3V 

RE3Z 

“S Q 

30 

0. 0000000E 00 

0. 0000000E 00 

0. 0000000E 00 

O § 





O T&> 

CD - 

RA3X 

RA3V 

RA3Z 

a/ r- 

31 •' 

0. 0000O0OE 00 

0. 0800000E 00 

0. 0000O0OE 00 

jO "0 





c 

CD 

. RE4X . 

RE4Y 

RE4Z 

> o 

32 

0. 0000000E 00 

0. 000O080E 00 

0. 000000QE 00. 

c ^ 





■5 r- 

CD 

RA4X 

RA4V 

RR4Z 

> C/3 

33 

0. 0000O00E 00 

0. 8000000E 00 

0. 0000080E 00 


CD 

RE5X ’ 

RE5V . 

RE5Z 


34 

0. 0080000E 00 

0. 0000000E 00 

0. 000O080E 00 


CD 

0MG1X 

0MG1V 

0MG12 


35 - 

0. 0000000E 00 

0. 0000000E 00 

0. 0800000E 00 


CD 

0MG2X 

0MG2V 

0MG2Z 


36 

0. 0000Q08E .00 

0. 0O00000E 00 

0 0080000E 00 


CD - 

011G3X 

0MG3Y ’ 

0MG3Z 


3 7 

0. 000O00OE 00’ 

0. O000000E 00 ' 

0. 000O000E 00 



II- 
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CD 

OMG4X 


OMG4V 


OMQ4Z 


38 

0. 0O00000E 

00 

0. 0000000E 

00 

0. 00O0000E 

00 

CD 

0MG5X 


0MG5V 


0MG5Z 


29 

0. 0000000E 

00 

0. 0000000E 

00 

0. 0000O00E 

00 

CD 

DPHI24 


DTHT24 


DPS I 24 


40 

'0. 000000OE 

00 

0. 00O0000E 

00 

0. 000000OE 

00 

CD 

DPHI35 


DTHT35 


DPS I 35 


41 

0. 0000000E 

00 

0. 0000000E 

00 

0. 000OO00E 

00 

CD 

DPHI12 


DTHT12 


DPS I 12 


42 

0. 0000000E 

00 

0. 0000000E 

00 

0. 000O000E 

00 

CD 

DPHI13 


DTHT13 


DPS I 12 


43 

8, 0000000E 

00 

0. 0000000E 

00 

0. 0000000E 

00 

CD 

PHI1 


THT1 


PSU 


44 

0. 0000000E 

00 

0. 0000000E 

00 

0. 2000000E- 

-01 

CD 

RDT1X 


RDT1V 


RDT1Z 


45 

0. 0000000E 

00 

0, 0000000E 

00 

0. 0000000E 

00 

CD 

R1X 


R1V 


R1Z 


46 

0, 0000000E 

00 

0. 0000000E 

00 

0. 0000000E 

00 

CD 

DT 


IPRT 


TEND 


4 7 

0. 5000000E- 

■01 

0. 2000000 E 

02 

0. 5008000E 

02 



APPENDIX B 


TRANSFER FUNCTIONS FOR SYMMETRIC 
Z-AXIS WITH STATE VECTOR 

X (t) = ® 2 ±’ U 42 ’ ^21 ' Y 42| 


Space Construction Base, 
Configuration 1 



ftUM COPS 

NUMBER OF CONTROLLED STATES ? S 

NUMBER OF CONTROLS AND OUTPUTS ? 2,5 

LIST INPUT MATRICES » <N0> VES 

ENTER INPUT MATRICES DATA FILE. X2AXIS2 

LIST CALCULATED NUI1. & DENOM. PAIRS ? <NO> VES 

name of generated transfer function file ? *222 

INPUT MATRIX A ■ 

O.O0OOE+OO 0.0000E+0O ' O.0O 00E+00 • 0.O000E+00 

. O.OO00E+OO -0.4765E-02 0.4231E-02 -0.8SSSE+00 

O.O000E+00 O.4739E-02 T0.8463E-02 0.8SS0E+00 

O.8000E+00 0.1000E+0I 0.O0DOE+0& 0.0000E+00 

O.000OE+00 0.0O0CE+O0 0. 1000E+O1 0.0000E+00 

w INPUT MATRIX E * 

|L 0.3724E-05 0.OO0OE+O0 

-0.3724E-O5 0. 1321E-02 

0.00O0E+00 -0. 1321E-02 

Q.O00OE-00 0.0000E+00 

0.00O0E+00 0.0000E*00 

INPUT MATRIX C - 

0.I000E+01 0.0000E+00 0.0000E+00 0.0000E+00 

0.0OO0E+00 . 0. 1000E+01 0.0000E+00 0.0000E+00 

0.00O0E+00 0. 0000E+80 0. 1000E+01 0.0000E+00 

0.0000E+00 0.0000E+00 0.0000E+00 0.1000E+01 

0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 

*** j/xjuh xxtx**,xt**xx**xxt*xxxx*txxtxxxxtxxtxxx*x 


OUTPUT 1 


INPUT 1 


----- NUMERATOR COEFFICIENTS 

0.37244E-0S 0.00000 0.41150E-05 0.14901E-0? 

0.38339E-0G 

........... DENOMINATOR coefficients 

O.13228E-01 1.1050 0.43524E-02 0.10295 

0.00000 

txtx*tx*ttt*txzxxtttxtx*ttt*tt*xxxtx*ixxxxxxxtt* 

OUTPUT 1 INPUT 2 

........... NUMERATOR COEFFICIENTS 

0.0000E+00 0.00000 0.00000 0.00000 0.00000 

0.1183E+00 0.00000 

-0.2365E+00 DENOMINATOR COEFFICIENTS ■•»*«•««• 

0.0800E+00 0..13228E-01 1.1050 0.42524E-O2 0.10295 

0.0000E+00 0.00000 


tt**xxxtxxtx**txxt*xt*xtttxxx*xxxx*%xxxttx*xtxxx 

• OUTPUT 2 INPUT 1 

........... NUMERATOR COEFFICIENTS •*»*«• «««“ 

-0 .37244Er05-0. 298O2E-C7-0 . 88073E-C6 0 . 00000 

0.00000 

........... DENOMINATOR COEFFICIENTS 

0.I3228E-01 1.1050 0.425H4E-02 0.’l0295 

0.00000 

0.0000E+00 

0.000CE+00 tttxxxtxxxxtxttxxtxxxxxxxtxxtxxtxxxxtxxtxtxtxxxx 

0.0000E+00 

0.000OE+00 

0.1000E+01 





OUTPUT' a INPUT a 

........... NUMERATOR COEFFICIENTS 

Q.1321OE-02 0.55879E-G5 0.15614E-G3 0.00000 

0.00000 

• DENOMINATOR COEFFICIENTS ........... 

0. 132a3E-0i 1.1050 G.42524E-02 0.10295 

0.00000 : 

txxxxixxxxxtxxxxttxxxxxtxxtxxxxxxxttttxxxxxxxxxx 

OUTPUT 3 INPUT 1 

........... NUMERATOR COEFFICIENTS ?«•.*■•« = > 

0.00000 0.00000 -0.3225OE-05 0.00000 

0-00000 

........... DENOMINATOR coefficients ........... 

0.1322SE-01 1.1O50 0. 42S24E-02 0.10235 

0.00000 

dd zxxxxxxxxxtxxxfttxxxtxtxxxxxtxttxtxxxxxxxxtxxxxx 

i 

bo OUTPUT 3 INPUT 2 

........... NUMERATOR COEFFICIENTS ......... 

-0 . 1 321 0E-02-0 .59605E-07-0. 33O26E-0S 0 . 00000 

0.00000 

........... DENOMINATOR COEFFICIENTS 

0.13228E-01 1.1O50 0.42524E-0E 0.10295 

0.00000 

XXxXttttXXXXXXtXXXtXtXXXtXXXXXXiXXXXXXXXXXXXXXX* 

OUTPUT 4 INPUT 1 

NUMERATOR COEFFICIENTS ........... 

-0 .37253E-05-0 . 31490E-07-0 .S3073E-06 


0.00000 

0.00000 


DENOMINATOR COEFFICIENTS 


0. 13223E-01 1.1050 

0.00000 


0.42524E-02 0.10295 


xxxxxtxtxxxxxtxxtxtxxxxxxxxxxxxxtxtxtxxxxxxxxxxt 

OUTPUT 4 INPUT 2 

........... NUMERATOR COEFFICIENTS ........... 

-0.1E62GE-08 0.13210E-02 0.S58S4E-05 0.15614E-03 

0.00000 

........... DENOMINATOR COEFFICIENTS . — 

0. 1322SE-01 1.1050 0.42524E-02 0.10295 

0.00000 

xt*x***xt*xx*titxttxt*xxt*xtx.xtxtxxxxxx*xxxxx %** 

OUTPUT 5 INPUT 1 

........... NUMERATOR COEFFICIENTS ■»«—*»* 

0.00000, 0.00000 -0. 18626E-07-0'.3H250E-05 

0.00000 

........... DENOMINATOR COEFFICIENTS 

0. 13228E-01 1.1050 O.42524E-02 0.1029S 

0.00000 

xxxtttxtxtxxxxxtxxxxxxxxxxxxxxxxtxxxxxttxxtxttxx 

OUTPUT 5 INPUT 2 

........... NUMERATOR COEFFICIENTS 

0.00000 -0 . 13210E-02-0 . 35G23E-O7-0 . 33024E-CS 

0.00000 

........... DENOMINATOR COEFFICIENTS =»- 



APPENDIX C 


DERIVATION AND ANALYSIS OF MINIMUM (TRACE 

SOLUTION FOR STATE CONTROLLER GAIN MATRIX C 



2 

The problem, minimize |Cx.-y. [ with respect to x where 

x. is the jth column of G and y. is the jth column of 
J J 

1 

A p + T I * -*- s e Q u;i - va l en 't to the problem, minimize 
T -b 

(trace E c E c ) with respect to G where: 

E c A A c + | I 


I A identity matrix of dimension n 


Expanding 




we have 


T T 

X . C Cx - 2y Cx. 
J 1 11 




Taking the gradient with respect to x , setting equal to 

1 

zero for a critical value and solving for x-, we obtain 

3 

Xj - C* Zj (C-l) 

where 

C*4(C T C) -1 C T (C-2) 

Therefore: 

G = C*(A p + if I) (C-3) 

Substituting equation (C-3) into eauation (7-23) yields: 
x p = A c x p + Dw (7-22) 


C-l 



where:- 


= A p- CC 


(A 


P 


rp I ) 


(C-4 ) 


It was desired that all of the eigenvalues of the con- 
trolled plant matrix, A c , be - jL. The c linearly 

T 

independent columns of C, that is, the vectors' :C6\ : 

—l ’ 

i = 1 to c, 


where: 




6 . ) 

ic' 


(C-5) 


=the Kronecker delta 
T 

are ..eigenvalues of A c corresponding to the eigenvalue, 
-T/T,- since 




(C-6) 


Therefore, at least c of the eigenvalues of A c . are -1/T 

since the number of linearly independent .eigenvectors 
corresponding to a given eigenvalue is not greater than 
the multiplicity of the eigenvalue and the transpose of 
a matrix has the same eigenvalues as the matrix. 

The rank of the' sum of two matrices is equal to or less 
than the sum' of the 1 ranks of the matrices. Furthermore, 
the rank of a nonzero scalar times a matrix is equal to 
the rank of the matrix. Hence, we have: 

rank A c _< rank R + rank CC* 


C-2 



where: 


R = (I - CC*)A p (C-7) 

If and Mg are pxr and rxq matrices, respectively, and 

each is of rank r, then M^Mg is of rank r. Hence, the 
matrix CC* has rank c. If the rank of a square matrix 
of dimension n is r, then at least n-r of the eigenvalues 
of the matrix are zero. Furthermore, a necessary con- 
dition for linear control system asymptotic stability, 
that is, for no disturbance, x-K) as V x(0)> is that 

all of the eigenvalues of the controlled plant matrix 
have negative real parts. Therefore, a necessary condition 
for R is: 

rank R > n-c ' 


The rank of the product of two matrices can be nb greater . 
than the rank of either matrix. Thus, if r^ of the ; 
columns of A are linearly independent and r _< n-c, 

P <v 


then the rank of R is less than or equal to n-c. If the 
jth column of A is in the range of C, that is, . if 

ir 

A = Cc where A •- is the jth column of A and c is a c 
~Pj “ ~Pj P “ 

component constant vector, then the jth column of R.is 

zero. Thus, if r of the columns of A p are linearly 


independent and r > n-c, then r+c-n of the columns of R 
will be zero and n-r will be zero or linearly dependent, 
that is, the rank of R will be n-c. Therefore, for 

T x 

the minimum (tr S E-) 2 solution for G to be asymptotically 

c c 

stable, it is necessary that the rank of R be n-c, that is, 
that n-c columns of the plant matrix must be linearly 
independent from the columns of C. 


C-3 



APPENDIX D 


Algorithm for Generating Optimum State 
Controller or Estimator Gain Matrix 



1. Compute K = C*(A + I I) or 

P T 

20 t Tx-1 

- (A + y 1 ) B (Bp B p > 

2. Set e k = n ” 4 (X, 1, , 1) T 

3. Set X, «= 0 

k 

4. Go to subroutine. 

5 . Set e = © k+1 / |® k+1 | 

6 . Set X = X fe+1 

7. Compute 6 = 0.01 largest) k | 

8. Set i ='l 

9. Set j = 1 

10. Set -k. . = ,k. . + 5 

11 . .Set e^ = b 

12. Set X^ = X 

13. Go to subroutine. 


D-l 



14. 

Compute = X k+1 - X 


15. 

Set k i0 ■ k ij - 5 


16. 

If j = n(m), go to step 

19. 

17. 

Set j = j+1 


18. 

Go to step 10. 


19. 

If i = c (dimension of x) 

, go to step 22 

20. 

Set i = i+1 


21. 

Go to step 9. 


22. 

Compute A = largest J 6 . . | 


23. 

get K = K - | 


24. 

Set e^ = e 


25. 

Set X k = X 


26. 

Go to subroutine. 


27. 

I£ l x k +i -M / x k+i >0 - 01 ’ 

go to step 5. 

28. 

Set G or H = K 



D-2 



SUBROUTINE 


1. 

Compute .E = W c 

(A p -fCK + ~I-)W- C * 


■ W e 

20 

(A - KB + — I)W * 

2. 

Compute e^. + ^ = 

E T E e k 

3. 

Compute = 

1 — k+1 1 

4. 

If X k+1 ~ 0 ’ go 

to step 6. 

5. 

H i ^ k+1 -* k i 

/\^ +1 <0.01, return 

6. 

Set e k = e k+1 

/| -k+1 

7. 

Set X k - X k+1 


8. 

Go to step 2. 



or 


D-3 



APPENDIX E 


T 

Derivation and Analysis of Minimum (trace E g E e ) 
Solution for State Estimation Gain Matrix H 



T 

The problem, minimize- |B ~^.±\ wi_ th respect to 
x . where : 


x T = ith row of H 

—l 

. T . . . . A . 20 _ 

— ith row of A + — I , 

T i 

is equivalent to the problem, minimize (tr E E y 
with respect to H where: 

20 

E = A HB + ™-I 
e t 

I = identity matrix of dimension A 


Expanding 


- Zi 


yields : 


m rp rp rp rp 

4 BB i x i -ByVx. +£ 1± 


Taking the gradient with respect to x. , setting equal 

1 T 

to zero for a critical value and solving for we 
have 



Z^B T (BB T ) 1 ' 


(E-l ) 


E-l 



where : 


B = [ Bp | 0 j 0 ] (7-44) 

Therefore: 

H = (A + “ I ) B T (BpBp'V 1 (E-2) 

Substitution of (E-2) into equation (7-42) yields: 



where : 


(E-3) 

H 1 " (A c T 

T T -1 

B P (b pV • 

(E-4) 

Since 

(B B T ) -:L B ' 0 | 0 1 T , 
p p' p i — i — J 7 


x e =[/(0) 

(E-5) 


the estimates for x and x and thus the estimates 

—v — w 

for the unknown part of the control and plant disturb- 
ances are zero for all time. 


E-2 



APPENDIX F 


EQUATIONS FOR THE FULL BEAM 
AND SPLIT-BEAM 



F.l 

(i) 

(ii) 

(iii) 

(iv) 

(v) 

(vi) 


EQUATIONS FOR THE FULL BEAM AND SPLIT-BEAM 

Translational equations for the full beam 

“*1 + (M 1 )u 1 -K iV2 -Ku 3 -(K 0+ K Pi )B 1 -K 13i S 2 -K |8 B3 = P^ 

mii 2 + (K+K 1 )K 2 -K 1 p 1 -Ky 4 -(K 0 -K )B 2 +K 0 Bl -K 0 B 4 = ?' 

XX <4 

“ii - < 2 K+ K 1 )V 1 -K Bl _ 2 -K lPl+1 -Ky 1+2 -K 01 

+K 0 »l-2 6. +1 S 6 i+2 ’ \ 

i = 3,5,...., 2n-3 

mi i +(2K+K 1 )u 1 -K Pi _ 2 -K V±+2 

+K 0 h-2 + \ h-1 - K 0 »i + 2 = \ 

i = 4,6,...., 2n-2 

my 2n-l + ( K+K l) v 2n-l“ K l w 2n" Kv 2n-3"( K p 1 “ K p* B 2n-l“ K 0 1 8 2n +K p B 2n-3 


= F 


e 


2n-l 


my 2n +(K+K l )u 2n K l u 2n-l“ Ky 2n-2 +(K j3 1 +K 0 ) B 2n +K 0 1 B 2n-l +K 0 B 2n-2 


= F, 


'2n 


where 


F = F 
e i e i 


m 


'Xi • 


„ F + m r . co 
M T e l o 


F-l 



Rotational equations for the full beam 

i 

(i) + K a e l +K D 1 B 2 +K b R 3 + ( V K P 1 )Wl_K S l i‘ 2 "V3 = 

V 

Cii) J Eg + K a e 2 +K b 1 B 1 +K b e 4 + ( V K 0 1 ), ‘2 +K 0l Vi "V4 = T e 2 

(ill) J6 1 +K d B 1 +K b 6 i . 2 +K bi S 1+1 +K b 6 1+2 +K 0i ll 1 +K 19 iJ 1 _ 2 

_K fl :L 1J i+l _K 0 u i+2 ~ T e ± 

i = 3,5, , 2n-3 

(iv) J & j. +K d B ± +K b e t _ 2 +K b 1 +K b 6 i+2~ K 0 ^ x +K 0 y i-2 

+K 0 1 u i-l“ K {5 y i+2 = T e j _ 

i = 4,6, , 2n-2 


<v) J »2n-l +K a 6 2n-l +K b 1 B 2n +K b e 2.i-3' F(K 0 1 - K !S ),i 2n-l- K ? 1 ' 1 2n +K 


= T. 


'2n-l 


( y i) J ^2n^ K a B 2n +K b 1 3 2n-l +K b e 2n-2 ' (K 0 1 +K 0 )y 2n +K !3 1 U 2ii-l +K 0 y 


= T 


'2n 


where K a = ^K^+Kj+K^ 

K d " VV K a- 
K b " V K e 

K bl = 

T e 1 * F »i ' J “o 


0 v 2n-3 


2n-2 


F-2 



p.3 Translational equations for the bodies involved in the 

connection of half-beams 

( i) m tf n _i +m Sv-i e n-l~°-‘ 5Ky n+3"° ’ 5K l u n+2 +0 * 5(K+K l )TI n+l 

-0 . 5X^+0 . SCK+K-l) y n _!-0 . 5Ku n _ g-0 . 5 ( K p +K 0i) S n+1 
-O..5K (3i B n -O.5(K 0i -K 0 )B tl _ 1 +O.5K 0 3 n _ 3 -O.5K 0 B I1+ 3-O.5K 0i 3 n+2 

t 

= F 

e n-l 

( ii ) my n +m£ n B n -0 . 5Ku n+4 +0 . 5 (K+K 1 ) V n+ 2~° • 5K l u n +l 

+0. 5(K+K 1 )u n -.0 . 5K 1 y n _ 1 ~0 . 5Ku n _ 2 +0 . SK^ 8 n+1 

+0 • 5 < 1 3 n + ° • 5 Vn-l + ° ' 5 V n-2 

-0 . 5KpB n+4 ~0. 5(K^-Kp^)B n+2 ” F e n 

( iii ) m U n+1 + me n+1 e n+1 -0 • 5K ^ n+ 3"° • 5K l v n+2 +0 ' 5 ( K+K l )u n+l"° * 5K l y n 

+O.5(K+K 1 )u n _ 1 -O.5Ku n _ 3 -O.5(K s5 +K 0i )B n+1 -O.5K^ i B n 
+0 • 5 ( K r K ^ } 6 n-l +0 ' 5K ^n-3-° • 5K ^ V2"° • 5K ^n + 3 = 


(iv) mU n+2 +me n+2 B n+2 -0.5Kq n+4 +0 . 5(K+K 1 )y n+2 


where : 


+0 . 5(K+K 1 )y n -0, 5K 1 u n _ 1 -0 . 5K]i n _ 2 -0 ■ 5(K 0 -K^) 3 n+2 

+0 * 5 Vn+l + ° * 5 ( W ® n+ ? • 5K 0! B n- 1 

+0 . 5K^ B n _ 2 - 0 • 5K p ^ n +4 


F* = *(F +F )- ~ 
e i e i e j Mr ] 


= F e 

n+2 

F +m(„ r . ~e • )ro 


% e 


% i 


oJl 


1=1 for i=n-l,n 
%=2 for i=n+l,n+2 
j=i+(-l;* +l 2 


F-3 



F.4 Rotational "equations for the bodies involved in the 

connection of half-beams 

(i) _K 0 (,, n+3-' 1 n+l )_d i K l (v n+2-> I n+l ) 

-d 2 K lC>V l W ) " K ? (w n-r l, n-3 ) 

+ <V°- 5K 0 )f! n+3 +( % +0 - 5K 9 + % 11 +Q • “s^n+l 

+c % 12 -°- 6K e 1 )|3 n +( %' 0 - 5K e )8 n-3 

= T" , (J<)% 
n-1 1 


(ii) 

J 3 n “ J 3 n+2 J w °2 J 

Ui 

°1 

(iii) 

J 3 n+l " J 3 n-1 = J W °1 r 

• 

' ^ w 02 


* *" 

* _ , 

(iv) 

C J + J 2 ) P n +2'" K p^n+4" y n+2^ 

~dg 


- d 4 K l - K p Cw n -V n _ 2 > 

+ ( V° • 5K s 5 e n + 4 +< V°- SK e +K ^ i3 +0 • 5K 9i > e n+2 

t ( K V , 1 3r 0 - 5K e 1 )6 n + l +( V 0;5 V K 4 - 14 +0 - 8K 9 1 )6 n 

+ (\ 4 '°- S \~> 8 n-l +( V° • 5K 8 )B .-2 


= \ +a - (J+J 2> “°2 


F-4 



where: 


J 1 = 


* 

J 2 = 


m d 


a. 


13 e n-l 


d. - -d. 


Oj 

i (d 


12 


a. 

- d 


31- 


^ 'Vi 'Vi 

d 2 = _d 3 != * ^ d 12 " d 13^ 


K. 


If 


11 


«* -d 1 K 1 


S 12 = % 


14 


K 


*12 = -d 2 


K h 


12 % 13 
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APPENDIX G 


MODE SHAPES OF THE 
NATURAL FREQUENCIES 
FOR THE ORIGINAL BEAM 



§5 mm 


NATURAL FRfJULNCY. .5821781*02 


TX 

) X 

TY 

TY . 

T 2 

U 

RX 

rx 

RY* ' 

RY 

R 2 ' " 

■r l 

real 

1 HAG 

REAL 

1 HAG 

real ' 

[ HAG 

REAL 

I HAG 

real 

1 HAG 

.. REAL 

I HAG 

-.3 3024E * 03 

0. 00000 

.23506E-12 

0.00000 

.448921-12 

0.00000 

.1/9001-13 

O.OOOOO, 

" - .19 128 E-12 

6. oo5oo " 

9 7 4*2 1 1- 01 

plqobpq“ 

~.33024C*00 

0.00000 

• 1050 7E - 12 

0.00000 

_ .430091-12 

0.00000 

-.12107E-14 

0.00000 

-. 19146E.-12 _ 

0.00 0 00” —.974211 — 01” 

o_._6oopq” 

-.1 36 12E *00 

0.00000 

. 10979b- 12 

' 0.00 000 

.41581E-12 

0.00000 

-.254461-13 

0.00000 

— .14 710E-12 

0.00000 ' 

- » Ooc -oi 

’o.oooob* * 

-.136l2E*00 

0.00000 

* 2 1 99Gt — 12 

O.OOOOO 

.425161-12 

0.00000 

116761-13 

0.00000 

-. 16659E— 12 

o. ooooo 

— ,. 92 8Q6E"-0*1 

'o'.poQoa* * 

.40347E-01 

0.00000 

.2234 76-12 

0.00000 

.426051-12 

0.00000 

-.231421-13 

0.00000 

— • 15103E— 12 

0.00300* 

'-.7aUi-)L-oi~ 

0. 00 9 00. * 

*.40347c-Gl 

o.ooooo 

.222196-12 * 

0.00000 
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RZ 
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RZ 
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ri5797E*00 
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-.1058AE-13 

a. ooooo 
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0.00000 " 

'7.312666-13.7 

76^00000 

-.9A9571-1A 
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0.00000 
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.A3371E-1A 

0.00000 
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. 12396E *00 
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. o.ooooq.. 

"7-9229bE — 13 

_o'. poo op 

77.2508p7.9T' 

“0". ooooo 

« 1322 bE *0 3 
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*0.000007 

7-.2717UE-1 3 
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7.. 732256-13 7 

0.0007007 
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0.00000 

* *‘*.539858-01 
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•18062E-13 
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•85277b- 12 
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;i9870E*00 

0.00000 

— . 8 2q58 E + 00 
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" '*115316—01 
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jp .'poo op: 
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6 : 00006 ’ 
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aToOOOO 
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o.oo'oob" 

— .69730E — 12" 

' O.OOOOO: 

•2165>Oc“l2. 

''b.'obooo 
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0.00000 
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rr 

TZ 
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RX 

RY 
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IMAG* 
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ImaG 
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IMAG 

REAL 

1 RAO 

REAL 

IMAG 

REAL 

IMAG, 

-.225056-13 

0 .000 00 

. 155136-12 

0.00000 

. 199926*00 

0.00000 

— . 52 dlOt *00 

0.00000 . 

-. 172 506-12 ,, 

0,00000 

,'.1508307,12 

'o.oooou. 

-.296756-1 3 

0.00000 

— .1305 8 L - 1 2 

0.00000 

, -» 199921*00 

0.00000 

• 528 LOE * 00 

0.00000 

".809516-13 

0.000007 

" r . 2 85 786 71 3 

'~a. ooooo' 

‘.56251E-1 3 

0.00000 

-.957081-13 

o.oooou 

-.1 L209E*00 

O.OOOOO 

— .105 70c* 00 

0.00000 

•23665E-12 

0.00000 

”.215581-13 

"o.ooooo] 

.. '. 16319E-13 

0.00000 

.910001-13 

0. 000 00 

.112091*00 

0.00000 
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-.209756-12 
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0.00000 
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o.oooou 
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0.00000 
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0.00000 

. 12053E— 12 

0.00000 
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0 . OuOOO 
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o.uoooo 
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o.ooooo” 
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0 . ooooo* 

* .379116-13, 
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-.53290L-13 

JJ.OOOOO 
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0.00000 

.323591-01 

0.00000 _ 
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o._ooo_oo__ 

7 .2939511^13' 
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“ .27950E-U 

o.ooooo 
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_"V. ooooo* 
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0.00000 

-. 323596-01 

0.00050, 

.*2 78 026 - 12 "" 

' 0.000007 

, 7 . 1856167.13 
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77.555256-13 

o.ooobo 

_ U5767E-12* 

o.oooou 
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0.00000 
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'o.ooooo 
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0 . 00000*7 
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-.15108E-12, 

0*00000 
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0.00000 
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_o. ooooo _ 

.15177E-12 

OjOOUOO '7 

150966-1' 3 

7 \."oddoq] 

7 -33072E-13 

o.ooobo 

" ^ 153231— 12 

' 0.00000 

.203251-02. 

0.00000 

.59o23E-02 

0.00000 

.22537E-12' 

, 0.000007 

4 7toc — 1 3 

77 P ." 06000 " 

’"-.196 70E-13 

o*qqqoo 

-.165 1 3c— 12 

0*00000 

“ ‘ .203256-02" 

o.ooooo 

-.5962 36-02 

0.00000 

-.27592E-L2*, 

0 . odyoo" 

7.815 5 31-13 

0 .* 00' 00 / 

.1076 71-12 

0.000 00 

-.5571 11-13 

. 0.00000 

.180351.-01 . 

0.000 00 

.323596-01 

0.00000 

-.278351-13., 

0.000 00 „ 

Id 15 71-13 

. 0.00000 

l^.tOAloE -12 

0.00000 

.27595L — 13 

. 0.00000 

I -.180351—01 . 

o.ooooo 

-.323591-01 

0.00000 

-.laiobcria 

0 . oaooo_ 

Y*‘l3669fc~LZ 

y o', oop 06 ” 

’ .133251-12 

O.OOJOO 

— . 192 1 31 — 12 

0.00003 

.8 629 31—0 1 

0.00000 

• 350971—01 

0.00000 

-.5 56591-12 

O.OOoOO 

-.81231E-T3 

0.00000 

_ -.95133E-13 

0.00000 

• 19581e-12 

O.UOQOJ 

-.862531-01 

0.00000 

-.350971-01 

Q.OCOOO 

-.896351-13 
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.1 55.1 7c — 12 

_■ 0.00 ood] 

7’_.52528L-12 

0.00000 ] 

-.118566-12 

... 0.00000 

112091* 00.. 

o.ooooo 

. -..10570b *00 

o.ooooo 

-.512181-12-' 

o.ooooo'. 

7*155531 -12" 

■ OaOSpOO” 

"-.37959E-12 

0.00000 

.10668L-12 

0.00000 

-.112091*00 

0.00000 

. IO570E *00 

0.00000 

-• 13197E-12 

o.ooooo' 

.171196-12 

o.ooooo* 

• 3739;,t-12 

0.03000 

. 13 d 131- 12 

o.ooooo 

199921*00 
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